Search

Background: Tubulins play crucial roles in numerous fundamental processes of plant development. In flowering plants, tubulins are grouped into α-, β- and γ-subfamilies, while α- and β-tubulins possess a large isotype diversity and gene number variations among different species. This circumstance leads to insufficient recognition of orthologous isotypes and significantly complicates extrapolation of obtained experimental results, and brings difficulties for the identification of particular tubulin isotype function. The aim of this research is to identify and characterize tubulins of an emerging biofuel crop Camelina sativa. Results: We report comprehensive identification and characterization of tubulin gene family in C. sativa, including analyses of exon-intron organization, duplicated genes comparison, proper isotype designation, phylogenetic analysis, and expression patterns in different tissues. 17 α-, 34 β- and 6 γ-tubulin genes were identified and assigned to a particular isotype. Recognition of orthologous tubulin isotypes was cross-referred, involving data of phylogeny, synteny analyses and genes allocation on reconstructed genomic blocks of Ancestral Crucifer Karyotype. An investigation of expression patterns of tubulin homeologs revealed the predominant role of N6 (A) and N7 (B) subgenomes in tubulin expression at various developmental stages, contrarily to general the dominance of transcripts of H7 (C) subgenome. Conclusions: For the first time a complete set of tubulin gene family members was identified and characterized for allohexaploid C. sativa species. The study demonstrates the comprehensive approach of precise inferring gene orthology. The applied technique allowed not only identifying C. sativa tubulin orthologs in model Arabidopsis species and tracking tubulin gene evolution, but also uncovered that A. thaliana is missing orthologs for several particular isotypes of α- and β-tubulins. [ABSTRACT FROM AUTHOR]

Background: The present availability of full genome sequences of a broad range of animal species across the whole range of evolutionary history enables one to ask questions as to the distribution of genes across the chromosomes. Do newly recruited genes, as new clades emerge, distribute at random or at non-random locations? Results: We extracted values for the ages of the human genes and for their current chromosome locations, from published sources. A quantitative analysis showed that the distribution of newly-added genes among and within the chromosomes appears to be increasingly non-random if one observes animals along the evolutionary series from the precursors of the tetrapoda through to the great apes, whereas the oldest genes are randomly distributed. Conclusions: Randomization will result from chromosome evolution, but less and less time is available for this process as evolution proceeds. Much of the bunching of recently-added genes arises from new gene formation as paralogues in gene families, near the location of genes that were recruited in the preceding phylostratum. As examples we cite the KRTAP, ZNF, OR and some minor gene families. We show that bunching can also result from the evolution of the chromosomes themselves when, as for the KRTAP genes, blocks of genes that had previously been on disparate chromosomes become linked together. [ABSTRACT FROM AUTHOR]

Background: Small nucleolar RNAs (snoRNAs) are an abundant class of noncoding RNAs present in all eukaryotes and best known for their involvement in ribosome biogenesis. In mammalian genomes, many snoRNAs exist in multiple copies, resulting from recombination and retrotransposition from an ancestral snoRNA. To gain insight into snoRNA copy regulation, we used Rfam classification and normal human tissue expression datasets generated using low structure bias RNA-seq to characterize snoRNA families. Results: We found that although box H/ACA families are on average larger than box C/D families, the number of expressed members is similar for both types. Family members can cover a wide range of average abundance values, but importantly, expression variability of individual members of a family is preferred over the total variability of the family, especially for box H/ACA snoRNAs, suggesting that while members are likely differentially regulated, mechanisms exist to ensure uniformity of the total family abundance across tissues. Box C/D snoRNA family members are mostly embedded in the same host gene while box H/ACA family members tend to be encoded in more than one different host, supporting a model in which box C/D snoRNA duplication occurred mostly by cis recombination while box H/ACA snoRNA families have gained copy members through retrotransposition. And unexpectedly, snoRNAs encoded in the same host gene can be regulated independently, as some snoRNAs within the same family vary in abundance in a divergent way between tissues. Conclusions: SnoRNA copy regulation affects family sizes, genomic location of the members and controls simultaneously member and total family abundance to respond to the needs of individual tissues. [ABSTRACT FROM AUTHOR]

Background: Plant pathogenesis related-1 (PR-1) proteins belong to the CAP superfamily and have been characterized as markers of induced defense against pathogens. Moniliophthora perniciosa and Moniliophthora roreri are hemibiotrophic fungi that respectively cause the witches' broom disease and frosty pod rot in Theobroma cacao. Interestingly, a large number of plant PR-1-like genes are present in the genomes of both species and many are up-regulated during the biotrophic interaction. In this study, we investigated the evolution of PR-1 proteins from 22 genomes of Moniliophthora isolates and 16 other Agaricales species, performing genomic investigation, phylogenetic reconstruction, positive selection search and gene expression analysis. Results: Phylogenetic analysis revealed conserved PR-1 genes (PR-1a, b, d, j), shared by many Agaricales saprotrophic species, that have diversified in new PR-1 genes putatively related to pathogenicity in Moniliophthora (PR-1f, g, h, i), as well as in recent specialization cases within M. perniciosa biotypes (PR-1c, k, l) and M. roreri (PR-1n). PR-1 families in Moniliophthora with higher evolutionary rates exhibit induced expression in the biotrophic interaction and positive selection clues, supporting the hypothesis that these proteins accumulated adaptive changes in response to host-pathogen arms race. Furthermore, although previous work showed that MpPR-1 can detoxify plant antifungal compounds in yeast, we found that in the presence of eugenol M. perniciosa differentially expresses only MpPR-1e, k, d, of which two are not linked to pathogenicity, suggesting that detoxification might not be the main function of most MpPR-1. Conclusions: Based on analyses of genomic and expression data, we provided evidence that the evolution of PR-1 in Moniliophthora was adaptive and potentially related to the emergence of the parasitic lifestyle in this genus. Additionally, we also discuss how fungal PR-1 proteins could have adapted from basal conserved functions to possible roles in fungal pathogenesis. [ABSTRACT FROM AUTHOR]

