Cross-species protein sequence and gene structure prediction with fine-tuned Webscipio 2.0 and Scipio

Exploring the genetic landscape of nitrogen uptake in durum wheat: genome-wide characterization and expression profiling of NPF and NRT2 gene families

Nitrate uptake by plants primarily relies on two gene families: Nitrate transporter 1/peptide transporter (NPF) and Nitrate transporter 2 (NRT2). Here, we extensively characterized the NPF and NRT2 families in the durum wheat genome, revealing 211 NPF and 20 NRT2 genes. The two families share many Cis Regulatory Elements (CREs) and Transcription Factor binding sites, highlighting a partially overlapping regulatory system and suggesting a coordinated response for nitrate transport and utilization. Analyzing RNA-seq data from 9 tissues and 20 cultivars, we explored expression profiles and co-expression relationships of both gene families. We observed a strong correlation between nucleotide variation and gene expression within the NRT2 gene family, implicating a shared selection mechanism operating on both coding and regulatory regions. Furthermore, NPF genes showed highly tissue-specific expression profiles, while NRT2s were mainly divided in two co-expression modules, one expressed in roots (NAR2/NRT3 dependent) and the other induced in anthers and/ovaries during maturation. Our evidences confirmed that the majority of these genes were retained after small-scale duplication events, suggesting a neo- or sub-functionalization of many NPFs and NRT2s. Altogether, these findings indicate that the expansion of these gene families in durum wheat could provide valuable genetic variability useful to identify NUE-related and candidate genes for future breeding programs in the context of low-impact and sustainable agriculture.

Toll like receptors and their evolution in the lymnaeid freshwater snail species Radix auricularia and Lymnaea stagnalis, key intermediate hosts for zoonotic trematodes

One of the major evolutionarily conserved pathways in innate immunity of invertebrates is the toll-like receptor (TLR) pathway. However, little is known of the TLR protein family in gastropod molluscs despite their role in the transmission of human diseases, especially the common lymnaeid freshwater snail species Radix auricularia and Lymnaea stagnalis, key intermediate hosts of zoonotic trematodes. Using comparative genomics and gene prediction approaches utilising the freshwater snail Biomphalaria glabrata genome as a reference ten putative TLR proteins were identified in both R. auricularia and L. stagnalis. Phylogenetic analyses revealed that unlike other molluscs the lymnaeid species also possessed class 1 TLRs, previously thought to be unique to B. glabrata. Gene duplication events were also seen across the TLR classes in the lymnaeids with several of the genes appearing to exist as potential tandem elements in R. auricularia. Each predicted TLR was shown to possess the typical the leucine-rich repeat extracellular and TIR intracellular domains and both single cysteine clusters and multiple cysteine clusters TLRs were identified in both lymnaeid species. Principle component analyses of 3D models of the predicted TLRs showed that class 1 and 5 proteins did not cluster based on similarity of structure, suggested to be potential adaptation to a range of pathogens. This study provides the first detailed account of TLRs in lymnaeids and affords a platform for further research into the role of these proteins into susceptibility and compatibility of these snails with trematodes and their role in transmission.

Evolutionary Origin of Insulin-Degrading Enzyme and Its Subcellular Localization and Secretion Mechanism: A Study in Microglial Cells

The insulin-degrading enzyme (IDE) is a zinc-dependent metalloendopeptidase that belongs to the M16A metalloprotease family. IDE is markedly expressed in the brain, where it is particularly relevant due to its in vitro amyloid beta (Aβ)-degrading activity. The subcellular localization of IDE, a paramount aspect to understand how this enzyme can perform its proteolytic functions in vivo, remains highly controversial. In this work, we addressed IDE subcellular localization from an evolutionary perspective. Phylogenetic analyses based on protein sequence and gene and protein structure were performed. An in silico analysis of IDE signal peptide suggests an evolutionary shift in IDE exportation at the prokaryote/eukaryote divide. Subcellular localization experiments in microglia revealed that IDE is mostly cytosolic. Furthermore, IDE associates to membranes by their cytoplasmatic side and further partitions between raft and non-raft domains. When stimulated, microglia change into a secretory active state, produces numerous multivesicular bodies and IDE associates with their membranes. The subsequent inward budding of such membranes internalizes IDE in intraluminal vesicles, which later allows IDE to be exported outside the cells in small extracellular vesicles. We further demonstrate that such an IDE exportation mechanism is regulated by stimuli relevant for microglia in physiological conditions and upon aging and neurodegeneration.

De Novo Gene Birth, Horizontal Gene Transfer, and Gene Duplication as Sources of New Gene Families Associated with the Origin of Symbiosis in Amanita

By introducing novel capacities and functions, new genes and gene families may play a crucial role in ecological transitions. Mechanisms generating new gene families include de novo gene birth, horizontal gene transfer, and neofunctionalization following a duplication event. The ectomycorrhizal (ECM) symbiosis is a ubiquitous mutualism and the association has evolved repeatedly and independently many times among the fungi, but the evolutionary dynamics enabling its emergence remain elusive. We developed a phylogenetic workflow to first understand if gene families unique to ECM Amanita fungi and absent from closely related asymbiotic species are functionally relevant to the symbiosis, and then to systematically infer their origins. We identified 109 gene families unique to ECM Amanita species. Genes belonging to unique gene families are under strong purifying selection and are upregulated during symbiosis, compared with genes of conserved or orphan gene families. The origins of seven of the unique gene families are strongly supported as either de novo gene birth (two gene families), horizontal gene transfer (four), or gene duplication (one). An additional 34 families appear new because of their selective retention within symbiotic species. Among the 109 unique gene families, the most upregulated gene in symbiotic cultures encodes a 1-aminocyclopropane-1-carboxylate deaminase, an enzyme capable of downregulating the synthesis of the plant hormone ethylene, a common negative regulator of plant-microbial mutualisms.

