Conserved Coding Sequences Research Articles

An Integrated Approach for Microprotein Identification and Sequence Analysis.

Next-generation sequencing (NGS) has propelled the field of genomics forward and produced whole genome sequences for numerous animal species and model organisms. However, despite this wealth of sequence information, comprehensive gene annotation efforts have proven challenging, especially for small proteins. Notably, conventional protein annotation methods were designed to intentionally exclude putative proteins encoded by short open reading frames (sORFs) less than 300 nucleotides in length to filter out the exponentially higher number of spurious noncoding sORFs throughout the genome. As a result, hundreds of functional small proteins called microproteins (<100 amino acids in length) have been incorrectly classified as noncoding RNAs or overlooked entirely. Here we provide a detailed protocol to leverage free, publicly available bioinformatic tools to query genomic regions for microprotein-coding potential based on evolutionary conservation. Specifically, we provide step-by-step instructions on how to examine sequence conservation and coding potential using Phylogenetic Codon Substitution Frequencies (PhyloCSF) on the user-friendly University of California Santa Cruz (UCSC) Genome Browser. Additionally, we detail steps to efficiently generate multiple species alignments of identified microprotein sequences to visualize amino acid sequence conservation and recommend resources to analyze microprotein characteristics, including predicted domain structures. These powerful tools can be used to help identify putative microprotein-coding sequences in noncanonical genomic regions or to rule out the presence of a conserved coding sequence with translational potential in a noncoding transcript of interest.

Functional characterization of Melampsora larici‐populina Hog1‐type MAPK in response to osmotic stress

AbstractThe obligate biotrophic pathogen Melampsora larici‐populina, responsible for poplar foliar rust disease, causes annual epidemics and severe damages to poplar plantations worldwide. Flexible responses to external osmotic pressure changes are important for the growth and survival of pathogens. In a previous study, we suggested that the Hog1‐type MAPK MlpHog1 in M. larici‐populina may play a role in infectious growth and responses to various environmental stresses. In the present study, we analysed its biological characteristics. MlpHog1 displayed a conserved coding sequence pattern among five strains from different regions of China. MlpHog1 restored the Hog1‐orthologue mutant defects responding to hyperosmotic stress in both Saccharomyces cerevisiae and Magnaporthe oryzae. Transient expression in wheat protoplasts revealed that MlpHog1 is localized in the cytoplasm and nucleus. These results indicate that MlpHog1 plays a positive role in responding to osmotic stress in the poplar rust M. larici‐populina.

The foraging Gene Is Involved in the Presence of Wings and Explorative Behaviours in Parthenogenetic Females of the Aphid Myzus persicae.

The foraging gene (for) encodes for a cyclic guanosine monophosphate (cGMP)-dependent protein kinase involved in behavioural plasticity in aphids and in other insects. In this paper, we analysed the complete for sequence in eight clones of the peach potato aphid Myzus persicae, reporting the presence of nonsense and frameshift mutations in three studied clones characterized by a reduced number of winged females and by the absence of exploratory behaviours. Quantitative PCR experiments evidenced similar results in clones possessing for genes with a conserved coding sequence, but low expression levels. The comparison of the for transcriptional level in Myzus persicae persicae and Myzus persicae nicotianae showed very different expression in the two studied M. p. nicotianae clones so that our data did not support a previous hypothesis suggesting that a differential for expression was related to ecological specialization of M. p. nicotianae. In view of its role in both the dispersal of winged females and exploratory behaviours, the screening of the for sequences could be useful for predicting invasions of cultivated areas by peach potato aphids.

Characterization of two newly isolated bacteriophages PW2 and PW4 and derived endolysins with lysis activity against Bacillus cereus group strains