Background: The SARS-CoV-2 pandemic is one of the greatest global medical and social challenges that have emerged in recent history. Human coronavirus strains discovered during previous SARS outbreaks have been hypothesized to pass from bats to humans using intermediate hosts, e.g. civets for SARS-CoV and camels for MERS-CoV. The discovery of an intermediate host of SARS-CoV-2 and the identification of specific mechanism of its emergence in humans are topics of primary evolutionary importance. In this study we investigate the evolutionary patterns of 11 main genes of SARS-CoV-2. Previous studies suggested that the genome of SARS-CoV-2 is highly similar to the horseshoe bat coronavirus RaTG13 for most of the genes and to some Malayan pangolin coronavirus (CoV) strains for the receptor binding (RB) domain of the spike protein. Results: We provide a detailed list of statistically significant horizontal gene transfer and recombination events (both intergenic and intragenic) inferred for each of 11 main genes of the SARS-CoV-2 genome. Our analysis reveals that two continuous regions of genes S and N of SARS-CoV-2 may result from intragenic recombination between RaTG13 and Guangdong (GD) Pangolin CoVs. Statistically significant gene transfer-recombination events between RaTG13 and GD Pangolin CoV have been identified in region [1215-1425] of gene S and region [534-727] of gene N. Moreover, some statistically significant recombination events between the ancestors of SARS-CoV-2, RaTG13, GD Pangolin CoV and bat CoV ZC45-ZXC21 coronaviruses have been identified in genes ORF1ab, S, ORF3a, ORF7a, ORF8 and N. Furthermore, topology-based clustering of gene trees inferred for 25 CoV organisms revealed a three-way evolution of coronavirus genes, with gene phylogenies of ORF1ab, S and N forming the first cluster, gene phylogenies of ORF3a, E, M, ORF6, ORF7a, ORF7b and ORF8 forming the second cluster, and phylogeny of gene ORF10 forming the third cluster. Conclusions: The results of our horizontal gene transfer and recombination analysis suggest that SARS-CoV-2 could not only be a chimera virus resulting from recombination of the bat RaTG13 and Guangdong pangolin coronaviruses but also a close relative of the bat CoV ZC45 and ZXC21 strains. They also indicate that a GD pangolin may be an intermediate host of this dangerous virus. [ABSTRACT FROM AUTHOR]

Background: hes genes are chordate homologs of Drosophila genes, hairy and enhancer of split, which encode a basic helix-loop-helix (bHLH) transcriptional repressor with a WRPW motif. Various developmental functions of hes genes, including early embryogenesis and neurogenesis, have been elucidated in vertebrates. However, their orthologous relationships remain unclear partly because of less conservation of relatively short amino acid sequences, the fact that the genome was not analyzed as it is today, and species-specific genome duplication. This results in complicated gene names in vertebrates, which are not consistent in orthologs. We previously revealed that Xenopus frogs have two clusters of hes5, named “the hes5.1 cluster” and “the hes5.3 cluster”, but the origin and the conservation have not yet been revealed. Results: Here, we elucidated the orthologous and paralogous relationships of all hes genes of human, mouse, chicken, gecko, zebrafish, medaka, coelacanth, spotted gar, elephant shark and three species of frogs, Xenopus tropicalis (X. tropicalis), X. laevis, Nanorana parkeri, by phylogenetic and synteny analyses. Any duplicated hes5 were not found in mammals, whereas hes5 clusters in teleost were conserved although not as many genes as the three frog species. In addition, hes5 cluster-like structure was found in the elephant shark genome, but not found in cyclostomata. Conclusion: These data suggest that the hes5 cluster existed in the gnathostome ancestor but became a single gene in mammals. The number of hes5 cluster genes were specifically large in frogs. [ABSTRACT FROM AUTHOR]

Background: Recent advances in genetics and genomics present unique opportunities for enhancing our understanding of mammalian biology and evolution through detailed multi-species comparative analysis of gene organization and expression. Yet, of the more than 20,000 protein coding genes found in mammalian genomes, fewer than 10% have been examined in any detail. Here we elucidate the power of data available in publicly-accessible genomic and genetic resources by querying them to evaluate Zmat2, a minimally studied gene whose human ortholog has been implicated in spliceosome function and in keratinocyte differentiation. Results: We find extensive conservation in coding regions and overall structure of Zmat2 in 18 mammals representing 13 orders and spanning ~ 165 million years of evolutionary development, and in their encoded proteins. We identify a tandem duplication in the Zmat2 gene and locus in opossum, but not in other monotremes, marsupials, or other mammals, indicating that this event occurred subsequent to the divergence of these species from one another. We also define a collection of Zmat2 pseudogenes in half of the mammals studied, and suggest based on phylogenetic analysis that they each arose independently in the recent evolutionary past. Conclusions: Mammalian Zmat2 genes and ZMAT2 proteins illustrate conservation of structure and sequence, along with the development and diversification of pseudogenes in a large fraction of species. Collectively, these observations also illustrate how the focused identification and interpretation of data found in public genomic and gene expression resources can be leveraged to reveal new insights of potentially high biological significance. [ABSTRACT FROM AUTHOR]