The evolving genetic landscape of congenital disorders of glycosylation

Congenital Disorders of Glycosylation (CDG) are an expanding and complex group of rare genetic disorders caused by defects in the glycosylation of proteins and lipids. The genetic spectrum of CDG is extremely broad with mutations in over 140 genes leading to a wide variety of symptoms ranging from mild to severe and life-threatening. There has been an expansion in the genetic complexity of CDG in recent years. More specifically several examples of alternate phenotypes in recessive forms of CDG and new types of CDG following an autosomal dominant inheritance pattern have been identified. In addition, novel genetic mechanisms such as expansion repeats have been reported and several already known disorders have been classified as CDG as their pathophysiology was better elucidated. Furthermore, we consider the future and outlook of CDG genetics, with a focus on exploration of the non-coding genome using whole genome sequencing, RNA-seq and multi-omics technology.

Diversity and evolution of the P450 family in arthropods

The P450 family (CYP genes) of arthropods encodes diverse enzymes involved in the metabolism of foreign compounds and in essential endocrine or ecophysiological functions. The P450 sequences (CYPome) from 40 arthropod species were manually curated, including 31 complete CYPomes, and a maximum likelihood phylogeny of nearly 3000 sequences is presented. Arthropod CYPomes are assembled from members of six CYP clans of variable size, the CYP2, CYP3, CYP4 and mitochondrial clans, as well as the CYP20 and CYP16 clans that are not found in Neoptera. CYPome sizes vary from two dozen genes in some parasitic species to over 200 in species as diverse as collembolans or ticks. CYPomes are comprised of few CYP families with many genes and many CYP families with few genes, and this distribution is the result of dynamic birth and death processes. Lineage-specific expansions or blooms are found throughout the phylogeny and often result in genomic clusters that appear to form a reservoir of catalytic diversity maintained as heritable units. Among the many P450s with physiological functions, six CYP families are involved in ecdysteroid metabolism. However, five so-called Halloween genes are not universally represented and do not constitute the unique pathway of ecdysteroid biosynthesis. The diversity of arthropod CYPomes has only partially been uncovered to date and many P450s with physiological functions regulating the synthesis and degradation of endogenous signal molecules (including ecdysteroids) and semiochemicals (including pheromones and defense chemicals) remain to be discovered. Sequence diversity of arthropod P450s is extreme, and P450 sequences lacking the universally conserved Cys ligand to the heme have evolved several times. A better understanding of P450 evolution is needed to discern the relative contributions of stochastic processes and adaptive processes in shaping the size and diversity of CYPomes.

The visual pigment xenopsin is widespread in protostome eyes and impacts the view on eye evolution

Photoreceptor cells in the eyes of Bilateria are often classified into microvillar cells with rhabdomeric opsin and ciliary cells with ciliary opsin, each type having specialized molecular components and physiology. First data on the recently discovered xenopsin point towards a more complex situation in protostomes. In this study, we provide clear evidence that xenopsin enters cilia in the eye of the larval bryozoan Tricellaria inopinata and triggers phototaxis. As reported from a mollusc, we find xenopsin coexpressed with rhabdomeric-opsin in eye photoreceptor cells bearing both microvilli and cilia in larva of the annelid Malacoceros fuliginosus. This is the first organism known to have both xenopsin and ciliary opsin, showing that these opsins are not necessarily mutually exclusive. Compiling existing data, we propose that xenopsin may play an important role in many protostome eyes and provides new insights into the function, evolution, and possible plasticity of animal eye photoreceptor cells.

Efflux pumps as an additional source of resistance to trichothecenes in Fusarium proliferatum and Fusarium oxysporum isolates

Role of efflux-mediated toxin resistance to trichothecenes is known in trichothecene-producing species. However, the role of trichothecene efflux pump homologues in non-producing fusaria such as F. oxysporum and F. proliferatum was not investigated in detail. Analysis of the homologues of trichothecene efflux pump from multiple fungal species allowed us to uncover and catalogue functional gene copies of conserved structure. Putative Tri12 candidates in Fusarium oxysporum and F. proliferatum were characterised via expression profiling in response to different trigger compounds, providing supporting evidence for role of Tri12 homologues in the resistance to trichothecenes. Our analysis of Tri12 phylogeny also suggests that efflux-mediated trichothecene resistance is likely to predate the divergence of Trichoderma and Fusarium species. On the regulatory level, we posit that the increased tolerance of trichothecenes by F. oxysporum is possibly related to the decoupling of Tri12 homologue expression from pH, due to the deletion of PACC/RIM101 transcription factor binding site in its promoter region.

Transcriptome-derived investigation of biosynthesis of quinolizidine alkaloids in narrow-leafed lupin (Lupinus angustifolius L.) highlights candidate genes linked to iucundus locus

Unravelling the biosynthetic pathway of quinolizidine alkaloids (QAs), regarded as antinutritional compounds of narrow-leafed lupin (NLL) seeds, is fundamental to best exploit NLL as food or feed. We investigated 12 candidate genes connected to QA biosynthesis, selecting them by transcriptomic and genomic approaches, from the landscape of genes differentially expressed in leaves of the high- and low-alkaloid NLL accessions. Linkage analysis enabled the assessment of the location of the candidate genes in relation to iucundus, a major locus of unknown identity, that confers reduced QA content in seeds. The key finding was the identification of APETALA2/ethylene response transcription factor, RAP2-7, cosegregating with the iucundus locus and located within a region with highly significant QTLs that affect QA composition. We additionally identified a 4-hydroxy-tetrahydrodipicolinate synthase (DHDPS) gene involved in L-lysine biosynthesis as being closely linked to iucundus. The distributed location of other remaining candidates (including previously known QA genes) across different linkage groups, also indirectly supports the transcription factor as a possible regulator of lupin alkaloid biosynthesis. Our findings provide crucial insight into QA biosynthesis in NLL. Additionally, we evaluated and selected appropriate reference genes for qRT-PCRs to analyse the expression levels of QA genes in NLL.