This study characterized two novel Siphoviridae phages, PW2 and PW4, that can infect 52% and 44% of the tested Bacillus cereus group isolates and display relatively high activity against four cereulide-producing isolates belonging to B. weihenstephanensis and B. paranthracis. The genome sequences of PW2 and PW4 are similar to six known phages infecting B. cereus group isolates, which can be classified into two conserved groups, with the PW2 genome harboring conserved coding sequences (CDSs) from both groups. Two phage-derived endolysins, LysPW2 and LysPW4, which are predicted to encode N-acetylmuramoyl-L-alanine amidase, and their enzymatically active domains (EADs), LysPW2-EAD and LysPW4-EAD, were heterologously expressed. Both LysPW2 and LysPW4, especially the former, show a much wider host range than the phages, albeit still limited to the B. cereus group for the tested bacteria. The optimal temperature and pH for LysPW2 ability is 37 °C and pH 8.0 and for LysPW4 is 50 °C and pH 9.0. Neither LysPW2-EAD nor LysPW4-EAD show any lytic activity against vegetative cells of the tested B. cereus group isolates but can inhibit germination in 66.3% and 65.7% of spores, respectively. In addition, both LysPW2-EAD and LysPW4-EAD exhibit spore-binding capabilities.

TALEN-mediated targeted mutagenesis of fatty acid desaturase 2 (FAD2) in peanut (Arachis hypogaea L.) promotes the accumulation of oleic acid

A first creation of high oleic acid peanut varieties by using transcription activator-like effecter nucleases (TALENs) mediated targeted mutagenesis of Fatty Acid Desaturase 2 (FAD2). Transcription activator like effector nucleases (TALENs), which allow the precise editing of DNA, have already been developed and applied for genome engineering in diverse organisms. However, they are scarcely used in higher plant study and crop improvement, especially in allopolyploid plants. In the present study, we aimed to create targeted mutagenesis by TALENs in peanut. Targeted mutations in the conserved coding sequence of Arachis hypogaea fatty acid desaturase 2 (AhFAD2) were created by TALENs. Genetic stability of AhFAD2 mutations was identified by DNA sequencing in up to 9.52 and 4.11% of the regeneration plants at two different targeted sites, respectively. Mutation frequencies among AhFAD2 mutant lines were significantly correlated to oleic acid accumulation. Genetically, stable individuals of positive mutant lines displayed a 0.5-2 fold increase in the oleic acid content compared with non-transgenic controls. This finding suggested that TALEN-mediated targeted mutagenesis could increase the oleic acid content in edible peanut oil. Furthermore, this was the first report on peanut genome editing event, and the obtained high oleic mutants could serve for peanut breeding project.

Development of single-copy nuclear intron markers for species-level phylogenetics: Case study with Paullinieae (Sapindaceae).

Premise of the study:We developed a bioinformatic pipeline that leverages a publicly available genome and published transcriptomes to design primers in conserved coding sequences flanking targeted introns of single-copy nuclear loci. Paullinieae (Sapindaceae) is used to demonstrate the pipeline.Methods and Results:Transcriptome reads phylogenetically closer to the lineage of interest are aligned to the closest genome. Single-nucleotide polymorphisms are called, generating a “pseudoreference” closer to the lineage of interest. Several filters are applied to meet the criteria of single-copy nuclear loci with introns of a desired size. Primers are designed in conserved coding sequences flanking introns. Using this pipeline, we developed nine single-copy nuclear intron markers for Paullinieae.Conclusions:This pipeline is highly flexible and can be used for any group with available genomic and transcriptomic resources. This pipeline led to the development of nine variable markers for phylogenetic study without generating sequence data de novo.

Genotypic and phenotypic characterization of multidrug resistant Salmonella Typhimurium and Salmonella Kentucky strains recovered from chicken carcasses.