Background: Wheat grain avenin-like proteins (ALPs) belong to a recently discovered class of wheat grain storage protein. ALPs in wheat grains not only have beneficial effects on dough quality but also display antifungal activities, which is a novel observation for wheat storage proteins. Previous studies have shown that ALPs are likely present in the albumin/globulin fractions of total protein extract from wheat flour. However, the accumulation characteristics of these ALPs in the mature wheat grain remains unknown. Results: In the present study, a total of 13 ALPs homologs were isolated and characterized in the albumin/globulin fractions of the wheat protein extract. A combination of multiple techniques including RP-HPLC, SDS-PAGE, MALDI-TOF and peptide sequencing were used for accurate separation and identification of individual ALP homolog. The C-terminal TaALP-by-4AL/7DS, TaALP-by-4AL/7AS/7DS, TaALP-bx/4AL/7AS/7DS, TaALP-ay-7DS, TaALP-ay-4AL, TaALP-ax-4AL, TaALP-ax-7AS, and TaALP-ax-7DS, were separated as individual protein bands from wheat flour for the first time. These unique ALPs peptides were mapped to the latest wheat genome assembly in the IWGSC database. The characteristic defence related proteins present in albumin and globulin fractions, such as protein disulfide-isomerase (PDI), grain softness protein (GSP), alpha-amylase inhibitors (AAIs) and endogenous alpha-amylase/subtilisin inhibitor were also found to co-segregate with these identified ALPs, avenin-3 and α-gliadins. The molecular weight range and the electrophoresis segregation properties of ALPs were characterised in comparison with the proteins containing the tryp_alpha_amyl domain (PF00234) and the gliadin domain (PF13016), which play a role in plant immunity and grain quality. We examined the phylogenetic relationships of the AAIs, GSP, avenin-3, α-gliadins and ALPs, based on the alignment of their functional domains. MALDI-TOF profiling indicated the occurrence of certain post-translations modifications (PTMs) in some ALP subunits. Conclusions: We reported for the first time the complete profiling of ALPs present in the albumin/globulin fractions of wheat grain protein extracts. We concluded that majority of the ALPs homologs are expressed in wheat grains. We found clear evidence of PTMs in several ALPs peptides. The identification of both gliadin domain (PF13016) and Tryp_alpha_amyl domain (PF00234) in the mature forms of ALPs highlighted the multiple functional properties of ALPs in grain quality and disease resistance. [ABSTRACT FROM AUTHOR]

Background: Sugarcane served as the model plant for discovery of the C4 photosynthetic pathway. Magnesium is the central atom of chlorophyll, and thus is considered as a critical nutrient for plant development and photosynthesis. In plants, the magnesium transporter (MGT) family is composed of a number of membrane proteins, which play crucial roles in maintaining Mg homeostasis. However, to date there is no information available on the genomics of MGTs in sugarcane due to the complexity of the Saccharum genome. Results: Here, we identified 10 MGTs from the Saccharum spontaneum genome. Phylogenetic analysis of MGTs suggested that the MGTs contained at least 5 last common ancestors before the origin of angiosperms. Gene structure analysis suggested that MGTs family of dicotyledon may be accompanied by intron loss and pseudoexon phenomena during evolution. The pairwise synonymous substitution rates corresponding to a divergence time ranged from 142.3 to 236.6 Mya, demonstrating that the MGTs are an ancient gene family in plants. Both the phylogeny and Ks analyses indicated that SsMGT1/SsMGT2 originated from the recent ρWGD, and SsMGT7/SsMGT8 originated from the recent σ WGD. These 4 recently duplicated genes were shown low expression levels and assumed to be functionally redundant. MGT6, MGT9 and MGT10 weredominant genes in the MGT family and werepredicted to be located inthe chloroplast. Of the 3 dominant MGTs, SsMGT6 expression levels were found to be induced in the light period, while SsMGT9 and SsMTG10 displayed high expression levels in the dark period. These results suggested that SsMGT6 may have a function complementary to SsMGT9 and SsMTG10 that follows thecircadian clock for MGT in the leaf tissues of S. spontaneum. MGT3, MGT7 and MGT10 had higher expression levels Insaccharum officinarum than in S. spontaneum, suggesting their functional divergence after the split of S. spontaneum and S. officinarum. Conclusions: This study of gene evolution and expression of MGTs in S. spontaneum provided basis for the comprehensive genomic study of the entire MGT genes family in Saccharum. The results are valuable for further functional analyses of MGT genes and utilization of the MGTs for Saccharum genetic improvement. [ABSTRACT FROM AUTHOR]

Background: Spanins are phage lysis proteins required to disrupt the outer membrane. Phages employ either two-component spanins or unimolecular spanins in this final step of Gram-negative host lysis. Two-component spanins like Rz-Rz1 from phage lambda consist of an integral inner membrane protein: i-spanin, and an outer membrane lipoprotein: o-spanin, that form a complex spanning the periplasm. Two-component spanins exist in three different genetic architectures; embedded, overlapped and separated. In contrast, the unimolecular spanins, like gp11 from phage T1, have an N-terminal lipoylation signal sequence and a C-terminal transmembrane domain to account for the topology requirements. Our proposed model for spanin function, for both spanin types, follows a common theme of the outer membrane getting fused with the inner membrane, effecting the release of progeny virions. Results: Here we present a SpaninDataBase which consists of 528 two-component spanins and 58 unimolecular spanins identified in this analysis. Primary analysis revealed significant differences in the secondary structure predictions for the periplasmic domains of the two-component and unimolecular spanin types, as well as within the three different genetic architectures of the two-component spanins. Using a threshold of 40% sequence identity over 40% sequence length, we were able to group the spanins into 143 i-spanin, 125 o-spanin and 13 u-spanin families. More than 40% of these families from each type were singletons, underlining the extreme diversity of this class of lysis proteins. Multiple sequence alignments of periplasmic domains demonstrated conserved secondary structure patterns and domain organization within family members. Furthermore, analysis of families with members from different architecture allowed us to interpret the evolutionary dynamics of spanin gene arrangement. Also, the potential universal role of intermolecular disulfide bonds in two-component spanin function was substantiated through bioinformatic and genetic approaches. Additionally, a novel lipobox motif, AWAC, was identified and experimentally verified. Conclusions: The findings from this bioinformatic approach gave us instructive insights into spanin function, evolution, domain organization and provide a platform for future spanin annotation, as well as biochemical and genetic experiments. They also establish that spanins, like viral membrane fusion proteins, adopt different strategies to achieve fusion of the inner and outer membranes. [ABSTRACT FROM AUTHOR]