Molecular evolution and functional characterisation of insulin related peptides in molluscs: Contributions of Crassostrea gigas genomic and transcriptomic-wide screening

Insulin Related Peptides (IRPs) belong to the insulin superfamily and possess a typical structure with two chains, B and A, linked by disulphide bonds. As the sequence conservation is usually low between members, IRPs are classified according to the number and position of their disulphide bonds. In molluscan species, the first IRPs identified, named Molluscan Insulin-related Peptides (MIPs), exhibit four disulphide bonds. The genomic and transcriptomic data screening in the Pacific oyster Crassostrea gigas (Mollusc, Bivalvia) allowed us to identify six IRP sequences belonging to three structural groups. Cg-MIP1 to 4 have the typical structure of MIPs with four disulphide bonds. Cg-ILP has three disulphide bonds like vertebrate Insulin-Like Peptides (ILPs). The last one, Cg-MILP7 has a significant homology with Drosophila ILP7 (DILP7) associated with two additional cysteines allowing the formation of a fourth disulphide bond. The phylogenetic analysis points out that ILPs may be the most ancestral form. Moreover, it appears that ILP7 orthologs are probably anterior to lophotrochozoa and ecdysozoa segregation. In order to investigate the diversity of physiological functions of the oyster IRPs, we combine in silico expression data, qPCR measurements and in situ hybridization. The Cg-ilp transcript, mainly detected in the digestive gland and in the gonadal area, is potentially involved in the control of digestion and gametogenesis. The expression of Cg-mip4 is mainly associated with the larval development. The Cg-mip transcript shared by the Cg-MIP1, 2 and 3, is mainly expressed in visceral ganglia but its expression was also observed in the gonads of mature males. This pattern suggested the key roles of IRPs in the control of sexual reproduction in molluscan species.

Unicellular ancestry and mechanisms of diversification of Goodpasture antigen-binding protein

The emergence of the basement membrane (BM), a specialized form of extracellular matrix, was essential in the unicellular transition to multicellularity. However, the mechanism is unknown. Goodpasture antigen-binding protein (GPBP), a BM protein, was uniquely poised to play diverse roles in this transition owing to its multiple isoforms (GPBP-1, -2, and -3) with varied intracellular and extracellular functions (ceramide trafficker and protein kinase). We sought to determine the evolutionary origin of GPBP isoforms. Our findings reveal the presence of GPBP in unicellular protists, with GPBP-2 as the most ancient isoform. In vertebrates, GPBP-1 assumed extracellular function that is further enhanced by membrane-bound GPBP-3 in mammalians, whereas GPBP-2 retained intracellular function. Moreover, GPBP-2 possesses a dual intracellular/extracellular function in cnidarians, an early nonbilaterian group. We conclude that GPBP functioning both inside and outside the cell was of fundamental importance for the evolutionary transition to animal multicellularity and tissue evolution.

In silico characterization and expression profiling of the diacylglycerol acyltransferase gene family (DGAT1, DGAT2, DGAT3 and WS/DGAT) from oil palm, Elaeis guineensis

The diacylglycerol acyltransferases (DGAT) (diacylglycerol:acyl-CoA acyltransferase, EC 2.3.1.20) are a key group of enzymes that catalyse the final and usually the most important rate-limiting step of triacylglycerol biosynthesis in plants and other organisms. Genes encoding four distinct functional families of DGAT enzymes have been characterised in the genome of the African oil palm, Elaeis guineensis. The contrasting features of the various isoforms within the four families of DGAT genes, namely DGAT1, DGAT2, DGAT3 and WS/DGAT are presented both in the oil palm itself and, for comparative purposes, in 12 other oil crop or model/related plants, namely Arabidopsis thaliana, Brachypodium distachyon, Brassica napus, Elaeis oleifera, Glycine max, Gossypium hirsutum, Helianthus annuus, Musa acuminata, Oryza sativa, Phoenix dactylifera, Sorghum bicolor, and Zea mays. The oil palm genome contains respectively three, two, two and two distinctly expressed functional copies of the DGAT1, DGAT2, DGAT3 and WS/DGAT genes. Phylogenetic analyses of the four DGAT families showed that the E. guineensis genes tend to cluster with sequences from P. dactylifera and M. acuminata rather than with other members of the Commelinid monocots group, such as the Poales which include the major cereal crops such as rice and maize. Comparison of the predicted DGAT protein sequences with other animal and plant DGATs was consistent with the E. guineensis DGAT1 being ER located with its active site facing the lumen while DGAT2, although also ER located, had a predicted cytosol-facing active site. In contrast, DGAT3 and some (but not all) WS/DGAT in E. guineensis are predicted to be soluble, cytosolic enzymes. Evaluation of E. guineensis DGAT gene expression in different tissues and developmental stages suggests that the four DGAT groups have distinctive physiological roles and are particularly prominent in developmental processes relating to reproduction, such as flowering, and in fruit/seed formation especially in the mesocarp and endosperm tissues.

An analysis of variability in genome organisation of intracellular calcium release channels across insect orders