Salmonella Typhimurium is the leading cause of human non-typhoidal gastroenteritis in the US. S. Kentucky is one the most commonly recovered serovars from commercially processed poultry carcasses. This study compared the genotypic and phenotypic properties of two Salmonella enterica strains Typhimurium (ST221_31B) and Kentucky (SK222_32B) recovered from commercially processed chicken carcasses using whole genome sequencing, phenotype characterizations and an intracellular killing assay. Illumina MiSeq platform was used for sequencing of two Salmonella genomes. Phylogenetic analysis employing homologous alignment of a 1,185 non-duplicated protein-coding gene in the Salmonella core genome demonstrated fully resolved bifurcating patterns with varying levels of diversity that separated ST221_31B and SK222_32B genomes into distinct monophyletic serovar clades. Single nucleotide polymorphism (SNP) analysis identified 2,432 (ST19) SNPs within 13 Typhimurium genomes including ST221_31B representing Sequence Type ST19 and 650 (ST152) SNPs were detected within 13 Kentucky genomes including SK222_32B representing Sequence Type ST152. In addition to serovar-specific conserved coding sequences, the genomes of ST221_31B and SK222_32B harbor several genomic regions with significant genetic differences. These included phage and phage-like elements, carbon utilization or transport operons, fimbriae operons, putative membrane associated protein-encoding genes, antibiotic resistance genes, siderophore operons, and numerous hypothetical protein-encoding genes. Phenotype microarray results demonstrated that ST221_31B is capable of utilizing certain carbon compounds more efficiently as compared to SK222_3B; namely, 1,2-propanediol, M-inositol, L-threonine, α-D-lactose, D-tagatose, adonitol, formic acid, acetoacetic acid, and L-tartaric acid. ST221_31B survived for 48 h in macrophages, while SK222_32B was mostly eliminated. Further, a 3-fold growth of ST221_31B was observed at 24 hours post-infection in chicken granulosa cells while SK222_32B was unable to replicate in these cells. These results suggest that Salmonella Typhimurium can survive host defenses better and could be more invasive than Salmonella Kentucky and provide some insights into the genomic determinants responsible for these differences.

Sperm Bindin Divergence under Sexual Selection and Concerted Evolution in Sea Stars.

Selection associated with competition among males or sexual conflict between mates can create positive selection for high rates of molecular evolution of gamete recognition genes and lead to reproductive isolation between species. We analyzed coding sequence and repetitive domain variation in the gene encoding the sperm acrosomal protein bindin in 13 diverse sea star species. We found that bindin has a conserved coding sequence domain structure in all 13 species, with several repeated motifs in a large central region that is similar among all sea stars in organization but highly divergent among genera in nucleotide and predicted amino acid sequence. More bindin codons and lineages showed positive selection for high relative rates of amino acid substitution in genera with gonochoric outcrossing adults (and greater expected strength of sexual selection) than in selfing hermaphrodites. That difference is consistent with the expectation that selfing (a highly derived mating system) may moderate the strength of sexual selection and limit the accumulation of bindin amino acid differences. The results implicate both positive selection on single codons and concerted evolution within the repetitive region in bindin divergence, and suggest that both single amino acid differences and repeat differences may affect sperm-egg binding and reproductive compatibility.

Comprehensive cytological characterization of the Gossypium hirsutum genome based on the development of a set of chromosome cytological markers

Cotton is the world's most important natural fiber crop. It is also a model system for studying polyploidization, genomic organization, and genome-size variation. Integrating the cytological characterization of cotton with its genetic map will be essential for understanding its genome structure and evolution, as well as for performing further genetic-map based mapping and cloning. In this study, we isolated a complete set of bacterial artificial chromosome clones anchored to each of the 52 chromosome arms of the tetraploid cotton Gossypium hirsutum. Combining these with telomere and centromere markers, we constructed a standard karyotype for the G. hirsutum inbred line TM-1. We dissected the chromosome arm localizations of the 45S and 5S rDNA and suggest a centromere repositioning event in the homoeologous chromosomes AT09 and DT09. By integrating a systematic karyotype analysis with the genetic linkage map, we observed different genome sizes and chromosomal structures between the subgenomes of the tetraploid cotton and those of its diploid ancestors. Using evidence of conserved coding sequences, we suggest that the different evolutionary paths of non-coding retrotransposons account for most of the variation in size between the subgenomes of tetraploid cotton and its diploid ancestors. These results provide insights into the cotton genome and will facilitate further genome studies in G. hirsutum.

Genome Sequences of Three Strains of the Pseudomonas aeruginosa PA7 Clade.

Draft genome sequences of three P. aeruginosa strains from the PA7 clade are presented here. Their lengths are 6.36 (EML528), 6.44 (EML545), and 6.33 Mb (EML548). Comparisons with the PA7 genome showed 5,113 conserved coding sequences (CDSs), and significant numbers of strain-specific CDSs. Their analysis will improve our understanding of this highly divergent clade.

Gene expression analysis of E. coli strains provides insights into the role of gene regulation in diversification.