Background: Protein disulfide isomerase (PDI) and PDI-like proteins contain thioredoxin domains that catalyze protein disulfide bond, inhibit aggregation of misfolded proteins, and function in isomerization during protein folding in endoplasmic reticulum and responses during abiotic stresses. Chinese cabbage is widely recognized as an economically important, nutritious vegetable, but its yield is severely hampered by various biotic and abiotic stresses. Because of, it is prime need to identify those genes whose are responsible for biotic and abiotic stress tolerance. PDI family genes are among of them. Results: We have identified 32 PDI genes from the Br135K microarray dataset, NCBI and BRAD database, and in silico characterized their sequences. Expression profiling of those genes was performed using cDNA of plant samples imposed to abiotic stresses; cold, salt, drought and ABA (Abscisic Acid) and biotic stress; Fusarium oxysporum f. sp. conglutinans infection. The Chinese cabbage PDI genes were clustered in eleven groups in phylogeny. Among them, 15 PDI genes were ubiquitously expressed in various organs, while 24 PDI genes were up-regulated under salt and drought stress. By contrast, cold and ABA stress responsive gene number were ten and nine, respectively. In case of F. oxysporum f. sp. conglutinans infection 14 BrPDI genes were highly up-regulated. Interestingly, BrPDI1-1 gene was identified as putative candidate against abiotic (salt and drought) and biotic stresses, BrPDI5-2 gene for ABA stress, and BrPDI1-4, 6-1 and 9-2 were putative candidate genes for both cold and chilling injury stresses. Conclusions: Our findings help to elucidate the involvement of PDI genes in stress responses, and they lay the foundation for functional genomics in future studies and molecular breeding of Brassica rapa crops. The stress-responsive PDI genes could be potential resources for molecular breeding of Brassica crops resistant to biotic and abiotic stresses. [ABSTRACT FROM AUTHOR]

Background: Identifying orthologous genes is an initial step required for phylogenetics, and it is also a common strategy employed in functional genetics to find candidates for functionally equivalent genes across multiple species. At the same time, in silico orthology prediction tools often require large computational resources only available on computing clusters. Here we present OrthoReD, an open-source orthology prediction tool with accuracy comparable to published tools that requires only a desktop computer. The low computational resource requirement of OrthoReD is achieved by repeating orthology searches on one gene of interest at a time, thereby generating a reduced dataset to limit the scope of orthology search for each gene of interest. Results: The output of OrthoReD was highly similar to the outputs of two other published orthology prediction tools, OrthologID and/or OrthoDB, for the three dataset tested, which represented three phyla with different ranges of species diversity and different number of genomes included. Median CPU time for ortholog prediction per gene by OrthoReD executed on a desktop computer was <15 min even for the largest dataset tested, which included all coding sequences of 100 bacterial species. Conclusions: With high-throughput sequencing, unprecedented numbers of genes from non-model organisms are available with increasing need for clear information about their orthologies and/or functional equivalents in model organisms. OrthoReD is not only fast and accurate as an orthology prediction tool, but also gives researchers flexibility in the number of genes analyzed at a time, without requiring a high-performance computing cluster. [ABSTRACT FROM AUTHOR]

Background: Sugarcane is an important sugar crop contributing up to about 80% of the world sugar production. Efforts to characterize the genes involved in sugar metabolism at the molecular level are growing since increasing sugar content is a major goal in the breeding of new sugarcane varieties. Fructokinases (FRK) are the main fructose phosphorylating enzymes with high substrate specificity and affinity. Results: In this study, by combining comparative genomics approaches with BAC resources, seven fructokinase genes were identified in S. spontaneum. Phylogenetic analysis based on representative monocotyledon and dicotyledon plant species suggested that the FRK gene family is ancient and its evolutionary history can be traced in duplicated descending order: SsFRK4, SsFRK6/SsFRK7,SsFRK5, SsFRK3 and SsFRK1/SsFRK2. Among the close orthologs, the number and position of exons in FRKs were conserved; in contrast, the size of introns varied among the paralogous FRKs in Saccharum. Genomic constraints were analyzed within the gene alleles and between S. spontaneum and Sorghum bicolor, and gene expression analysis was performed under drought stress and with exogenous applications of plant hormones. FRK1, which was under strong functional constraint selection, was conserved among the gene allelic haplotypes, and displayed dominant expression levels among the gene families in the control conditions, suggesting that FRK1 plays a major role in the phosphorylation of fructose. FRK3 and FRK5 were dramatically induced under drought stress, and FRK5 was also found to increase its expression levels in the mature stage of Saccharum. Similarly, FRK3 and FRK5 were induced in response to drought stress in Saccharum. FRK2 and FRK7 displayed lower expression levels than the other FRK family members; FRK2 was under strong genomic selection constraints whereas FRK7 was under neutral selection. FRK7 may have become functionally redundant in Saccharum through pseudogenization. FRK4 and FRK6 shared the most similar expression pattern: FRK4 was revealed to have higher expression levels in mature tissues than in premature tissues of Saccharum, and FRK6 presented a slight increase under drought stress. Conclusions: Our study presents a comprehensive genomic study of the entire FRK gene family in Saccharum, providing the foundations for approaches to characterize the molecular mechanism regulated by the SsFRK family in sugarcane. [ABSTRACT FROM AUTHOR]