Using publicly available genomic data, combined with RT-PCR validation, we explore structural genomic variation for two major ion channels across insect classes. We have manually curated ryanodine receptor (RyR) and inositol 1,4,5-trisphosphate receptor (IP₃R) ORFs and their corresponding genomic structures from 26 different insects covering major insect orders. We found that, despite high protein identity for both RyRs (>75%) and IP₃Rs (~67%), the overall complexity of the gene structure varies greatly between different insect orders with the simplest genes (fewest introns) found in Diptera and the most complex in Lepidoptera. Analysis of intron conservation patterns indicated that the majority of conserved introns are found close to the 5' end of the channels and in RyR around the highly conserved mutually exclusive splice site. Of the two channels the IP₃Rs appear to have a less well conserved organisation with a greater overall number of unique introns seen between insect orders. We experimentally validated two of the manually curated ORFs for IP₃Rs and confirmed an atypical (3799aa) IP₃R receptor in Myzus persicae, which is approximately 1000 amino acids larger than previously reported for IP₃Rs.

Evolutionary and genomic analysis of the caleosin/peroxygenase (CLO/PXG) gene/protein families in the Viridiplantae

Bioinformatics analyses of caleosin/peroxygenases (CLO/PXG) demonstrated that these genes are present in the vast majority of Viridiplantae taxa for which sequence data are available. Functionally active CLO/PXG proteins with roles in abiotic stress tolerance and lipid droplet storage are present in some Trebouxiophycean and Chlorophycean green algae but are absent from the small number of sequenced Prasinophyceaen genomes. CLO/PXG-like genes are expressed during dehydration stress in Charophyte algae, a sister clade of the land plants (Embryophyta). CLO/PXG-like sequences are also present in all of the >300 sequenced Embryophyte genomes, where some species contain as many as 10–12 genes that have arisen via selective gene duplication. Angiosperm genomes harbour at least one copy each of two distinct CLO/PX isoforms, termed H (high) and L (low), where H-forms contain an additional C-terminal motif of about 30–50 residues that is absent from L-forms. In contrast, species in other Viridiplantae taxa, including green algae, non-vascular plants, ferns and gymnosperms, contain only one (or occasionally both) of these isoforms per genome. Transcriptome and biochemical data show that CLO/PXG-like genes have complex patterns of developmental and tissue-specific expression. CLO/PXG proteins can associate with cytosolic lipid droplets and/or bilayer membranes. Many of the analysed isoforms also have peroxygenase activity and are involved in oxylipin metabolism. The distribution of CLO/PXG-like genes is consistent with an origin >1 billion years ago in at least two of the earliest diverging groups of the Viridiplantae, namely the Chlorophyta and the Streptophyta, after the Viridiplantae had already diverged from other Archaeplastidal groups such as the Rhodophyta and Glaucophyta. While algal CLO/PXGs have roles in lipid packaging and stress responses, the Embryophyte proteins have a much wider spectrum of roles and may have been instrumental in the colonisation of terrestrial habitats and the subsequent diversification as the major land flora.

Co-expression of xenopsin and rhabdomeric opsin in photoreceptors bearing microvilli and cilia

Ciliary and rhabdomeric opsins are employed by different kinds of photoreceptor cells, such as ciliary vertebrate rods and cones or protostome microvillar eye photoreceptors, that have specialized structures and molecular physiologies. We report unprecedented cellular co-expression of rhabdomeric opsin and a visual pigment of the recently described xenopsins in larval eyes of a mollusk. The photoreceptors bear both microvilli and cilia and express proteins that are orthologous to transporters in microvillar and ciliary opsin trafficking. Highly conserved but distinct gene structures suggest that xenopsins and ciliary opsins are of independent origin, irrespective of their mutually exclusive distribution in animals. Furthermore, we propose that frequent opsin gene loss had a large influence on the evolution, organization and function of brain and eye photoreceptor cells in bilaterian animals. The presence of xenopsin in eyes of even different design might be due to a common origin and initial employment of this protein in a highly plastic photoreceptor cell type of mixed microvillar/ciliary organization.

Animal eyes have photoreceptor cells that contain light-sensitive molecules called opsins. Although all animal photoreceptor cells of this kind share a common origin, the cells found in different organisms can differ considerably. The photoreceptor cells in flies, squids and other invertebrates store a type of opsin called r-opsin in thin projections on the surface known as microvilli. On the other hand, the visual photoreceptor cells in human and other vertebrate eyes transport another type of opsin (known as c-opsin) into more prominent extensions called cilia.

It has been suggested that the fly and vertebrate photoreceptor cells represent clearly distinct evolutionary lineages of cells, which diverged early in animal evolution. However, several organisms that are more closely related to flies than to vertebrates have eye photoreceptor cells with cilia. Do all eye photoreceptors with cilia have a common origin in evolution or did they emerge independently in vertebrates and certain invertebrates?

The photoreceptor cells of a marine mollusc called Leptochiton asellus, are unusual because they bear both microvilli and cilia, suggesting they have intermediate characteristics between the two well-known types of photoreceptor cells. Previous studies have shown that these photoreceptor cells use r-opsin, but Vöcking et al. have now detected the presence of an additional opsin in the cells. This opsin is a member of the recently discovered xenopsin family of molecules. Further analyses support the findings of previous studies that suggested this type of opsin emerged early on in animal evolution, independently from c-opsin. Other invertebrates that have cilia on their eye photoreceptors also use xenopsin and not c-opsin.

The findings of Vöcking et al. suggest that, in addition to c-opsin and r-opsin, xenopsin has also driven the evolution of photoreceptor cells in animals. Eye photoreceptor cells in invertebrates with cilia probably share a common origin with the microvilli photoreceptor cells that is distinct from that of vertebrate visual cells. The observation that two very different types of opsin can be produced within a single cell suggests that the molecular processes that respond to light in photoreceptor cells may be much more complex than previously anticipated. Further work on these processes may help us to understand how animal eyes work and how they are affected by disease.