Escherichia coli spans a genetic continuum from enteric strains to several phylogenetically distinct, atypical lineages that are rare in humans, but more common in extra-intestinal environments. To investigate the link between gene regulation, phylogeny and diversification in this species, we analyzed global gene expression profiles of four strains representing distinct evolutionary lineages, including a well-studied laboratory strain, a typical commensal (enteric) strain and two environmental strains. RNA-Seq was employed to compare the whole transcriptomes of strains grown under batch, chemostat and starvation conditions. Highly differentially expressed genes showed a significantly lower nucleotide sequence identity compared with other genes, indicating that gene regulation and coding sequence conservation are directly connected. Overall, distances between the strains based on gene expression profiles were largely dependent on the culture condition and did not reflect phylogenetic relatedness. Expression differences of commonly shared genes (all four strains) and E. coli core genes were consistently smaller between strains characterized by more similar primary habitats. For instance, environmental strains exhibited increased expression of stress defense genes under carbon-limited growth and entered a more pronounced survival-like phenotype during starvation compared with other strains, which stayed more alert for substrate scavenging and catabolism during no-growth conditions. Since those environmental strains show similar genetic distance to each other and to the other two strains, these findings cannot be simply attributed to genetic relatedness but suggest physiological adaptations. Our study provides new insights into ecologically relevant gene-expression and underscores the role of (differential) gene regulation for the diversification of the model bacterial species.

Identification of small ORFs in vertebrates using ribosome footprinting and evolutionary conservation

Identification of the coding elements in the genome is a fundamental step to understanding the building blocks of living systems. Short peptides (< 100 aa) have emerged as important regulators of development and physiology, but their identification has been limited by their size. We have leveraged the periodicity of ribosome movement on the mRNA to define actively translated ORFs by ribosome footprinting. This approach identifies several hundred translated small ORFs in zebrafish and human. Computational prediction of small ORFs from codon conservation patterns corroborates and extends these findings and identifies conserved sequences in zebrafish and human, suggesting functional peptide products (micropeptides). These results identify micropeptide-encoding genes in vertebrates, providing an entry point to define their function in vivo.

Effect of duplicate genes on mouse genetic robustness: an update.

In contrast to S. cerevisiae and C. elegans, analyses based on the current knockout (KO) mouse phenotypes led to the conclusion that duplicate genes had almost no role in mouse genetic robustness. It has been suggested that the bias of mouse KO database toward ancient duplicates may possibly cause this knockout duplicate puzzle, that is, a very similar proportion of essential genes (P E) between duplicate genes and singletons. In this paper, we conducted an extensive and careful analysis for the mouse KO phenotype data and corroborated a strong effect of duplicate genes on mouse genetics robustness. Moreover, the effect of duplicate genes on mouse genetic robustness is duplication-age dependent, which holds after ruling out the potential confounding effect from coding-sequence conservation, protein-protein connectivity, functional bias, or the bias of duplicates generated by whole genome duplication (WGD). Our findings suggest that two factors, the sampling bias toward ancient duplicates and very ancient duplicates with a proportion of essential genes higher than that of singletons, have caused the mouse knockout duplicate puzzle; meanwhile, the effect of genetic buffering may be correlated with sequence conservation as well as protein-protein interactivity.

In silico identification of conserved intercoding sequences in Leishmania genomes: Unraveling putative cis-regulatory elements

In silico analyses of Leishmania spp. genome data are a powerful resource to improve the understanding of these pathogens’ biology. Trypanosomatids such as Leishmania spp. have their protein-coding genes grouped in long polycistronic units of functionally unrelated genes. The control of gene expression happens by a variety of posttranscriptional mechanisms. The high degree of synteny among Leishmania species is accompanied by highly conserved coding sequences (CDS) and poorly conserved intercoding untranslated sequences. To identify the elements involved in the control of gene expression, we conducted an in silico investigation to find conserved intercoding sequences (CICS) in the genomes of L. major, L. infantum, and L. braziliensis.We used a combination of computational tools, such as Linux-Shell, PERL and R languages, BLAST, MSPcrunch, SSAKE, and Pred-A-Term algorithms to construct a pipeline which was able to: (i) search for conservation in target-regions, (ii) eliminate CICS redundancy and mask repeat elements, (iii) predict the mRNA's extremities, (iv) analyze the distribution of orthologous genes within the generated LeishCICS-clusters, (v) assign GO terms to the LeishCICS-clusters, and (vi) provide statistical support for the gene-enrichment annotation. We associated the LeishCICS-cluster data, generated at the end of the pipeline, with the expression profile of L. donovani genes during promastigote–amastigote differentiation, as previously evaluated by others (GEO accession: GSE21936). A Pearson's correlation coefficient greater than 0.5 was observed for 730 LeishCICS-clusters containing from 2 to 17 genes. The designed computational pipeline is a useful tool and its application identified potential regulatory cis elements and putative regulons in Leishmania.