Background: Body plan development in multi-cellular organisms is largely determined by homeotic genes. Expression of homeotic genes, in turn, is partially regulated by insulator binding proteins (IBPs). While only a few enhancer blocking IBPs have been identified in vertebrates, the common fruit fly Drosophila melanogaster harbors at least twelve different enhancer blocking IBPs. We screened recently compiled insect transcriptomes from the 1KITE project and genomic and transcriptomic data from public databases, aiming to trace the origin of IBPs in insects and other arthropods. Results: Our study shows that the last common ancestor of insects (Hexapoda) already possessed a substantial number of IBPs. Specifically, of the known twelve insect IBPs, at least three (i.e., CP190, Su(Hw), and CTCF) already existed prior to the evolution of insects. Furthermore we found GAF orthologs in early branching insect orders, including Zygentoma (silverfish and firebrats) and Diplura (two-pronged bristletails). Mod(mdg4) is most likely a derived feature of Neoptera, while Pita is likely an evolutionary novelty of holometabolous insects. Zw5 appears to be restricted to schizophoran flies, whereas BEAF-32, ZIPIC and the Elba complex, are probably unique to the genus Drosophila. Selection models indicate that insect IBPs evolved under neutral or purifying selection. Conclusions: Our results suggest that a substantial number of IBPs either pre-date the evolution of insects or evolved early during insect evolution. This suggests an evolutionary history of insulator binding proteins in insects different to that previously thought. Moreover, our study demonstrates the versatility of the 1KITE transcriptomic data for comparative analyses in insects and other arthropods. [ABSTRACT FROM AUTHOR]

Background: Feline immunodeficiency virus (FIV) is a global pathogen of Felidae species and a model system for Human immunodeficiency virus (HIV)-induced AIDS. In felids such as the domestic cat (Felis catus), APOBEC3 (A3) genes encode for single-domain A3Z2s, A3Z3 and double-domain A3Z2Z3 anti-viral cytidine deaminases. The feline A3Z2Z3 is expressed following read-through transcription and alternative splicing, introducing a previously untranslated exon in frame, encoding a domain insertion called linker. Only A3Z3 and A3Z2Z3 inhibit Vif-deficient FIV. Feline A3s also are restriction factors for HIV and Simian immunodeficiency viruses (SIV). Surprisingly, HIV-2/SIV Vifs can counteract feline A3Z2Z3. Results: To identify residues in feline A3s that Vifs need for interaction and degradation, chimeric human-feline A3s were tested. Here we describe the molecular direct interaction of feline A3s with Vif proteins from cat FIV and present the first structural A3 model locating these interaction regions. In the Z3 domain we have identified residues involved in binding of FIV Vif, and their mutation blocked Vif-induced A3Z3 degradation. We further identified additional essential residues for FIV Vif interaction in the A3Z2 domain, allowing the generation of FIV Vif resistant A3Z2Z3. Mutated feline A3s also showed resistance to the Vif of a lion-specific FIV, indicating an evolutionary conserved Vif-A3 binding. Comparative modelling of feline A3Z2Z3 suggests that the residues interacting with FIV Vif have, unlike Vif-interacting residues in human A3s, a unique location at the domain interface of Z2 and Z3 and that the linker forms a homeoboxlike domain protruding of the Z2Z3 core. HIV-2/SIV Vifs efficiently degrade feline A3Z2Z3 by possible targeting the linker stretch connecting both Z-domains. Conclusions: Here we identified in feline A3s residues important for binding of FIV Vif and a unique protein domain insertion (linker). To understand Vif evolution, a structural model of the feline A3 was developed. Our results show that HIV Vif binds human A3s differently than FIV Vif feline A3s. The linker insertion is suggested to form a homeo-box domain, which is unique to A3s of cats and related species, and not found in human and mouse A3s. Together, these findings indicate a specific and different A3 evolution in cats and human. [ABSTRACT FROM AUTHOR]

Background: Given a gene and a species tree, reconciliation methods attempt to retrieve the macro-evolutionary events that best explain the discrepancies between the two tree topologies. The DTL parsimonious approach searches for a most parsimonious reconciliation between a gene tree and a (dated) species tree, considering four possible macro-evolutionary events (speciation, duplication, transfer, and loss) with specific costs. Unfortunately, many events are erroneously predicted due to errors in the input trees, inappropriate input cost values or because of the existence of several equally parsimonious scenarios. It is thus crucial to provide a measure of the reliability for predicted events. It has been recently proposed that the reliability of an event can be estimated via its frequency in the set of most parsimonious reconciliations obtained using a variety of reasonable input cost vectors. To compute such a support, a straightforward but time-consuming approach is to generate the costs slightly departing from the original ones, independently compute the set of all most parsimonious reconciliations for each vector, and combine these sets a posteriori. Another proposed approach uses Pareto-optimality to partition cost values into regions which induce reconciliations with the same number of DTL events. The support of an event is then defined as its frequency in the set of regions. However, often, the number of regions is not large enough to provide reliable supports. Results: We present here a method to compute efficiently event supports via a polynomial-sized graph, which can represent all reconciliations for several different costs. Moreover, two methods are proposed to take into account alternative input costs: either explicitly providing an input cost range or allowing a tolerance for the over cost of a reconciliation. Our methods are faster than the region based method, substantially faster than the sampling-costs approach, and have a higher event-prediction accuracy on simulated data. Conclusions: We propose a new approach to improve the accuracy of event supports for parsimonious reconciliation methods to account for uncertainty in the input costs. Furthermore, because of their speed, our methods can be used on large gene families. Our algorithms are implemented in the ecceTERA program, freely available from http://mbb.univ-montp2.fr/MBB/. [ABSTRACT FROM AUTHOR]