Genome-wide DNA Methylation Analysis Reveals GABBR2 as a Novel Epigenetic Target for EGFR 19 Deletion Lung Adenocarcinoma with Induction Erlotinib Treatment

Purpose: The past decade has witnessed the rapid development of personalized targeted therapies in lung cancer. It is still unclear whether epigenetic changes are involved in the response to tyrosine kinase inhibitor (TKI) treatment in epidermal growth factor receptor (EGFR)-mutated lung cancer.Experimental Design: Methyl-sensitive cut counting sequencing (MSCC) was applied to investigate the methylation changes in paired tissues before and after erlotinib treatment for 42 days with partial response (PR) from stage IIIa (N2) lung adenocarcinoma patients (N = 2) with EGFR 19 deletion. The Sequenom EpiTYPER assay was used to validate the changed methylated candidate genes. Up- or downregulation of the candidate gene was performed to elucidate the potential mechanism in the regulation of erlotinib treatment response.Results: Sixty aberrant methylated genes were screened using MSCC sequencing. Two aberrant methylated genes, CBFA2T3 and GABBR2, were clearly validated. A same differential methylated region (DMR) between exon 2 and exon 3 of GABBR2 gene was confirmed consistently in both patients. GABBR2 was significantly downregulated in EGFR 19 deletion cells, HCC4006 and HCC827, but remained conserved in EGFR wild-type A549 cells after erlotinib treatment. Upregulation of GABBR2 expression significantly rescued erlotinib-induced apoptosis in HCC827 cells. GABBR2 was significantly downregulated, along with the reduction of S6, p-p70 S6, and p-ERK1/2, demonstrating that GABBR2 may play an important role in EGFR signaling through the ERK1/2 pathway.Conclusions: We demonstrated that GABBR2 gene might be a novel potential epigenetic treatment target with induction erlotinib treatment for stage IIIa (N2) EGFR 19 deletion lung adenocarcinoma. Clin Cancer Res; 23(17); 5003-14. ©2017 AACR.

Alternative exon-encoding regions of Locusta migratoria muscle myosin modulate the pH dependence of ATPase activity

Whereas the vertebrate muscle myosin heavy chains (MHCs) are encoded by a family of Mhc genes, most insects examined to date contain a single Mhc gene and produce all of the different MHC isoforms by alternative RNA splicing. Here, we found that the migratory locust, Locusta migratoria, has one Mhc gene, which contains 41 exons, including five alternative exclusive exons and one differently included penultimate exon, and potentially encodes 360 MHC isoforms. From the adult L. migratoria, we identified 14 MHC isoforms (including two identical isoforms): four from flight muscle (the thorax dorsal longitudinal muscle), three from jump muscle (the hind leg extensor tibiae muscle) and seven from the abdominal intersegmental muscle. We purified myosins from flight muscle and jump muscle and characterized their motor activities. At neutral pH, the flight and the jump muscle myosins displayed similar levels of in vitro actin-gliding activity, whereas the former had a slightly higher actin-activated ATPase activity than the latter. Interestingly, the pH dependences of the actin-activated ATPase activity of these two myosins are different. Because the dominant MHC isoforms in these two muscles are identical except for the two alternative exon-encoding regions, we propose that these two alternative regions modulate the pH dependence of L. migratoria muscle myosin.

Alternative haplotypes of antigen processing genes in zebrafish diverged early in vertebrate evolution

Antigen processing and presentation genes found within the MHC are among the most highly polymorphic genes of vertebrate genomes, providing populations with diverse immune responses to a wide array of pathogens. Here, we describe transcriptome, exome, and whole-genome sequencing of clonal zebrafish, uncovering the most extensive diversity within the antigen processing and presentation genes of any species yet examined. Our CG2 clonal zebrafish assembly provides genomic context within a remarkably divergent haplotype of the core MHC region on chromosome 19 for six expressed genes not found in the zebrafish reference genome: mhc1uga, proteasome-β 9b (psmb9b), psmb8f, and previously unknown genes psmb13b, tap2d, and tap2e We identify ancient lineages for Psmb13 within a proteasome branch previously thought to be monomorphic and provide evidence of substantial lineage diversity within each of three major trifurcations of catalytic-type proteasome subunits in vertebrates: Psmb5/Psmb8/Psmb11, Psmb6/Psmb9/Psmb12, and Psmb7/Psmb10/Psmb13. Strikingly, nearby tap2 and MHC class I genes also retain ancient sequence lineages, indicating that alternative lineages may have been preserved throughout the entire MHC pathway since early diversification of the adaptive immune system ∼500 Mya. Furthermore, polymorphisms within the three MHC pathway steps (antigen cleavage, transport, and presentation) are each predicted to alter peptide specificity. Lastly, comparative analysis shows that antigen processing gene diversity is far more extensive than previously realized (with ancient coelacanth psmb8 lineages, shark psmb13, and tap2t and psmb10 outside the teleost MHC), implying distinct immune functions and conserved roles in shaping MHC pathway evolution throughout vertebrates.

The paralog-to-contig assignment problem: high quality gene models from fragmented assemblies

Background

The accurate annotation of genes in newly sequenced genomes remains a challenge. Although sophisticated comparative pipelines are available, computationally derived gene models are often less than perfect. This is particularly true when multiple similar paralogs are present. The issue is aggravated further when genomes are assembled only at a preliminary draft level to contigs or short scaffolds. However, these genomes deliver valuable information for studying gene families. High accuracy models of protein coding genes are needed in particular for phylogenetics and for the analysis of gene family histories.