Exploration for Functional Nucleotide Sequence Candidates within Coding Regions of Mammalian Genes

The primary role of a protein coding gene is to encode amino acids. Therefore, synonymous sites of codons, which do not change the encoded amino acid, are regarded as evolving neutrally. However, if a certain region of a protein coding gene contains a functional nucleotide element (e.g. splicing signals), synonymous sites in the region may have selective pressure. The existence of such elements would be detected by searching regions of low nucleotide substitution. We explored invariant nucleotide sequences in 10 790 orthologous genes of six mammalian species (Homo sapiens, Macaca mulatta, Mus musculus, Rattus norvegicus, Bos taurus, and Canis familiaris), and extracted 4150 sequences whose conservation is significantly stronger than other regions of the gene and named them significantly conserved coding sequences (SCCSs). SCCSs are observed in 2273 genes. The genes are mainly involved with development, transcriptional regulation, and the neurons, and are expressed in the nervous system and the head and neck organs. No strong influence of conventional factors that affect synonymous substitution was observed in SCCSs. These results imply that SCCSs may have double function as nucleotide element and protein coding sequence and retained in the course of mammalian evolution.

Cloning, sequence identification and expression profile analysis of α-l-fucosidase gene from the Mediterranean fruit fly Ceratitis capitata

The Mediterranean fruit fly Ceratitis capitata (Diptera: Tephritidae) is one of the most destructive agricultural pests, a polyphagus insect of relevant economic importance and is widespread in many regions around the world. It is the best-studied fruit fly pest at genetic and molecular level and much has been learned on its ecology and behaviour. An α-l-fucosidase has been recently hypothesized to be involved in sperm–egg interactions in Drosophila melanogaster and in other Drosophila species. Here, a complete cDNA encoding a putative α-l-fucosidase of the medfly was amplified using the reverse polymerase chain reaction (RT-PCR) with degenerate based on the conserved coding sequence information of several insect α-l-fucosidases, cloned and sequenced (GenBank accession no. FJ177429). The coding region consisted of 1482bp which encoded a 485-residues protein (named CcFUCA) with a predicted molecular mass of 56.1kDa. The deduced protein sequence showed 75% amino acid identity to D. melanogaster α-l-fucosidase, and in fact the phylogenetic tree analysis revealed that CcFUCA had closer relationships with the α-l-fucosidases of drosophilid species. The tissue expression analysis indicated that CcFuca was expressed in a single transcript in all tissues, suggesting a ubiquitous localization pattern of the encoded protein. Our findings provide novel insights on a gene encoding a protein potentially involved in primary gamete interactions in C. capitata.

Some novel intron positions in conserved Drosophila genes are caused by intron sliding or tandem duplication

BackgroundPositions of spliceosomal introns are often conserved between remotely related genes. Introns that reside in non-conserved positions are either novel or remnants of frequent losses of introns in some evolutionary lineages. A recent gain of such introns is difficult to prove. However, introns verified as novel are needed to evaluate contemporary processes of intron gain.ResultsWe identified 25 unambiguous cases of novel intron positions in 31 Drosophila genes that exhibit near intron pairs (NIPs). Here, a NIP consists of an ancient and a novel intron position that are separated by less than 32 nt. Within a single gene, such closely-spaced introns are very unlikely to have coexisted. In most cases, therefore, the ancient intron position must have disappeared in favour of the novel one. A survey for NIPs among 12 Drosophila genomes identifies intron sliding (migration) as one of the more frequent causes of novel intron positions. Other novel introns seem to have been gained by regional tandem duplications of coding sequences containing a proto-splice site.ConclusionsRecent intron gains sometimes appear to have arisen by duplication of exonic sequences and subsequent intronization of one of the copies. Intron migration and exon duplication together may account for a significant amount of novel intron positions in conserved coding sequences.