Background Methionine adenosyltransferase (MAT) is a ubiquitous essential enzyme that, in eukaryotes, occurs in two relatively divergent paralogues: MAT and MATX. MATX has a punctate distribution across the tree of eukaryotes and, except for a few cases, is mutually exclusive with MAT. This phylogenetic pattern could have arisen by either differential loss of old paralogues or the spread of one of these paralogues by horizontal gene transfer. Our aim was to map the distribution of MAT/MATX genes within the Euglenida in order to more comprehensively characterize the evolutionary history of MATX. Results We generated 26 new sequences from 23 different lineages of euglenids and one prasinophyte alga Pyramimonas parkeae. MATX was present only in photoautotrophic euglenids. The mixotroph Rapaza viridis and the prasinophyte alga Pyramimonas parkeae, which harbors chloroplasts that are most closely related to the chloroplasts in photoautotrophic euglenids, both possessed only the MAT paralogue. We found both the MAT and MATX paralogues in two photoautotrophic species (Phacus orbicularis and Monomorphina pyrum). The significant conflict between eukaryotic phylogenies inferred from MATX and SSU rDNA data represents strong evidence that MATX paralogues have undergone horizontal gene transfer across the tree of eukaryotes. Conclusions Our results suggest that MATX entered the euglenid lineage in a single horizontal gene transfer event that took place after the secondary endosymbiotic origin of the euglenid chloroplast. The origin of the MATX paralogue is unclear, and it cannot be excluded that it arose by a gene duplication event before the most recent common ancestor of eukaryotes. [ABSTRACT FROM AUTHOR]

Background Hydrophobins are small secreted cysteine-rich proteins that play diverse roles during different phases of fungal life cycle. In basidiomycetes, hydrophobin-encoding genes often form large multigene families with up to 40 members. The evolutionary forces driving hydrophobin gene expansion and diversification in basidiomycetes are poorly understood. The functional roles of individual genes within such gene families also remain unclear. The relationship between the hydrophobin gene number, the genome size and the lifestyle of respective fungal species has not yet been thoroughly investigated. Here, we present results of our survey of hydrophobin gene families in two species of wood-degrading basidiomycetes, Phlebia brevispora and Heterobasidion annosum s.l. We have also investigated the regulatory pattern of hydrophobin-encoding genes from H. annosum s.s. during saprotrophic growth on pine wood as well as on culture filtrate from Phlebiopsis gigantea using micro-arrays. These data are supplemented by results of the protein structure modeling for a representative set of hydrophobins. Results We have identified hydrophobin genes from the genomes of two wood-degrading species of basidiomycetes, Heterobasidion irregulare, representing one of the microspecies within the aggregate H. annosum s.l., and Phlebia brevispora. Although a high number of hydrophobinencoding genes were observed in H. irregulare (16 copies), a remarkable expansion of these genes was recorded in P. brevispora (26 copies). A significant expansion of hydrophobinencoding genes in other analyzed basidiomycetes was also documented (1–40 copies), whereas contraction through gene loss was observed among the analyzed ascomycetes (1–11 copies). Our phylogenetic analysis confirmed the important role of gene duplication events in the evolution of hydrophobins in basidiomycetes. Increased number of hydrophobin-encoding genes appears to have been linked to the species' ecological strategy, with the non-pathogenic fungi having increased numbers of hydrophobins compared with their pathogenic counterparts. However, there was no significant relationship between the number of hydrophobin-encoding genes and genome size. Furthermore, our results revealed significant differences in the expression levels of the 16 H. annosum s.s.hydrophobin-encoding genes which suggest possible differences in their regulatory patterns. Conclusions A considerable expansion of the hydrophobin-encoding genes in basidiomycetes has been observed. The distribution and number of hydrophobin-encoding genes in the analyzed species may be connected to their ecological preferences. Results of our analysis also have shown that H. annosum s.l hydrophobin-encoding genes may be under positive selection. Our gene expression analysis revealed differential expression of H. annosum s.s. hydrophobin genes under different growth conditions, indicating their possible functional diversification. [ABSTRACT FROM AUTHOR]