Results

We present a pipeline, ExonMatchSolver, that is designed to help the user to produce and curate high quality models of the protein-coding part of genes. The tool in particular tackles the problem of identifying those coding exon groups that belong to the same paralogous genes in a fragmented genome assembly. This paralog-to-contig assignment problem is shown to be NP-complete. It is phrased and solved as an Integer Linear Programming problem.

Conclusions

The ExonMatchSolver-pipeline can be employed to build highly accurate models of protein coding genes even when spanning several genomic fragments. This sets the stage for a better understanding of the evolutionary history within particular gene families which possess a large number of paralogs and in which frequent gene duplication events occurred.

Electronic supplementary material

The online version of this article (doi:10.1186/s13015-016-0063-y) contains supplementary material, which is available to authorized users.

Tricholoma vaccinum host communication during ectomycorrhiza formation

The genome sequence of Tricholoma vaccinum was obtained to predict its secretome in order to elucidate communication of T. vaccinum with its host tree spruce (Picea abies) in interkingdom signaling. The most prominent protein domains within the 206 predicted secreted proteins belong to energy and nutrition (52%), cell wall degradation (19%) and mycorrhiza establishment (9%). Additionally, we found small secreted proteins that show typical features of effectors potentially involved in host communication. From the secretome, 22 proteins could be identified, two of which showed higher protein abundances after spruce root exudate exposure, while five were downregulated in this treatment. The changes in T. vaccinum protein excretion with first recognition of the partner were used to identify small secreted proteins with the potential to act as effectors in the mutually beneficial symbiosis. Our observations support the hypothesis of a complex communication network including a cocktail of communication molecules induced long before physical contact of the partners.

Polyphyly of nuclear lamin genes indicates an early eukaryotic origin of the metazoan-type intermediate filament proteins

The nuclear lamina is a protein meshwork associated with the inner side of the nuclear envelope contributing structural, signalling and regulatory functions. Here, I report on the evolution of an important component of the lamina, the lamin intermediate filament proteins, across the eukaryotic tree of life. The lamins show a variety of protein domain and sequence motif architectures beyond the classical α-helical rod, nuclear localisation signal, immunoglobulin domain and CaaX motif organisation, suggesting extension and adaptation of functions in many species. I identified lamin genes not only in metazoa and Amoebozoa as previously described, but also in other opisthokonts including Ichthyosporea and choanoflagellates, in oomycetes, a sub-family of Stramenopiles, and in Rhizaria, implying that they must have been present very early in eukaryotic evolution if not even the last common ancestor of all extant eukaryotes. These data considerably extend the current perception of lamin evolution and have important implications with regard to the evolution of the nuclear envelope.

Horizontal transfer of carbohydrate metabolism genes into ectomycorrhizal Amanita

The genus Amanita encompasses both symbiotic, ectomycorrhizal fungi and asymbiotic litter decomposers; all species are derived from asymbiotic ancestors. Symbiotic species are no longer able to degrade plant cell walls. The carbohydrate esterases family 1 (CE1s) is a diverse group of enzymes involved in carbon metabolism, including decomposition and carbon storage. CE1 genes of the ectomycorrhizal A. muscaria appear diverged from all other fungal homologues, and more similar to CE1s of bacteria, suggesting a horizontal gene transfer (HGT) event. In order to test whether AmanitaCE1s were acquired horizontally, we built a phylogeny of CE1s collected from across the tree of life, and describe the evolution of CE1 genes among Amanita and relevant lineages of bacteria. CE1s of symbiotic Amanita were very different from CE1s of asymbiotic Amanita, and are more similar to bacterial CE1s. The protein structure of one CE1 gene of A. muscaria matched a depolymerase that degrades the carbon storage molecule poly((R)-3-hydroxybutyrate) (PHB). Asymbiotic Amanita do not carry sequence or structural homologues of these genes. The CE1s acquired through HGT may enable novel metabolisms, or play roles in signaling or defense. This is the first evidence for the horizontal transfer of carbohydrate metabolism genes into ectomycorrhizal fungi.

Six subgroups and extensive recent duplications characterize the evolution of the eukaryotic tubulin protein family

Tubulins belong to the most abundant proteins in eukaryotes providing the backbone for many cellular substructures like the mitotic and meiotic spindles, the intracellular cytoskeletal network, and the axonemes of cilia and flagella. Homologs have even been reported for archaea and bacteria. However, a taxonomically broad and whole-genome-based analysis of the tubulin protein family has never been performed, and thus, the number of subfamilies, their taxonomic distribution, and the exact grouping of the supposed archaeal and bacterial homologs are unknown. Here, we present the analysis of 3,524 tubulins from 504 species. The tubulins formed six major subfamilies, α to ζ. Species of all major kingdoms of the eukaryotes encode members of these subfamilies implying that they must have already been present in the last common eukaryotic ancestor. The proposed archaeal homologs grouped together with the bacterial TubZ proteins as sister clade to the FtsZ proteins indicating that tubulins are unique to eukaryotes. Most species contained α- and/or β-tubulin gene duplicates resulting from recent branch- and species-specific duplication events. This shows that tubulins cannot be used for constructing species phylogenies without resolving their ortholog–paralog relationships. The many gene duplicates and also the independent loss of the δ-, ε-, or ζ-tubulins, which have been shown to be part of the triplet microtubules in basal bodies, suggest that tubulins can functionally substitute each other.

Molecular phylogeny of sequenced Saccharomycetes reveals polyphyly of the alternative yeast codon usage