Introducing the plant RNA editing prediction and analysis computer tool PREPACT and an update on RNA editing site nomenclature

Transcripts in mitochondria and chloroplasts of land plants are modified through RNA editing, the exchanges of pyrimidines-a post-transcriptional process that may affect more than 1,000 sites in the mitochondrial transcriptomes of some plant species. RNA editing mainly acts as a correcting mechanism to re-create evolutionary conserved coding sequences on mRNA level and can be reasonably well predicted in new plant organelle gene sequence data. Identification and annotation of RNA editing sites is cumbersome and error-prone for larger data sets or organelle sequences subject to highly frequent RNA editing. We here present PREPACT, WWW-accessible at http://www.prepact.de , which allows prediction, analysis, annotation and graphical display of RNA editing sites for both directions of pyrimidine exchanges, using the recently proposed RNA editing nomenclature. PREPACT offers prediction of RNA editing, analysis of partial editing in cDNA pools and a BLASTX mode for simultaneous prediction of genes and RNA editing sites in novel sequences. Output options include (i) lists and annotations of RNA editing sites, (ii) sequence alignments with user-controlled color highlighting of editings, (iii) graphical displays of RNA editing in sequences and alignments. Finally, binary matrices of editing positions can be produced that may be used for downstream (e.g. phylogenetic) analyses.

Accurate discrimination of conserved coding and non-coding regions through multiple indicators of evolutionary dynamics

BackgroundThe conservation of sequences between related genomes has long been recognised as an indication of functional significance and recognition of sequence homology is one of the principal approaches used in the annotation of newly sequenced genomes. In the context of recent findings that the number non-coding transcripts in higher organisms is likely to be much higher than previously imagined, discrimination between conserved coding and non-coding sequences is a topic of considerable interest. Additionally, it should be considered desirable to discriminate between coding and non-coding conserved sequences without recourse to the use of sequence similarity searches of protein databases as such approaches exclude the identification of novel conserved proteins without characterized homologs and may be influenced by the presence in databases of sequences which are erroneously annotated as coding.ResultsHere we present a machine learning-based approach for the discrimination of conserved coding sequences. Our method calculates various statistics related to the evolutionary dynamics of two aligned sequences. These features are considered by a Support Vector Machine which designates the alignment coding or non-coding with an associated probability score.ConclusionWe show that our approach is both sensitive and accurate with respect to comparable methods and illustrate several situations in which it may be applied, including the identification of conserved coding regions in genome sequences and the discrimination of coding from non-coding cDNA sequences.

Cryptic Sex in the Smallest Eukaryotic Marine Green Alga

Ostreococcus spp. are common worldwide oceanic picoeukaryotic pelagic algae. The complete genomes of three strains from different ecological niches revealed them to represent biologically distinct species despite their identical cellular morphologies (cryptic species). Their tiny genomes (13 Mb), with approximately 20 chromosomes, are colinear and densely packed with coding sequences, but no sexual life cycle has been described. Seventeen new strains of one of these species, Ostreococcus tauri, were isolated from 98 seawater samplings from the NW Mediterranean by filtering, culturing, cloning, and plating for single colonies and identification by sequencing their ribosomal 18S gene. In order to find the genetic markers for detection of polymorphisms and sexual recombination, we used an in silico approach to screen available genomic data. Intergenic regions of DNA likely to evolve neutrally were analyzed following polymerase chain reaction amplification of sequences using flanking primers from adjacent conserved coding sequences that were present as syntenic pairs in two different species of Ostreococcus. Analyses of such DNA regions from eight marker loci on two chromosomes from each strain revealed that the isolated O. tauri clones were haploid and that the overall level of polymorphism was approximately 0.01. Four different genetic tests for recombination showed that sexual exchanges must be inferred to account for the between-locus and between-chromosome marker combinations observed. However, our data suggest that sexual encounters are infrequent because we estimate the frequency of meioses/mitoses among the sampled strains to be 10(-6). Ostreococcus tauri and related species encode and express core genes for mitosis and meiosis, but their mechanisms of cell division and recombination, nevertheless, remain enigmatic because a classical eukaryotic spindle with 40 canonical microtubules would be much too large for the available approximately 0.9-microm(3) cellular volume.