Background: Staphylococcus aureus is a widely distributed human pathogen capable of infecting almost every ecological niche of the host. As a result, it is responsible for causing many different diseases. S. aureus has a vast array of virulence determinants whose expression is modulated by an intricate regulatory network, where transcriptional factors (TFs) are the primary elements. In this work, using diverse sequence analysis, we evaluated the repertoire of TFs and sigma factors in the community-associated methicillin resistant S. aureus (CA-MRSA) strain USA300-FPR3757. Results: A total of 135 TFs and sigma factors were identified and classified into 36 regulatory families. From these around 43% have been experimentally characterized to date, which demonstrates the significant work still at hand to unravel the regulatory network in place for this important pathogen. A comparison of the TF repertoire of S. aureus against 1209 sequenced bacterial genomes was carried out allowing us to identify a core set of orthologous TFs for the Staphylococacceae, and also allowing us to assign potential functions to previously uncharacterized TFs. Finally, the USA300 TFs were compared to those in eleven other S. aureus strains including: Newman, COL, JH1, JH9, MW2, Mu3, Mu50, N315, RF122, MRSA252 and MSSA476. We identify conserved TFs among these strains and suggest possible regulatory interactions. Conclusions: The analysis presented herein highlights the complexity of regulatory networks in S. aureus strains, identifies key conserved TFs among the Staphylococacceae, and offers unique insights into several as yet uncharacterized TFs. [ABSTRACT FROM AUTHOR]

Background: Rapeseed is the third largest oil seed crop in the world. The seeds of this plant store lipids in oil bodies, and oleosin is the most important structural protein in oil bodies. However, the function of oleosin in oil crops has received little attention. Results: In the present study, 48 oleosin sequences from the Brassica napus genome were identified and divided into four lineages (T, U, SH, SL). Synteny analysis revealed that most of the oleosin genes were conserved, and all of these genes experienced purifying selection during evolution. Three and four important oleosin genes from Arabidopsis and B. napus, respectively, were cloned and analyzed for function in Arabidopsis. Overexpression of these oleosin genes in Arabidopsis increased the seed oil content slightly, except for BnaOLE3. Further analysis revealed that the average oil body size of the transgenic seeds was slightly larger than that of the wild type (WT), except for BnaOLE1. The fatty acid profiles showed that the linoleic acid content (13.3% at most) increased and the peanut acid content (11% at most) decreased in the transgenic lines. In addition, the seed size and thousand-seed weight (TSW) also increased in the transgenic lines, which could lead to increased total lipid production. Conclusion: We identified oleosin genes in the B. napus genome, and overexpression of oleosin in Arabidopsis seeds increased the seed weight and linoleic acid content (13.3% at most). [ABSTRACT FROM AUTHOR]

Background: The SWEET (Sugars Will Eventually be Exported Transporters) gene family is a recently identified group of sugar transporters that play an indispensable role in sugar efflux, phloem loading, plant-pathogen interaction, nectar secretion, and reproductive tissue development. However, little information on Saccharum SWEET is available for this crop with a complex genetic background. Results: In this study, 22 SWEET genes were identified from Saccharum spontaneum Bacterial Artificial Chromosome libraries sequences. Phylogenetic analyses of SWEETs from 11 representative plant species showed that gene expansions of the SWEET family were mainly caused by the recent gene duplication in dicot plants, while these gene expansions were attributed to the ancient whole genome duplication (WGD) in monocot plant species. Gene expression profiles were obtained from RNA-seq analysis. SWEET1a and SWEET2s had higher expression levels in the transitional zone and maturing zone than in the other analyzed zones. SWEET1b was mainly expressed in the leaf tissues and the mature zone of the leaf of both S. spontaneum and S. officinarum, and displayed a peak in the morning and was undetectable in both sclerenchyma and parenchyma cells from the mature stalks of S. officinarum. SsSWEET4a\4b had higher expression levels than SWEET4c and were mainly expressed in the stems of seedlings and mature plants. SWEET13s are recently duplicated genes, and the expression of SWEET13s dramatically increased from the maturing to mature zones. SWEET16b's expression was not detected in S. officinarum, but displayed a rhythmic diurnal expression pattern. Conclusions: Our study revealed the gene evolutionary history of SWEETs in Saccharum and SWEET1b was found to be a sucrose starvation-induced gene involved in the sugar transportation in the high photosynthetic zones. SWEET13c was identified as the key player in the efflux of sugar transportation in mature photosynthetic tissues. SWEET4a\4b were found to be mainly involved in sugar transportation in the stalk. SWEET1a\2a\4a\4b\13a\16b were suggested to be the genes contributing to the differences in sugar contents between S. spontaneum and S. officinarum. Our results are valuable for further functional analysis of SWEET genes and utilization of the SWEET genes for genetic improvement of Saccharum for biofuel production. [ABSTRACT FROM AUTHOR]

Background: Protein disulfide isomerase (PDI) and PDI-like proteins contain thioredoxin domains that catalyze protein disulfide bond, inhibit aggregation of misfolded proteins, and function in isomerization during protein folding in endoplasmic reticulum and responses during abiotic stresses.Chinese cabbage is widely recognized as an economically important, nutritious vegetable, but its yield is severely hampered by various biotic and abiotic stresses. Because of, it is prime need to identify those genes whose are responsible for biotic and abiotic stress tolerance. PDI family genes are among of them., Results: We have identified 32 PDI genes from the Br135K microarray dataset, NCBI and BRAD database, and in silico characterized their sequences. Expression profiling of those genes was performed using cDNA of plant samples imposed to abiotic stresses; cold, salt, drought and ABA (Abscisic Acid) and biotic stress; Fusarium oxysporum f. sp. conglutinans infection. The Chinese cabbage PDI genes were clustered in eleven groups in phylogeny. Among them, 15 PDI genes were ubiquitously expressed in various organs, while 24 PDI genes were up-regulated under salt and drought stress. By contrast, cold and ABA stress responsive gene number were ten and nine, respectively. In case of F. oxysporum f. sp. conglutinans infection 14 BrPDI genes were highly up-regulated. Interestingly, BrPDI1-1 gene was identified as putative candidate against abiotic (salt and drought) and biotic stresses, BrPDI5-2 gene for ABA stress, and BrPDI1-4, 6-1 and 9-2 were putative candidate genes for both cold and chilling injury stresses., Conclusions: Our findings help to elucidate the involvement of PDI genes in stress responses, and they lay the foundation for functional genomics in future studies and molecular breeding of Brassica rapa crops. The stress-responsive PDI genes could be potential resources for molecular breeding of Brassica crops resistant to biotic and abiotic stresses.