The universal genetic code defines the translation of nucleotide triplets, called codons, into amino acids. In many Saccharomycetes a unique alteration of this code affects the translation of the CUG codon, which is normally translated as leucine. Most of the species encoding CUG alternatively as serine belong to the Candida genus and were grouped into a so-called CTG clade. However, the “Candida genus” is not a monophyletic group and several Candida species are known to use the standard CUG translation. The codon identity could have been changed in a single branch, the ancestor of the Candida, or to several branches independently leading to a polyphyletic alternative yeast codon usage (AYCU). In order to resolve the monophyly or polyphyly of the AYCU, we performed a phylogenomics analysis of 26 motor and cytoskeletal proteins from 60 sequenced yeast species. By investigating the CUG codon positions with respect to sequence conservation at the respective alignment positions, we were able to unambiguously assign the standard code or AYCU. Quantitative analysis of the highly conserved leucine and serine alignment positions showed that 61.1% and 17% of the CUG codons coding for leucine and serine, respectively, are at highly conserved positions, whereas only 0.6% and 2.3% of the CUG codons, respectively, are at positions conserved in the respective other amino acid. Plotting the codon usage onto the phylogenetic tree revealed the polyphyly of the AYCU with Pachysolen tannophilus and the CTG clade branching independently within a time span of 30–100 Ma.

Kassiopeia: a database and web application for the analysis of mutually exclusive exomes of eukaryotes

BACKGROUND

Alternative splicing is an important process in higher eukaryotes that allows obtaining several transcripts from one gene. A specific case of alternative splicing is mutually exclusive splicing, in which exactly one exon out of a cluster of neighbouring exons is spliced into the mature transcript. Recently, a new algorithm for the prediction of these exons has been developed based on the preconditions that the exons of the cluster have similar lengths, sequence homology, and conserved splice sites, and that they are translated in the same reading frame.

DESCRIPTION

In this contribution we introduce Kassiopeia, a database and web application for the generation, storage, and presentation of genome-wide analyses of mutually exclusive exomes. Currently, Kassiopeia provides access to the mutually exclusive exomes of twelve Drosophila species, the thale cress Arabidopsis thaliana, the flatworm Caenorhabditis elegans, and human. Mutually exclusive spliced exons (MXEs) were predicted based on gene reconstructions from Scipio. Based on the standard prediction values, with which 83.5% of the annotated MXEs of Drosophila melanogaster were reconstructed, the exomes contain surprisingly more MXEs than previously supposed and identified. The user can search Kassiopeia using BLAST or browse the genes of each species optionally adjusting the parameters used for the prediction to reveal more divergent or only very similar exon candidates.

CONCLUSIONS

We developed a pipeline to predict MXEs in the genomes of several model organisms and a web interface, Kassiopeia, for their visualization. For each gene Kassiopeia provides a comprehensive gene structure scheme, the sequences and predicted secondary structures of the MXEs, and, if available, further evidence for MXE candidates from cDNA/EST data, predictions of MXEs in homologous genes of closely related species, and RNA secondary structure predictions. Kassiopeia can be accessed at http://www.motorprotein.de/kassiopeia.

Shared gene structures and clusters of mutually exclusive spliced exons within the metazoan muscle myosin heavy chain genes

Multicellular animals possess two to three different types of muscle tissues. Striated muscles have considerable ultrastructural similarity and contain a core set of proteins including the muscle myosin heavy chain (Mhc) protein. The ATPase activity of this myosin motor protein largely dictates muscle performance at the molecular level. Two different solutions to adjusting myosin properties to different muscle subtypes have been identified so far: Vertebrates and nematodes contain many independent differentially expressed Mhc genes while arthropods have single Mhc genes with clusters of mutually exclusive spliced exons (MXEs). The availability of hundreds of metazoan genomes now allowed us to study whether the ancient bilateria already contained MXEs, how MXE complexity subsequently evolved, and whether additional scenarios to control contractile properties in different muscles could be proposed, By reconstructing the Mhc genes from 116 metazoans we showed that all intron positions within the motor domain coding regions are conserved in all bilateria analysed. The last common ancestor of the bilateria already contained a cluster of MXEs coding for part of the loop-2 actin-binding sequence. Subsequently the protostomes and later the arthropods gained many further clusters while MXEs got completely lost independently in several branches (vertebrates and nematodes) and species (for example the annelid Helobdella robusta and the salmon louse Lepeophtheirus salmonis). Several bilateria have been found to encode multiple Mhc genes that might all or in part contain clusters of MXEs. Notable examples are a cluster of six tandemly arrayed Mhc genes, of which two contain MXEs, in the owl limpet Lottia gigantea and four Mhc genes with three encoding MXEs in the predatory mite Metaseiulus occidentalis. Our analysis showed that similar solutions to provide different myosin isoforms (multiple genes or clusters of MXEs or both) have independently been developed several times within bilaterian evolution.

Whole genome duplication events in plant evolution reconstructed and predicted using myosin motor proteins

Background

The evolution of land plants is characterized by whole genome duplications (WGD), which drove species diversification and evolutionary novelties. Detecting these events is especially difficult if they date back to the origin of the plant kingdom. Established methods for reconstructing WGDs include intra- and inter-genome comparisons, K_S age distribution analyses, and phylogenetic tree constructions.

Results

By analysing 67 completely sequenced plant genomes 775 myosins were identified and manually assembled. Phylogenetic trees of the myosin motor domains revealed orthologous and paralogous relationships and were consistent with recent species trees. Based on the myosin inventories and the phylogenetic trees, we have identified duplications of the entire myosin motor protein family at timings consistent with 23 WGDs, that had been reported before. We also predict 6 WGDs based on further protein family duplications. Notably, the myosin data support the two recently reported WGDs in the common ancestor of all extant angiosperms. We predict single WGDs in the Manihot esculenta and Nicotiana benthamiana lineages, two WGDs for Linum usitatissimum and Phoenix dactylifera, and a triplication or two WGDs for Gossypium raimondii. Our data show another myosin duplication in the ancestor of the angiosperms that could be either the result of a single gene duplication or a remnant of a WGD.