Background: Sugarcane is an important sugar crop contributing up to about 80% of the world sugar production. Efforts to characterize the genes involved in sugar metabolism at the molecular level are growing since increasing sugar content is a major goal in the breeding of new sugarcane varieties. Fructokinases (FRK) are the main fructose phosphorylating enzymes with high substrate specificity and affinity., Results: In this study, by combining comparative genomics approaches with BAC resources, seven fructokinase genes were identified in S. spontaneum. Phylogenetic analysis based on representative monocotyledon and dicotyledon plant species suggested that the FRK gene family is ancient and its evolutionary history can be traced in duplicated descending order: SsFRK4, SsFRK6/SsFRK7,SsFRK5, SsFRK3 and SsFRK1/SsFRK2. Among the close orthologs, the number and position of exons in FRKs were conserved; in contrast, the size of introns varied among the paralogous FRKs in Saccharum. Genomic constraints were analyzed within the gene alleles and between S. spontaneum and Sorghum bicolor, and gene expression analysis was performed under drought stress and with exogenous applications of plant hormones. FRK1, which was under strong functional constraint selection, was conserved among the gene allelic haplotypes, and displayed dominant expression levels among the gene families in the control conditions, suggesting that FRK1 plays a major role in the phosphorylation of fructose. FRK3 and FRK5 were dramatically induced under drought stress, and FRK5 was also found to increase its expression levels in the mature stage of Saccharum. Similarly, FRK3 and FRK5 were induced in response to drought stress in Saccharum. FRK2 and FRK7 displayed lower expression levels than the other FRK family members; FRK2 was under strong genomic selection constraints whereas FRK7 was under neutral selection. FRK7 may have become functionally redundant in Saccharum through pseudogenization. FRK4 and FRK6 shared the most similar expression pattern: FRK4 was revealed to have higher expression levels in mature tissues than in premature tissues of Saccharum, and FRK6 presented a slight increase under drought stress., Conclusions: Our study presents a comprehensive genomic study of the entire FRK gene family in Saccharum, providing the foundations for approaches to characterize the molecular mechanism regulated by the SsFRK family in sugarcane.

Background: Feline immunodeficiency virus (FIV) is a global pathogen of Felidae species and a model system for Human immunodeficiency virus (HIV)-induced AIDS. In felids such as the domestic cat (Felis catus), APOBEC3 (A3) genes encode for single-domain A3Z2s, A3Z3 and double-domain A3Z2Z3 anti-viral cytidine deaminases. The feline A3Z2Z3 is expressed following read-through transcription and alternative splicing, introducing a previously untranslated exon in frame, encoding a domain insertion called linker. Only A3Z3 and A3Z2Z3 inhibit Vif-deficient FIV. Feline A3s also are restriction factors for HIV and Simian immunodeficiency viruses (SIV). Surprisingly, HIV-2/SIV Vifs can counteract feline A3Z2Z3., Results: To identify residues in feline A3s that Vifs need for interaction and degradation, chimeric human-feline A3s were tested. Here we describe the molecular direct interaction of feline A3s with Vif proteins from cat FIV and present the first structural A3 model locating these interaction regions. In the Z3 domain we have identified residues involved in binding of FIV Vif, and their mutation blocked Vif-induced A3Z3 degradation. We further identified additional essential residues for FIV Vif interaction in the A3Z2 domain, allowing the generation of FIV Vif resistant A3Z2Z3. Mutated feline A3s also showed resistance to the Vif of a lion-specific FIV, indicating an evolutionary conserved Vif-A3 binding. Comparative modelling of feline A3Z2Z3 suggests that the residues interacting with FIV Vif have, unlike Vif-interacting residues in human A3s, a unique location at the domain interface of Z2 and Z3 and that the linker forms a homeobox-like domain protruding of the Z2Z3 core. HIV-2/SIV Vifs efficiently degrade feline A3Z2Z3 by possible targeting the linker stretch connecting both Z-domains., Conclusions: Here we identified in feline A3s residues important for binding of FIV Vif and a unique protein domain insertion (linker). To understand Vif evolution, a structural model of the feline A3 was developed. Our results show that HIV Vif binds human A3s differently than FIV Vif feline A3s. The linker insertion is suggested to form a homeo-box domain, which is unique to A3s of cats and related species, and not found in human and mouse A3s. Together, these findings indicate a specific and different A3 evolution in cats and human.

Background: Studying gene evolution in non-model species by PCR-based approaches is limited to highly conserved genes. The plummeting cost of next generation sequencing enables the application of de novo transcriptomics to any species., Results: Here we describe how to apply de novo transcriptomics to pursue the evolution of a single gene of interest. We follow a rapidly evolving seed protein that encodes small, stable peptides. We use software that needs limited bioinformatics background and assemble four de novo seed transcriptomes. To demonstrate the quality of the assemblies, we confirm the predicted genes at the peptide level on one species which has over ten copies of our gene of interest. We explain strategies that favour assembly of low abundance genes, what assembly parameters help capture the maximum number of transcripts, how to develop a suite of control genes to test assembly quality and we compare several sequence depths to optimise cost and data volume., Conclusions: De novo transcriptomics is an effective approach for studying gene evolution in species for which genome support is lacking.