Conclusions

We have shown that the myosin inventories in angiosperms retain evidence of numerous WGDs that happened throughout plant evolution. In contrast to other protein families, many myosins are still present in extant species. They are closely related and have similar domain architectures, and their phylogenetic grouping follows the genome duplications. Because of its broad taxonomic sampling the dataset provides the basis for reliable future identification of further whole genome duplications.

WebScipio: Reconstructing alternative splice variants of eukaryotic proteins

Accurate exon–intron structures are essential prerequisites in genomics, proteomics and for many protein family and single gene studies. We originally developed Scipio and the corresponding web service WebScipio for the reconstruction of gene structures based on protein sequences and available genome assemblies. WebScipio also allows predicting mutually exclusive spliced exons and tandemly arrayed gene duplicates. The obtained gene structures are illustrated in graphical schemes and can be analysed down to the nucleotide level. The set of eukaryotic genomes available at the WebScipio server is updated on a daily basis. The current version of the web server provides access to ∼3400 genome assembly files of >1100 sequenced eukaryotic species. Here, we have also extended the functionality by adding a module with which expressed sequence tag (EST) and cDNA data can be mapped to the reconstructed gene structure for the identification of all types of alternative splice variants. WebScipio has a user-friendly web interface, and we believe that the improved web server will provide better service to biologists interested in the gene structure corresponding to their protein of interest, including all types of alternative splice forms and tandem gene duplicates. WebScipio is freely available at http://www.webscipio.org.

GenePainter: a fast tool for aligning gene structures of eukaryotic protein families, visualizing the alignments and mapping gene structures onto protein structures

Background

All sequenced eukaryotic genomes have been shown to possess at least a few introns. This includes those unicellular organisms, which were previously suspected to be intron-less. Therefore, gene splicing must have been present at least in the last common ancestor of the eukaryotes. To explain the evolution of introns, basically two mutually exclusive concepts have been developed. The introns-early hypothesis says that already the very first protein-coding genes contained introns while the introns-late concept asserts that eukaryotic genes gained introns only after the emergence of the eukaryotic lineage. A very important aspect in this respect is the conservation of intron positions within homologous genes of different taxa.

Results

GenePainter is a standalone application for mapping gene structure information onto protein multiple sequence alignments. Based on the multiple sequence alignments the gene structures are aligned down to single nucleotides. GenePainter accounts for variable lengths in exons and introns, respects split codons at intron junctions and is able to handle sequencing and assembly errors, which are possible reasons for frame-shifts in exons and gaps in genome assemblies. Thus, even gene structures of considerably divergent proteins can properly be compared, as it is needed in phylogenetic analyses. Conserved intron positions can also be mapped to user-provided protein structures. For their visualization GenePainter provides scripts for the molecular graphics system PyMol.

Conclusions

GenePainter is a tool to analyse gene structure conservation providing various visualization options. A stable version of GenePainter for all operating systems as well as documentation and example data are available at http://www.motorprotein.de/genepainter.html.

Calcitonin-like diuretic hormones in insects

Insect neuropeptides control various biological processes including growth, development, homeostasis and reproduction. The calcitonin-like diuretic hormone (CT/DH) is one such neuropeptide that has been shown to affect salt and water transport by Malpighian tubules of several insects. With an increase in the number of sequenced insect genomes, CT/DHs have been predicted in several insect species, making it easier to characterize the gene encoding this hormone and determine its function in the species in question. This mini review summarizes the current knowledge on insect CT/DHs, focusing on mRNA and peptide structures, distribution patterns, physiological roles, and receptors in insects.

Evolution of the eukaryotic dynactin complex, the activator of cytoplasmic dynein

Background

Dynactin is a large multisubunit protein complex that enhances the processivity of cytoplasmic dynein and acts as an adapter between dynein and the cargo. It is composed of eleven different polypeptides of which eight are unique to this complex, namely dynactin1 (p150^Glued), dynactin2 (p50 or dynamitin), dynactin3 (p24), dynactin4 (p62), dynactin5 (p25), dynactin6 (p27), and the actin-related proteins Arp1 and Arp10 (Arp11).

Results

To reveal the evolution of dynactin across the eukaryotic tree the presence or absence of all dynactin subunits was determined in most of the available eukaryotic genome assemblies. Altogether, 3061 dynactin sequences from 478 organisms have been annotated. Phylogenetic trees of the various subunit sequences were used to reveal sub-family relationships and to reconstruct gene duplication events. Especially in the metazoan lineage, several of the dynactin subunits were duplicated independently in different branches. The largest subunit repertoire is found in vertebrates. Dynactin diversity in vertebrates is further increased by alternative splicing of several subunits. The most prominent example is the dynactin1 gene, which may code for up to 36 different isoforms due to three different transcription start sites and four exons that are spliced as differentially included exons.

Conclusions

The dynactin complex is a very ancient complex that most likely included all subunits in the last common ancestor of extant eukaryotes. The absence of dynactin in certain species coincides with that of the cytoplasmic dynein heavy chain: Organisms that do not encode cytoplasmic dynein like plants and diplomonads also do not encode the unique dynactin subunits. The conserved core of dynactin consists of dynactin1, dynactin2, dynactin4, dynactin5, Arp1, and the heterodimeric actin capping protein. The evolution of the remaining subunits dynactin3, dynactin6, and Arp10 is characterized by many branch- and species-specific gene loss events.