Comparison of DNA sequences with protein sequences

Pairwise comparison of orthologous olfactory receptor genes between two sympatric sibling sea kraits of the genus Laticauda in Vanuatu

Olfaction-based reproductive isolation is widely observed in animals, but little is known about the genetic basis of such isolation mechanisms. Two species of sibling amphibious sea snakes, Laticauda colubrina and L. frontalis live in Vanuatu sympatrically and syntopically, but no natural hybrids have been reported. Adult females of both taxa possess distinctive lipids in the skin, and male L. frontalis distinguishes conspecific females based on olfactory cues. To shed light on the molecular basis of the evolution of olfaction-based isolation mechanisms, olfactory receptor (OR) gene repertoires of both taxa were identified using pyrosequencing-based technology, and orthologous OR gene sets were identified. Few species-specific gene duplications or species-specific gene losses were found. However, the nonsynonymous-to-synonymous substitution rate ratio was relatively higher between orthologous OR genes of L. frontalis and L. colubrina, indicating that L. frontalis and L. colubrina have evolved to possess different olfactory senses. We suggest that L. frontalis and L. colubrina have evolved allopatrically, and this may be a byproduct of the allopatric evolution, and that this dissimilarity may function as a premating isolation barrier, since L. frontalis has returned to the ancestral range (Vanuatu).

Identifying Mendelian disease genes with the variant effect scoring tool

Background

Whole exome sequencing studies identify hundreds to thousands of rare protein coding variants of ambiguous significance for human health. Computational tools are needed to accelerate the identification of specific variants and genes that contribute to human disease.

Results

We have developed the Variant Effect Scoring Tool (VEST), a supervised machine learning-based classifier, to prioritize rare missense variants with likely involvement in human disease. The VEST classifier training set comprised ~ 45,000 disease mutations from the latest Human Gene Mutation Database release and another ~45,000 high frequency (allele frequency >1%) putatively neutral missense variants from the Exome Sequencing Project. VEST outperforms some of the most popular methods for prioritizing missense variants in carefully designed holdout benchmarking experiments (VEST ROC AUC = 0.91, PolyPhen2 ROC AUC = 0.86, SIFT4.0 ROC AUC = 0.84). VEST estimates variant score p-values against a null distribution of VEST scores for neutral variants not included in the VEST training set. These p-values can be aggregated at the gene level across multiple disease exomes to rank genes for probable disease involvement. We tested the ability of an aggregate VEST gene score to identify candidate Mendelian disease genes, based on whole-exome sequencing of a small number of disease cases. We used whole-exome data for two Mendelian disorders for which the causal gene is known. Considering only genes that contained variants in all cases, the VEST gene score ranked dihydroorotate dehydrogenase (DHODH) number 2 of 2253 genes in four cases of Miller syndrome, and myosin-3 (MYH3) number 2 of 2313 genes in three cases of Freeman Sheldon syndrome.

Conclusions

Our results demonstrate the potential power gain of aggregating bioinformatics variant scores into gene-level scores and the general utility of bioinformatics in assisting the search for disease genes in large-scale exome sequencing studies. VEST is available as a stand-alone software package at http://wiki.chasmsoftware.org and is hosted by the CRAVAT web server at http://www.cravat.us

Serum proprotein convertase subtilisin/kexin type 9 and cell surface low-density lipoprotein receptor: evidence for a reciprocal regulation

BACKGROUND

Proprotein convertase subtilisin/kexin type 9 (PCSK9) modulates low-density lipoprotein (LDL) receptor (LDLR) degradation, thus influencing serum cholesterol levels. However, dysfunctional LDLR causes hypercholesterolemia without affecting PCSK9 clearance from the circulation.

METHODS AND RESULTS

To study the reciprocal effects of PCSK9 and LDLR and the resultant effects on serum cholesterol, we produced transgenic mice expressing human (h) PCSK9. Although hPCSK9 was expressed mainly in the kidney, LDLR degradation was more evident in the liver. Adrenal LDLR levels were not affected, likely because of the impaired PCSK9 retention in this tissue. In addition, hPCSK9 expression increased hepatic secretion of apolipoprotein B-containing lipoproteins in an LDLR-independent fashion. Expression of hPCSK9 raised serum murine PCSK9 levels by 4.3-fold in wild-type mice and not at all in LDLR(-/-) mice, in which murine PCSK9 levels were already 10-fold higher than in wild-type mice. In addition, LDLR(+/-) mice had a 2.7-fold elevation in murine PCSK9 levels and no elevation in cholesterol levels. Conversely, acute expression of human LDLR in transgenic mice caused a 70% decrease in serum murine PCSK9 levels. Turnover studies using physiological levels of hPCSK9 showed rapid clearance in wild-type mice (half-life, 5.2 minutes), faster clearance in human LDLR transgenics (2.9 minutes), and much slower clearance in LDLR(-/-) recipients (50.5 minutes). Supportive results were obtained with an in vitro system. Finally, up to 30% of serum hPCSK9 was associated with LDL regardless of LDLR expression.

CONCLUSIONS

Our results support a scenario in which LDLR represents the main route of elimination of PCSK9 and a reciprocal regulation between these 2 proteins controls serum PCSK9 levels, hepatic LDLR expression, and serum LDL levels.

Accurate timekeeping is controlled by a cycling activator in Arabidopsis

Transcriptional feedback loops are key to circadian clock function in many organisms. Current models of the Arabidopsis circadian network consist of several coupled feedback loops composed almost exclusively of transcriptional repressors. Indeed, a central regulatory mechanism is the repression of evening-phased clock genes via the binding of morning-phased Myb-like repressors to evening element (EE) promoter motifs. We now demonstrate that a related Myb-like protein, REVEILLE8 (RVE8), is a direct transcriptional activator of EE-containing clock and output genes. Loss of RVE8 and its close homologs causes a delay and reduction in levels of evening-phased clock gene transcripts and significant lengthening of clock pace. Our data suggest a substantially revised model of the circadian oscillator, with a clock-regulated activator essential both for clock progression and control of clock outputs. Further, our work suggests that the plant clock consists of a highly interconnected, complex regulatory network rather than of coupled morning and evening feedback loops. DOI:http://dx.doi.org/10.7554/eLife.00473.001.

Accurate timekeeping is controlled by a cycling activator in Arabidopsis

Transcriptional feedback loops are key to circadian clock function in many organisms. Current models of the Arabidopsis circadian network consist of several coupled feedback loops composed almost exclusively of transcriptional repressors. Indeed, a central regulatory mechanism is the repression of evening-phased clock genes via the binding of morning-phased Myb-like repressors to evening element (EE) promoter motifs. We now demonstrate that a related Myb-like protein, REVEILLE8 (RVE8), is a direct transcriptional activator of EE-containing clock and output genes. Loss of RVE8 and its close homologs causes a delay and reduction in levels of evening-phased clock gene transcripts and significant lengthening of clock pace. Our data suggest a substantially revised model of the circadian oscillator, with a clock-regulated activator essential both for clock progression and control of clock outputs. Further, our work suggests that the plant clock consists of a highly interconnected, complex regulatory network rather than of coupled morning and evening feedback loops.

DOI:http://dx.doi.org/10.7554/eLife.00473.001

We live in a world with a 24-hr cycle in which day follows night follows day with complete predictability. Life on earth has evolved to take advantage of this predictability by using circadian clocks to prepare for the coming of night (or day), and plants are no exception. Even in constant darkness, characteristics such as leaf movements show a constant cycle of around 24 hr.

Most circadian clocks rely on negative feedback loops involving various genes and proteins to keep track of time. In one of these feedback loops, certain genes—called morning-phased genes—are expressed as proteins during the day, and these proteins prevent other genes—called evening-phased genes—from producing proteins. As night approaches, however, a second feedback loop acts to stop the morning-phased genes being expressed, thus allowing the evening-phased genes to produce proteins. And as day approaches, expression of these genes is stopped and the whole cycle starts again.

Many of the genes and proteins involved in the circadian system of Arabidopsis thaliana, a small flowering plant that is widely used as a model organism, have been identified, and its circadian clock was thought to rely almost entirely on proteins called repressors that block the transcription of genes. Now, Hsu et al. have shown that the Arabidopsis clock also involves proteins that increase the expression of certain genes at specific times of the day.

Hsu et al. focused on the promoter regions of evening-phased genes: these regions are stretches of DNA that proteins called transcription factors bind to and either encourage the expression of a gene (if the protein is a transcriptional activator) or block its expression (as a transcriptional repressor). In particular, they focused on a protein called RVE8 that is most strongly expressed in the afternoon and, based on previous research, is thought to activate the transcription of genes. Using genetically modified plants in which the gene for RVE8 can be turned on and off, they found that this protein led to increases in the expression of some genes, and reductions in the expression of others.

Further analysis showed that RVE8 was able to activate the expression of evening-phased genes directly, without requiring that new proteins be made first. By contrast, morning-expressed genes were likely to be suppressed by RVE8 via an indirect mechanism that involved other proteins that had previously been activated by RVE8. The expression of RVE8 itself is regulated by other clock genes and also by an undefined post-transcriptional process. Therefore rather than consisting of a morning feedback loop coupled to an evening feedback loop, with both loops being based on repressors, the plant clock is instead better viewed as a highly connected network of activators and repressors. Further research is clearly necessary to understand this unexpected complexity in the circadian clock of Arabidopsis.

DOI:http://dx.doi.org/10.7554/eLife.00473.002

Maternal separation with early weaning: a rodent model providing novel insights into neglect associated developmental deficits

Child neglect is the most prevalent form of child maltreatment in the United States, and poses a serious public health concern. Children who survive such episodes go on to experience long-lasting psychological and behavioral problems, including higher rates of post-traumatic stress disorder symptoms, depression, alcohol and drug abuse, attention-deficit/hyperactivity disorder, and cognitive deficits. To date, most research into the causes of these life-long problems has focused on well-established targets such as stress responsive systems, including the hypothalamus-pituitary-adrenal axis. Using the maternal separation and early weaning model, we have attempted to provide comprehensive molecular profiling of a model of early-life neglect in an organism amenable to genomic manipulation: the mouse. In this article, we report new findings generated with this model using chromatin immunoprecipitation sequencing, diffuse tensor magnetic resonance imaging, and behavioral analyses. We also review the validity of the maternal separation and early weaning model, which reflects behavioral deficits observed in neglected humans including hyperactivity, anxiety, and attentional deficits. Finally, we summarize the molecular characterization of these animals, including RNA profiling and label-free proteomics, which highlight protein translation and myelination as novel pathways of interest.

The cis-acting CTTC-P1BS module is indicative for gene function of LjVTI12, a Qb-SNARE protein gene that is required for arbuscule formation in Lotus japonicus

The majority of land plants live in symbiosis with arbuscular mycorrhizal fungi from the phylum Glomeromycota. This symbiosis improves acquisition of phosphorus (P) by the host plant in exchange for carbohydrates, especially under low-P availability. The symbiosome, constituted by root cortex cells accommodating arbuscular mycorrhizal fungal hyphae, is the site at which bi-directional exchange of nutrients and metabolites takes place. Uptake of orthophosphate (Pi) in the symbiosome is facilitated by mycorrhiza-specific plant Pi transporters. Modifications of the potato Pi transporter 3 (StPT3) promoter were analysed in transgenic mycorrhizal roots, and it was found that the CTTC cis-regulatory element is necessary and sufficient for a transcriptional response to fungal colonization under low-Pi conditions. Phylogenetic footprinting also revealed binary combination of the CTTC element with the Pi starvation response-associated PHR1-binding site (P1BS) in the promoters of several mycorrhiza-specific Pi transporter genes. Scanning of the Lotus japonicus genome for gene promoters containing both cis-regulatory elements revealed a strong over-representation of genes involved in transport processes. One of these, LjVTI12, encoding a member of the SNARE family of proteins involved in membrane transport, exhibited enhanced transcript levels in Lotus roots colonized with the arbuscular mycorrhizal fungus Glomus intraradices. Down-regulation of LjVTI12 by RNA interference resulted in a mycorrhiza-specific phenotype characterized by distorted arbuscule morphology. The results highlight cooperative cis-regulation which integrates mycorrhiza and Pi starvation signaling with vesicle trafficking in symbiosome development.

Genomic evaluation of Thermoanaerobacter spp. for the construction of designer co-cultures to improve lignocellulosic biofuel production

The microbial production of ethanol from lignocellulosic biomass is a multi-component process that involves biomass hydrolysis, carbohydrate transport and utilization, and finally, the production of ethanol. Strains of the genus Thermoanaerobacter have been studied for decades due to their innate abilities to produce comparatively high ethanol yields from hemicellulose constituent sugars. However, their inability to hydrolyze cellulose, limits their usefulness in lignocellulosic biofuel production. As such, co-culturing Thermoanaerobacter spp. with cellulolytic organisms is a plausible approach to improving lignocellulose conversion efficiencies and yields of biofuels. To evaluate native lignocellulosic ethanol production capacities relative to competing fermentative end-products, comparative genomic analysis of 11 sequenced Thermoanaerobacter strains, including a de novo genome, Thermoanaerobacter thermohydrosulfuricus WC1, was conducted. Analysis was specifically focused on the genomic potential for each strain to address all aspects of ethanol production mentioned through a consolidated bioprocessing approach. Whole genome functional annotation analysis identified three distinct clades within the genus. The genomes of Clade 1 strains encode the fewest extracellular carbohydrate active enzymes and also show the least diversity in terms of lignocellulose relevant carbohydrate utilization pathways. However, these same strains reportedly are capable of directing a higher proportion of their total carbon flux towards ethanol, rather than non-biofuel end-products, than other Thermoanaerobacter strains. Strains in Clade 2 show the greatest diversity in terms of lignocellulose hydrolysis and utilization, but proportionately produce more non-ethanol end-products than Clade 1 strains. Strains in Clade 3, in which T. thermohydrosulfuricus WC1 is included, show mid-range potential for lignocellulose hydrolysis and utilization, but also exhibit extensive divergence from both Clade 1 and Clade 2 strains in terms of cellular energetics. The potential implications regarding strain selection and suitability for industrial ethanol production through a consolidated bioprocessing co-culturing approach are examined throughout the manuscript.

The host-range, genomics and proteomics of Escherichia coli O157:H7 bacteriophage rV5

BACKGROUND

Bacteriophages (phages) have been used extensively as analytical tools to type bacterial cultures and recently for control of zoonotic foodborne pathogens in foods and in animal reservoirs.

METHODS

We examined the host range, morphology, genome and proteome of the lytic E. coli O157 phage rV5, derived from phage V5, which is a member of an Escherichia coli O157:H7 phage typing set.

RESULTS

Phage rV5 is a member of the Myoviridae family possessing an icosahedral head of 91 nm between opposite apices. The extended tail measures 121 x 17 nm and has a sheath of 44 x 20 nm and a 7 nm-wide core in the contracted state. It possesses a 137,947 bp genome (43.6 mol%GC) which encodes 233 ORFs and six tRNAs. Until recently this virus appeared to be phylogenetically isolated with almost 70% of its gene products ORFans. rV5 is closely related to coliphages Delta and vB-EcoM-FY3, and more distantly related to Salmonella phages PVP-SE1 and SSE-121, Cronobacter sakazakii phage vB_CsaM_GAP31, and coliphages phAPEC8 and phi92. A complete shotgun proteomic analysis was carried out on rV5, extending what had been gleaned from the genomic analyses. Host range studies revealed that rV5 is active against several other E. coli.

Comparative genomic analyses of Mycoplasma hyopneumoniae pathogenic 168 strain and its high-passaged attenuated strain

Background

Mycoplasma hyopneumoniae is the causative agent of porcine enzootic pneumonia (EP), a mild, chronic pneumonia of swine. Despite presenting with low direct mortality, EP is responsible for major economic losses in the pig industry. To identify the virulence-associated determinants of M. hyopneumoniae, we determined the whole genome sequence of M. hyopneumoniae strain 168 and its attenuated high-passage strain 168-L and carried out comparative genomic analyses.

Results

We performed the first comprehensive analysis of M. hyopneumoniae strain 168 and its attenuated strain and made a preliminary survey of coding sequences (CDSs) that may be related to virulence. The 168-L genome has a highly similar gene content and order to that of 168, but is 4,483 bp smaller because there are 60 insertions and 43 deletions in 168-L. Besides these indels, 227 single nucleotide variations (SNVs) were identified. We further investigated the variants that affected CDSs, and compared them to reported virulence determinants. Notably, almost all of the reported virulence determinants are included in these variants affected CDSs. In addition to variations previously described in mycoplasma adhesins (P97, P102, P146, P159, P216, and LppT), cell envelope proteins (P95), cell surface antigens (P36), secreted proteins and chaperone protein (DnaK), mutations in genes related to metabolism and growth may also contribute to the attenuated virulence in 168-L. Furthermore, many mutations were located in the previously described repeat motif, which may be of primary importance for virulence.

Conclusions

We studied the virulence attenuation mechanism of M. hyopneumoniae by comparative genomic analysis of virulent strain 168 and its attenuated high-passage strain 168-L. Our findings provide a preliminary survey of CDSs that may be related to virulence. While these include reported virulence-related genes, other novel virulence determinants were also detected. This new information will form the foundation of future investigations into the pathogenesis of M. hyopneumoniae and facilitate the design of new vaccines.

Coordinate regulation of DNA methylation and H3K27me3 in mouse embryonic stem cells

Chromatin is separated into functional domains distinguished by combinatorial patterns of post-translational histone modifications and DNA methylation. Recent studies examining multiple histone modifications have found numerous chromatin states with distinct profiles of chromatin marks and functional enrichments. There are data showing coordinate regulation between DNAme and H3K27me3, which are both involved in the establishment and maintenance of epigenetic gene silencing, but the data are conflicting. Multiple studies have presented evidence to support the theory that PRC2 and DNAme cooperate to achieve silencing, or alternatively that H3K27me3 and DNAme act antagonistically. Here we examine the effect loss of either PRC2 or DNA methyltransferase activity has on the placement of the reciprocal mark in mouse ES cells. We find that DNAme is acting globally to antagonize the placement of H3K27me3, in accordance with recently published results. At least 471,011 domains in the mouse genome acquire H3K27me3 when DNAme is diminished. Of these 466,563 have been shown to be fully methylated in wildtype ES cells, indicating the effects of DNAme on H3K27me3 are direct. In a reciprocal experiment, we examine the effect loss of PRC2 has on the placement of DNAme. In contrast to the global antagonism DNAme has on the placement of H3K27me3, loss of H3K27me3 has a modest effect on DNAme, with only 4% of genes undergoing changes in DNAme, including 861 showing increases and 552 showing losses of overall DNAme. We anticipate that integrating genomic datasets where the effect of loss of a particular epigenetic mark has on the placement of other marks will help elucidate the rules governing epigenetic regulation and what role coordinate regulation of epigenetic marks plays in development and disease.

An efficient rRNA removal method for RNA sequencing in GC-rich bacteria

Background

Next generation sequencing (NGS) technologies have revolutionized gene expression studies and functional genomics analysis. However, further improvement of RNA sequencing protocols is still desirable, in order to reduce NGS costs and to increase its accuracy. In bacteria, a major problem in RNA sequencing is the abundance of ribosomal RNA (rRNA), which accounts for 95-98% of total RNA and can therefore hinder sufficient coverage of mRNA, the main focus of transcriptomic studies. Thus, efficient removal of rRNA is necessary to achieve optimal coverage, good detection sensitivity and reliable results. An additional challenge is presented by microorganisms with GC-rich genomes, in which rRNA removal is less efficient.

Results

In this work, we tested two commercial kits for rRNA removal, either alone or in combination, on Burkholderia thailandensis. This bacterium, chosen as representative of the important Burkholderia genus, which includes both pathogenic and environmental bacteria, has a rather large (6.72 Mb) and GC-rich (67.7%) genome. Each enriched mRNA sample was sequenced through paired-end Illumina GAIIx run in duplicate, yielding between 10 and 40 million reads. We show that combined treatment with both kits allows an mRNA enrichment of more than 238-fold, enabling the sequencing of almost all (more than 90%) B. thailandensis transcripts from less than 10 million reads, without introducing any bias in mRNA relative abundance, thus preserving differential expression profile.

Conclusions

The mRNA enrichment protocol presented in this work leads to an increase in detection sensitivity up to 770% compared to total RNA; such increased sensitivity allows for a corresponding reduction in the number of sequencing reads necessary for the complete analysis of whole transcriptome expression profiling. Thus we can conclude that the MICROBExpress/Ovation combined rRNA removal method could be suitable for RNA sequencing of whole transcriptomes of microorganisms with high GC content and complex genomes enabling at the same time an important scaling down of sequencing costs.

HORDE: comprehensive resource for olfactory receptor genomics

Olfactory receptors (ORs) constitute the largest gene family in the mammalian genome. The existence of these proteins underlies the nature of, and variability in, odorant perception. The Human Olfactory Receptor Data Explorer (HORDE, http://genome.weizmann.ac.il/horde/ ) is a free online resource, which presents a complete compendium of all OR genes and pseudogenes in the genome of human and four other vertebrates. HORDE includes three parts: (1) an automated pipeline, which mines OR gene and pseudogene sequences out of complete genomes, and generates gene symbols based on sequence similarity; (2) a card generator that obtains and displays annotative information on individual ORs retrieved from external databases and relevant studies; and (3) a search engine that allows user retrieval of OR information. For human ORs, HORDE specifically addresses the universe of interindividual variation, as obtained from several sources, including whole genome sequences made possible by next-generation sequencing. This encompasses single nucleotide polymorphisms (SNP) and copy number variation (CNV), including deleterious mutational events. HORDE also hosts a number of tools designed specifically to assist in the study of OR evolution and function. In this chapter, we describe the status of HORDE (build #43). We also discuss plans for future enhancements and a road map for HORDE to become a better community-based bioinformatics tool. We highlight HORDE's role as a major research tool in the study of an expanding cohort of OR repertoires.

FastUniq: a fast de novo duplicates removal tool for paired short reads

The presence of duplicates introduced by PCR amplification is a major issue in paired short reads from next-generation sequencing platforms. These duplicates might have a serious impact on research applications, such as scaffolding in whole-genome sequencing and discovering large-scale genome variations, and are usually removed. We present FastUniq as a fast de novo tool for removal of duplicates in paired short reads. FastUniq identifies duplicates by comparing sequences between read pairs and does not require complete genome sequences as prerequisites. FastUniq is capable of simultaneously handling reads with different lengths and results in highly efficient running time, which increases linearly at an average speed of 87 million reads per 10 minutes. FastUniq is freely available at http://sourceforge.net/projects/fastuniq/.

Enhancement of carotenoid production by disrupting the C22-sterol desaturase gene (CYP61) in Xanthophyllomyces dendrorhous

BACKGROUND

Xanthophyllomyces dendrorhous is a basidiomycetous yeast that synthesizes astaxanthin, which is a carotenoid with a great biotechnological impact. The ergosterol and carotenoid synthesis pathways are derived from the mevalonate pathway, and in both pathways, cytochrome P450 enzymes are involved.

RESULTS

In this study, we isolated and described the X. dendrorhous CYP61 gene, which encodes a cytochrome P450 involved in ergosterol biosynthesis. This gene is composed of nine exons and encodes a 526 amino acid polypeptide that shares significant percentages of identity and similitude with the C22-sterol desaturase, CYP61, from other fungi. Mutants derived from different parental strains were obtained by disrupting the CYP61 gene with an antibiotic selection marker. These mutants were not able to produce ergosterol and accumulated ergosta-5,8,22-trien-3-ol and ergosta-5,8-dien-3-ol. Interestingly, all of the mutants had a more intense red color phenotype than their respective parental strains. The carotenoid composition was qualitatively and quantitatively analyzed by RP-HPLC, revealing that the carotenoid content was higher in the mutant strains without major changes in their composition. The expression of the HMGR gene, which encodes an enzyme involved in the mevalonate pathway (3-hydroxy-3-methylglutaryl-CoA reductase), was analyzed by RT-qPCR showing that its transcript levels are higher in the CYP61 mutants.

CONCLUSIONS

These results suggest that in X. dendrorhous, ergosterol regulates HMGR gene expression by a negative feedback mechanism and in this way; it contributes in the regulation of the carotenoid biosynthesis.

Transcriptional profiles of mating-responsive genes from testes and male accessory glands of the Mediterranean fruit fly, Ceratitis capitata

BACKGROUND

Insect seminal fluid is a complex mixture of proteins, carbohydrates and lipids, produced in the male reproductive tract. This seminal fluid is transferred together with the spermatozoa during mating and induces post-mating changes in the female. Molecular characterization of seminal fluid proteins in the Mediterranean fruit fly, Ceratitis capitata, is limited, although studies suggest that some of these proteins are biologically active.

METHODOLOGY/PRINCIPAL FINDINGS

We report on the functional annotation of 5914 high quality expressed sequence tags (ESTs) from the testes and male accessory glands, to identify transcripts encoding putative secreted peptides that might elicit post-mating responses in females. The ESTs were assembled into 3344 contigs, of which over 33% produced no hits against the nr database, and thus may represent novel or rapidly evolving sequences. Extraction of the coding sequences resulted in a total of 3371 putative peptides. The annotated dataset is available as a hyperlinked spreadsheet. Four hundred peptides were identified with putative secretory activity, including odorant binding proteins, protease inhibitor domain-containing peptides, antigen 5 proteins, mucins, and immunity-related sequences. Quantitative RT-PCR-based analyses of a subset of putative secretory protein-encoding transcripts from accessory glands indicated changes in their abundance after one or more copulations when compared to virgin males of the same age. These changes in abundance, particularly evident after the third mating, may be related to the requirement to replenish proteins to be transferred to the female.

CONCLUSIONS/SIGNIFICANCE

We have developed the first large-scale dataset for novel studies on functions and processes associated with the reproductive biology of Ceratitis capitata. The identified genes may help study genome evolution, in light of the high adaptive potential of the medfly. In addition, studies of male recovery dynamics in terms of accessory gland gene expression profiles and correlated remating inhibition mechanisms may permit the improvement of pest management approaches.

Expression of endogenous retroviruses is negatively regulated by the pluripotency marker Rex1/Zfp42

Rex1/Zfp42 is a Yy1-related zinc-finger protein whose expression is frequently used to identify pluripotent stem cells. We show that depletion of Rex1 levels notably affected self-renewal of mouse embryonic stem (ES) cells in clonal assays, in the absence of evident differences in expression of marker genes for pluripotency or differentiation. By contrast, marked differences in expression of several endogenous retroviral elements (ERVs) were evident upon Rex1 depletion. We demonstrate association of REX1 to specific elements in chromatin-immunoprecipitation assays, most strongly to muERV-L and to a lower extent to IAP and musD elements. Rex1 regulates muERV-L expression in vivo, as we show altered levels upon transient gain-and-loss of Rex1 function in pre-implantation embryos. We also find REX1 can associate with the lysine-demethylase LSD1/KDM1A, suggesting they act in concert. Similar to REX1 binding to retrotransposable elements (REs) in ES cells, we also detected binding of the REX1 related proteins YY1 and YY2 to REs, although the binding preferences of the two proteins were slightly different. Altogether, we show that Rex1 regulates ERV expression in mouse ES cells and during pre-implantation development and suggest that Rex1 and its relatives have evolved as regulators of endogenous retroviral transcription.

A new heterogeneous family of telomerically encoded Cryptosporidium proteins

Cryptosporidiosis is predominantly caused by two closely related species of protozoan parasites the zoonotic Cryptosporidium parvum and anthroponotic Cryptosporidium hominis which diverge phenotypically in respect to host range and virulence. Using comparative genomics we identified two genes displaying overt heterogeneity between species. Although initial work suggested both were species specific, Cops-1 for C. parvum and Chos-1 for C. hominis, subsequent study identified an abridged ortholog of Cops-1 in C. hominis. Cops-1 and Chos-1 showed limited, but significant, similarity to each other and share common features: (i) telomeric location: Cops-1 is the last gene on chromosome 2, whilst Chos-1 is the first gene on chromosome 5, (ii) encode circa 50-kDa secreted proteins with isoelectric points above 10, (iii) are serine rich, and (iv) contain internal nucleotide repeats. Importantly, Cops-1 sequence contains specific SNPs with good discriminatory power useful epidemiologically. C. parvum-infected patient sera recognized a 50-kDa protein in antigen preparations of C. parvum but not C. hominis, consistent with Cops-1 being antigenic for patients. Interestingly, anti-Cops-1 monoclonal antibody (9E1) stained oocyst content and sporozoite surface of C. parvum only. This study provides a new example of protozoan telomeres as rapidly evolving contingency loci encoding putative virulence factors.

Reduction of false-negative papillary thyroid carcinomas by the routine analysis of BRAF(T1799A) mutation on fine-needle aspiration biopsy specimens: a prospective study of 814 thyroid FNAB patients

OBJECTIVES

To evaluate prospectively the usefulness of the routine determination of BRAF(T1799A) mutation on thyroid fine-needle aspiration biopsy (FNAB) to detect cytopathology false negative papillary thyroid carcinomas (PTC) and, therefore, as a tool to improve the sensitivity of the preoperative cytopathological diagnosis of PTC in thyroid nodules.

BACKGROUND

FNAB is the most reliable diagnostic test to discriminate between malignant and benign thyroid nodules, but nondiagnostic results remain a clinical management dilemma. BRAF(T1799A) mutation is the most prevalent genetic alteration in thyroid cancers and is specific for PTC, characteristics that make it the most potentially helpful genetic tool to improve the diagnostic accuracy of FNAB.

METHODS

An exhaustive recruitment of all patients subjected to thyroid FNAB in our institution during 4 years was performed. BRAF(T1799A) mutation was determined on thyroid FNAB specimens by PCR and restriction fragment length polymorphism, plus direct sequencing in positive samples.

RESULTS

BRAF(T1799A) mutation on FNAB detected 47.2% (17/36) of PTC cases. It confirmed preoperatively 45.5% (5/11) of the PTC cases in the indeterminate category and decreased the rate of cytopathology false-negatives in 33.3% (6/18), improving the combined (BRAF(T1799A) mutation + cytopathological analysis) sensitivity of the detection of PTC on FNAB in 16.7%.

CONCLUSIONS

BRAF(T1799A) mutation improves the diagnosis of PTC on FNAB, mainly because of the detection of cytopathology false-negatives, and it can be helpful in the routine analysis of thyroid nodules, especially in clinical settings with moderate sensitivity to detect PTC on FNAB.

The contribution of CLIP2 haploinsufficiency to the clinical manifestations of the Williams-Beuren syndrome

Williams-Beuren syndrome is a rare contiguous gene syndrome, characterized by intellectual disability, facial dysmorphisms, connective-tissue abnormalities, cardiac defects, structural brain abnormalities, and transient infantile hypercalcemia. Genes lying telomeric to RFC2, including CLIP2, GTF2I and GTF2IRD1, are currently thought to be the most likely major contributors to the typical Williams syndrome cognitive profile, characterized by a better-than-expected auditory rote-memory ability, a relative sparing of language capabilities, and a severe visual-spatial constructive impairment. Atypical deletions in the region have helped to establish genotype-phenotype correlations. So far, however, hardly any deletions affecting only a single gene in the disease region have been described. We present here two healthy siblings with a pure, hemizygous deletion of CLIP2. A putative role in the cognitive and behavioral abnormalities seen in Williams-Beuren patients has been suggested for this gene on the basis of observations in a knock-out mouse model. The presented siblings did not show any of the clinical features associated with the syndrome. Cognitive testing showed an average IQ for both and no indication of the Williams syndrome cognitive profile. This shows that CLIP2 haploinsufficiency by itself does not lead to the physical or cognitive characteristics of the Williams-Beuren syndrome, nor does it lead to the Williams syndrome cognitive profile. Although contribution of CLIP2 to the phenotype cannot be excluded when it is deleted in combination with other genes, our results support the hypothesis that GTF2IRD1 and GTF2I are the main genes causing the cognitive defects associated with Williams-Beuren syndrome.

SRY protein function in sex determination: thinking outside the box

Even though the mammalian sex-determining gene Sry has been intensively studied for the two decades since its discovery, the regions outside the conserved HMG box DNA-binding domain have received less attention due to a lack of sequence conservation and of obvious structural/functional motifs. Here, we summarize the available evidence for function beyond the HMG box, identify the known and postulated biochemical functions of the non-HMG-box domains in sex determination, and present possible explanations for the puzzling diversity of these non-HMG-box domains.

Extensive variation in surface lipoprotein gene content and genomic changes associated with virulence during evolution of a novel North American house finch epizootic strain of Mycoplasma gallisepticum

Mycoplasma gallisepticum, a significant respiratory and reproductive pathogen of domestic poultry, has since 1994 been recognized as an emergent pathogen of the American house finch (Carpodacus mexicanus). Epizootic spread and pathognomonic characteristics of house finch-associated Mycoplasma gallisepticum (HFMG) have been studied as a model of an emergent to endemic pathogen in a novel host. Here we present comparative analysis of eight HFMG genomes, including one from an index isolate and seven isolates separated spatially and temporally (1994-2008) across the epizootic, and notably having differences in virulence. HFMG represented a monophyletic clade relative to sequenced poultry isolates, with genomic changes indicating a novel M. gallisepticum lineage and including unique deletions of coding sequence. Though most of the HFMG genome was highly conserved among isolates, genetic distances correlated with temporal-spatial distance from the index. The most dramatic genomic differences among HFMG involved phase-variable and immunodominant VlhA lipoprotein genes, including those variable in presence and genomic location. Other genomic differences included tandem copy number variation of a 5 kbp repeat, changes in and adjacent to the clustered regularly interspaced short palindromic repeats, and small-scale changes affecting coding potential and association of genes with virulence. Divergence of monophyletic isolates from similar time/space in the epizootic indicated local diversification of distinct HFMG sublineages. Overall, these data identify candidate virulence genes and reveal the importance of phase-variable lipoproteins during the evolution of M. gallisepticum during its emergence and dissemination in a novel host in nature, likely mediating an important role at the interface between pathogen virulence and host immunity.

A widespread occurrence of extra open reading frames in plant Ty3/gypsy retrotransposons

Long terminal repeat (LTR) retrotransposons make up substantial parts of most higher plant genomes where they accumulate due to their replicative mode of transposition. Although the transposition is facilitated by proteins encoded within the gag-pol region which is common to all autonomous elements, some LTR retrotransposons were found to potentially carry an additional protein coding capacity represented by extra open reading frames located upstream or downstream of gag-pol. In this study, we performed a comprehensive in silico survey and comparative analysis of these extra open reading frames (ORFs) in the group of Ty3/gypsy LTR retrotransposons as the first step towards our understanding of their origin and function. We found that extra ORFs occur in all three major lineages of plant Ty3/gypsy elements, being the most frequent in the Tat lineage where most (77 %) of identified elements contained extra ORFs. This lineage was also characterized by the highest diversity of extra ORF arrangement (position and orientation) within the elements. On the other hand, all of these ORFs could be classified into only two broad groups based on their mutual similarities or the presence of short conserved motifs in their inferred protein sequences. In the Athila lineage, the extra ORFs were confined to the element 3' regions but they displayed much higher sequence diversity compared to those found in Tat. In the lineage of Chromoviruses the extra ORFs were relatively rare, occurring only in 5' regions of a group of elements present in a single plant family (Poaceae). In all three lineages, most extra ORFs lacked sequence similarities to characterized gene sequences or functional protein domains, except for two Athila-like elements with similarities to LOGL4 gene and part of the Chromoviruses extra ORFs that displayed partial similarity to histone H3 gene. Thus, in these cases the extra ORFs most likely originated by transduction or recombination of cellular gene sequences. In addition, the protein domain which is otherwise associated with DNA transposons have been detected in part of the Tat-like extra ORFs, pointing to their origin from an insertion event of a mobile element.

Discrimination of juvenile yellowfin (Thunnus albacares) and bigeye (T. obesus) Tunas using mitochondrial DNA control region and liver morphology

Yellowfin tuna, Thunnus albacares (Bonnaterre, 1788) and bigeye tuna, Thunnus obesus (Lowe, 1839) are two of the most economically important tuna species in the world. However, identification of their juveniles, especially at sizes less than 40 cm, is very difficult, often leading to misidentification and miscalculation of their catch estimates. Here, we applied the mitochondrial DNA control region D-loop, a recently validated genetic marker used for identifying tuna species (Genus Thunnus), to discriminate juvenile tunas caught by purse seine and ringnet sets around fish aggregating devices (FADs) off the Southern Iloilo Peninsula in Central Philippines. We checked individual identifications using the Neighbor-Joining Method and compared results with morphometric analyses and the liver phenotype. We tested 48 specimens ranging from 13 to 31 cm fork length. Morpho-meristic analyses suggested that 12 specimens (25%) were bigeye tuna and 36 specimens (75%) were yellowfin tuna. In contrast, the genetic and liver analyses both showed that 5 specimens (10%) were bigeye tuna and 43 (90%) yellowfin tuna. This suggests that misidentification can occur even with highly stringent morpho-meristic characters and that the mtDNA control region and liver phenotype are excellent markers to discriminate juveniles of yellowfin and bigeye tunas.

An approach for searching insertions in bacterial genes leading to the phase shift of triplet periodicity

The concept of the phase shift of triplet periodicity (TP) was used for searching potential DNA insertions in genes from 17 bacterial genomes. A mathematical algorithm for detection of these insertions has been developed. This approach can detect potential insertions and deletions with lengths that are not multiples of three bases, especially insertions of relatively large DNA fragments (>100 bases). New similarity measure between triplet matrixes was employed to improve the sensitivity for detecting the TP phase shift. Sequences of 17,220 bacterial genes with each consisting of more than 1,200 bases were analyzed, and the presence of a TP phase shift has been shown in ~16% of analysed genes (2,809 genes), which is about 4 times more than that detected in our previous work. We propose that shifts of the TP phase may indicate the shifts of reading frame in genes after insertions of the DNA fragments with lengths that are not multiples of three bases. A relationship between the phase shifts of TP and the frame shifts in genes is discussed.

Identification of a Mycoplasma ovis-like organism in a herd of farmed white-tailed deer (Odocoileus virginianus) in rural Indiana

Mycoplasma ovis is a hemoplasma parasite of sheep, goats, and reindeer; however, natural hemoplasma infection in white-tailed deer has not previously been reported. Subsequent to finding many coccoid, bacillary, and ring-shaped organisms, consistent with hemotropic mycoplasmas, on RBCs from a 72-day-old female white-tailed fawn, we sought to (1) identify the putative hemoplasma observed in blood from the fawn, (2) evaluate others in the herd for hemoplasma infection, and (3) identify clinicopathologic characteristics of hemoplasma-infected white-tailed deer. EDTA-anticoagulated whole blood was collected from the fawn and 8 apparently healthy does in the same herd. CBCs were performed on 7 nonclotted samples from the fawn and 6 does. DNA was extracted from all samples, followed by PCR amplification of bacterial (16S rDNA) and protozoal (18S rDNA) genes. The nearly complete 16S rDNA product from the fawn's sample was directly sequenced and compared with known sequences in the GenBank database. Samples from the fawn and 7 of 8 does were PCR-positive using hemoplasma-specific and M ovis-specific protocols. The fawn was PCR-negative for Anaplasma spp., Babesia spp., and Theileria spp. The 16S rDNA sequence from the fawn (GenBank accession number, FJ824847) was most closely related to M ovis (AF338268), having 98.5% sequence identity. The fawn had a mild nonregenerative anemia, a neutrophilic left-shift with toxic change, aspiration bronchopneumonia, and gastrointestinal disease. Hematologic values, including blood film evaluation, in infected does were unremarkable. The M ovis-like organism may have acted as either an opportunistic or primary pathogen in the fawn. The high occurrence of subclinical infections in the does suggests that white-tailed deer may act as wildlife reservoirs for M ovis.

Next generation sequencing-based analysis of repetitive DNA in the model dioecious [corrected] plant Silene latifolia

Background

Silene latifolia is a dioceous plant with well distinguished X and Y chromosomes that is used as a model to study sex determination and sex chromosome evolution in plants. However, efficient utilization of this species has been hampered by the lack of large-scale sequencing resources and detailed analysis of its genome composition, especially with respect to repetitive DNA, which makes up the majority of the genome.

Methodology/Principal Findings

We performed low-pass 454 sequencing followed by similarity-based clustering of 454 reads in order to identify and characterize sequences of all major groups of S. latifolia repeats. Illumina sequencing data from male and female genomes were also generated and employed to quantify the genomic proportions of individual repeat families. The majority of identified repeats belonged to LTR-retrotransposons, constituting about 50% of genomic DNA, with Ty3/gypsy elements being more frequent than Ty1/copia. While there were differences between the male and female genome in the abundance of several repeat families, their overall repeat composition was highly similar. Specific localization patterns on sex chromosomes were found for several satellite repeats using in situ hybridization with probes based on k-mer frequency analysis of Illumina sequencing data.

Conclusions/Significance

This study provides comprehensive information about the sequence composition and abundance of repeats representing over 60% of the S. latifolia genome. The results revealed generally low divergence in repeat composition between the sex chromosomes, which is consistent with their relatively recent origin. In addition, the study generated various data resources that are available for future exploration of the S. latifolia genome.

Retrogenes in rice (Oryza sativa L. ssp. japonica) exhibit correlated expression with their source genes

Gene duplication occurs by either DNA- or RNA-based processes; the latter duplicates single genes via retroposition of messenger RNA. The expression of a retroposed gene copy (retrocopy) is expected to be uncorrelated with its source gene because upstream promoter regions are usually not part of the retroposition process. In contrast, DNA-based duplication often encompasses both the coding and the intergenic (promoter) regions; hence, expression is often correlated, at least initially, between DNA-based duplicates. In this study, we identified 150 retrocopies in rice (Oryza sativa L. ssp japonica), most of which represent ancient retroposition events. We measured their expression from high-throughput RNA sequencing (RNAseq) data generated from seven tissues. At least 66% of the retrocopies were expressed but at lower levels than their source genes. However, the tissue specificity of retrogenes was similar to their source genes, and expression between retrocopies and source genes was correlated across tissues. The level of correlation was similar between RNA- and DNA-based duplicates, and they decreased over time at statistically indistinguishable rates. We extended these observations to previously identified retrocopies in Arabidopsis thaliana, suggesting they may be general features of the process of retention of plant retrogenes.

Functional characterization of Fdx1: evidence for an evolutionary relationship between P450-type and ISC-type ferredoxins

Ferredoxins are ubiquitous proteins with electron transfer activity involved in a variety of biological processes. In this work, we investigated the characteristics and function of Fdx1 from Sorangium cellulosum So ce56 by using a combination of bioinformatics and of biochemical/biophysical approaches. We were able to experimentally confirm a role of Fdx1 in the iron-sulfur cluster biosynthesis by in vitro reduction studies with cluster-loaded So ce56 IscU and by transfer studies of the cluster from the latter protein to apo-aconitase A. Moreover, we found that Fdx1 can replace mammalian adrenodoxin in supporting the activity of bovine CYP11A1. This makes S. cellulosum Fdx1 the first prokaryotic ferredoxin reported to functionally interact with this mammalian enzyme. Although the interaction with CYP11A1 is non-physiological, this is-to the best of our knowledge-the first study to experimentally prove the activity of a postulated ISC-type ferredoxin in both the ISC assembly and a cytochrome P450 system. This proves that a single ferredoxin can be structurally able to provide electrons to both cytochromes P450 and IscU and thus support different biochemical processes. Combining this finding with phylogenetic and evolutionary trace analyses led us to propose the evolution of eukaryotic mitochondrial P450-type ferredoxins and ISC-type ferredoxins from a common prokaryotic ISC-type ancestor.

Remarkable expression in the colon adenocarcinoma of Hmat-Xa, a human mannosyltransferase-like gene, that is homologous to drosophila gene GC15914

We cloned a novel human mannosyltransferase-like gene, designated Hmat-Xa, as a gene homologous to the Drosophila GC15914 gene encoding the 9QVXN0 protein: see "Project Report for FY2002 on the 'Construction of Libraries of Human Genes Participating in Glycosylation' project" 43-45 (2003), New Energy and Industrial Technology Development Organization (NEDO), NEDO and Research Association for Biotechnology, Tokyo, Japan (in Japanese). After that, the GTDC1 gene, as reported by Zhao et al., DNA Cell Biol., 23, 183-187 (2004), was found to be the same as the Hmat-Xa gene. Domain EXFGI/L/VX(2)L/VE in the Hmat-Xa protein, also present in both human mannosyltransferase II/III and mannosyltransferase IV/V, which are involved in the synthesis of lipid-linked oligosaccharides, and some bacterial mannosyltransferases. A real-time PCR study of Hmat-Xa mRNA expression in human normal and tumor multiple tissue cDNA identified its tissue-specific expression and its remarkable expression in colon adenocarcinoma as compared to the normal counterpart. Thus the elevated expression of Hmat-Xa might serve as a candidate marker for colon adenocarcinoma.

Massive changes in genome architecture accompany the transition to self-fertility in the filamentous fungus Neurospora tetrasperma

A large region of suppressed recombination surrounds the sex-determining locus of the self-fertile fungus Neurospora tetrasperma. This region encompasses nearly one-fifth of the N. tetrasperma genome and suppression of recombination is necessary for self-fertility. The similarity of the N. tetrasperma mating chromosome to plant and animal sex chromosomes and its recent origin (<5 MYA), combined with a long history of genetic and cytological research, make this fungus an ideal model for studying the evolutionary consequences of suppressed recombination. Here we compare genome sequences from two N. tetrasperma strains of opposite mating type to determine whether structural rearrangements are associated with the nonrecombining region and to examine the effect of suppressed recombination for the evolution of the genes within it. We find a series of three inversions encompassing the majority of the region of suppressed recombination and provide evidence for two different types of rearrangement mechanisms: the recently proposed mechanism of inversion via staggered single-strand breaks as well as ectopic recombination between transposable elements. In addition, we show that the N. tetrasperma mat a mating-type region appears to be accumulating deleterious substitutions at a faster rate than the other mating type (mat A) and thus may be in the early stages of degeneration.

Evolutionary patterns of RNA-based duplication in non-mammalian chordates

The role of RNA-based duplication, or retroposition, in the evolution of new gene functions in mammals, plants, and Drosophila has been widely reported. However, little is known about RNA-based duplication in non-mammalian chordates. In this study, we screened ten non-mammalian chordate genomes for retrocopies and investigated their evolutionary patterns. We identified numerous retrocopies in these species. Examination of the age distribution of these retrocopies revealed no burst of young retrocopies in ancient chordate species. Upon comparing these non-mammalian chordate species to the mammalian species, we observed that a larger fraction of the non-mammalian retrocopies was under strong evolutionary constraints than mammalian retrocopies are, as evidenced by signals of purifying selection and expression profiles. For the Western clawed frog, Medaka, and Sea squirt, many retrogenes have evolved gonad and brain expression patterns, similar to what was observed in human. Testing of retrogene movement in the Medaka genome, where the nascent sex chrosomes have been well assembled, did not reveal any significant gene movement. Taken together, our analyses demonstrate that RNA-based duplication generates many functional genes and can make a significant contribution to the evolution of non-mammalian genomes.

A quantitative TaqMan PCR assay for the detection of Mycoplasma suis

AIM

To develop a TaqMan probe-based, highly sensitive and specific quantitative PCR (qPCR) assay for the detection and quantification of Mycoplasma suis in the blood of pigs.

METHODS AND RESULTS

Primers and probes specific to Myc. suis 16S rRNA gene were designed. The qPCR assay's specificity, detection limit, intra- and inter-assay variability were evaluated and its performance was compared with a Myc. suis conventional PCR assay (cPCR). Blood of two experimentally infected pigs, 40 Indiana pigs, 40 Brazilian sows and 28 peccaries were tested. The assay detected as few as ten copies of Myc. suis plasmids and was 100-fold more sensitive than the cPCR. No cross-reactivity with nontarget pig mycoplasmas was observed. An average of 1·62 × 10(11) and 2·75 × 10(8) target copies ml(-1) of blood were detected in the acutely and chronically infected pigs, respectively. Three (7·5%) pigs and 32 (80·0%) sows were positive while all peccaries were negative for Myc. suis.

CONCLUSION

The developed qPCR assay is highly sensitive and specific for Myc. suis detection and quantification.

SIGNIFICANCE AND IMPACT OF THE STUDY

TaqMan qPCR is an accurate and quick test for detection of Myc. suis infected pigs, which can be used on varied instrumentation platforms.

The chloroplast genome of the green alga Schizomeris leibleinii (Chlorophyceae) provides evidence for bidirectional DNA replication from a single origin in the chaetophorales

In the Chlorophyceae, the chloroplast genome is extraordinarily fluid in architecture and displays unique features relative to other groups of green algae. For the Chaetophorales, 1 of the 5 major lineages of the Chlorophyceae, it has been shown that the distinctive architecture of the 223,902-bp genome of Stigeoclonium helveticum is consistent with bidirectional DNA replication from a single origin. Here, we report the 182,759-bp chloroplast genome sequence of Schizomeris leibleinii, a member of the earliest diverging lineage of the Chaetophorales. Like its Stigeoclonium homolog, the Schizomeris genome lacks a large inverted repeat encoding the rRNA operon and displays a striking bias in coding regions that is associated with a bias in base composition along each strand. Our results support the notion that these two chaetophoralean genomes replicate bidirectionally from a putative origin located in the vicinity of the small subunit ribosomal RNA gene. Their shared structural characteristics were most probably inherited from the common ancestor of all chaetophoralean algae. Short dispersed repeats account for most of the 41-kb size variation between the Schizomeris and Stigeoclonium genomes, and there is no indication that homologous recombination between these repeated elements led to the observed gene rearrangements. A comparison of the extent of variation sustained by the Stigeoclonium and Schizomeris chloroplast DNAs (cpDNAs) with that observed for the cpDNAs of the chlamydomonadalean Chlamydomonas and Volvox suggests that gene rearrangements as well as changes in the abundance of intergenic and intron sequences occurred at a slower pace in the Chaetophorales than in the Chlamydomonadales.

Taking organelles apart, putting them back together and creating new ones: lessons from the endoplasmic reticulum

The endoplasmic reticulum (ER) is a highly dynamic organelle. It is composed of four subcompartments including nuclear envelope (NE), rough ER (rER), smooth ER (sER) and transitional ER (tER). The subcompartments are interconnected, can fragment and dissociate and are able to reassemble again. They coordinate with cell function by way of protein regulators in the surrounding cytosol. The activity of the many associated molecular machines of the ER as well as the fluid nature of the limiting membrane of the ER contribute extensively to the dynamics of the ER. This review examines the properties of the ER that permit its isolation and purification and the physiological conditions that permit reconstitution both in vitro and in vivo in normal and in disease conditions.

The thyroid hormone responsive protein (THRP) has a critical role in the embryogenesis of Xenopus laevis

The human and mouse homologs of the rat thyroid hormone responsive protein (THRP), c-abl-interacting protein 2 (Abi-2), are critically involved in neurological development. The Abi-2 gene is evolutionarily conserved in vertebrates, and is also found in Xenopus laevis and Drosophila melanogaster. The THRP gene is one of the few genes regulated by thyroid hormone in adult animals. Sequence analysis of the 5'-flanking region of the THRP gene identified a putative thyroid hormone response element (TRE) that is conserved between rat and human. To determine whether or not THRP regulates neural growth and development, THRP was constitutively expressed in transgenic X. laevis. Growth of most animals was halted in early neurulation while the few animals that survived the process developed into grossly malformed tadpoles. In contrast, control animals reached late embryonic stage 25. These observations suggest that THRP over-expression in early development is not compatible with completion of early embryogenesis and that a different strategy needs to be employed to investigate THRP function in this model.

Patterns of codon usage bias in Silene latifolia

Patterns of codon usage bias (CUB) convey useful information about the selection on synonymous codons induced by gene expression and contribute to an understanding of substitution patterns observed at synonymous sites. They can also be informative about the distinctive evolutionary properties of sex chromosomes such as genetic degeneration of the Y chromosome, dosage compensation, and hemizygosity of the X chromosome in males, which can affect the selection on codon usage. Here, we study CUB in Silene latifolia, a species of interest for studying the early stages of sex chromosome evolution. We have obtained a large expressed sequence tag data set containing more than 1,608 sequence fragments by 454 sequencing. Using three different methods, we conservatively define 21 preferred codons. Interestingly, the preferred codons in S. latifolia are almost identical to those in Arabidopsis thaliana, despite their long divergence time (we estimate average nonsynonymous site divergence to be 0.216, and synonymous sites are saturated). The agreement suggests that the nature of selection on codon usage has not changed significantly during the long evolutionary time separating the two species. As in many other organisms, the frequency of preferred codons is negatively correlated with protein length. For the 43 genes with both exon and intron sequences, we find a positive correlation between gene expression levels and GC content at third codon positions, but a strong negative correlation between expression and intron GC content, suggesting that the CUB we detect in S. latifolia is more likely to be due to natural selection than to mutational bias. Using polymorphism data, we detect evidence of ongoing natural selection on CUB, but we find little support for effects of biased gene conversion. An analysis of ten sex-linked genes reveals that the X chromosome has experienced significantly more unpreferred to preferred than preferred to unpreferred substitutions, suggesting that it may be evolving higher CUB. In contrast, numbers of substitutions between preferred and unpreferred codons are similar in both directions in the Y-linked genes, contrary to the expectation of genetic degeneration.

Development and application of bovine and porcine oligonucleotide arrays with protein-based annotation

The design of oligonucleotide sequences for the detection of gene expression in species with disparate volumes of genome and EST sequence information has been broadly studied. However, a congruous strategy has yet to emerge to allow the design of sensitive and specific gene expression detection probes. This study explores the use of a phylogenomic approach to align transcribed sequences to vertebrate protein sequences for the detection of gene families to design genomewide 70-mer oligonucleotide probe sequences for bovine and porcine. The bovine array contains 23,580 probes that target the transcripts of 16,341 genes, about 72% of the total number of bovine genes. The porcine array contains 19,980 probes targeting 15,204 genes, about 76% of the genes in the Ensembl annotation of the pig genome. An initial experiment using the bovine array demonstrates the specificity and sensitivity of the array.

Combining target enrichment with barcode multiplexing for high throughput SNP discovery

Background

The primary goal of genetic linkage analysis is to identify genes affecting a phenotypic trait. After localisation of the linkage region, efficient genetic dissection of the disease linked loci requires that functional variants are identified across the loci. These functional variations are difficult to detect due to extent of genetic diversity and, to date, incomplete cataloguing of the large number of variants present both within and between populations. Massively parallel sequencing platforms offer unprecedented capacity for variant discovery, however the number of samples analysed are still limited by cost per sample. Some progress has been made in reducing the cost of resequencing using either multiplexing methodologies or through the utilisation of targeted enrichment technologies which provide the ability to resequence genomic areas of interest rather that full genome sequencing.

Results

We developed a method that combines current multiplexing methodologies with a solution-based target enrichment method to further reduce the cost of resequencing where region-specific sequencing is required. Our multiplex/enrichment strategy produced high quality data with nominal reduction of sequencing depth. We undertook a genotyping study and were successful in the discovery of novel SNP alleles in all samples at uniplex, duplex and pentaplex levels.

Conclusion

Our work describes the successful combination of a targeted enrichment method and index barcode multiplexing to reduce costs, time and labour associated with processing large sample sets. Furthermore, we have shown that the sequencing depth obtained is adequate for credible SNP genotyping analysis at uniplex, duplex and pentaplex levels.

Genetack: frameshift identification in protein-coding sequences by the Viterbi algorithm

We describe a new program for ab initio frameshift detection in protein-coding nucleotide sequences. The task is to distinguish the same strand overlapping ORFs that occur in the sequence due to a presence of a frameshifted gene from the same strand overlapping ORFs that encompass true overlapping or adjacent genes. The GeneTack program uses a hidden Markov model (HMM) of genomic sequence with possibly frameshifted protein-coding regions. The Viterbi algorithm finds the maximum likelihood path that discriminates between true adjacent genes and those adjacent protein-coding regions that just appear to be separate entities due to frameshifts. Therefore, the program can identify spurious predictions made by a conventional gene-finding program misled by a frameshift. We tested GeneTack as well as two earlier developed programs FrameD and FSFind on 17 prokaryotic genomes with frameshifts introduced randomly into known genes. We observed that the average frameshift prediction accuracy of GeneTack, in terms of (Sn + Sp)/2 values, was higher by a significant margin than the accuracy of two other programs. In addition, we observed that the average accuracy of GeneTack is favorably compared with the accuracy of the FSFind-BLAST program that uses protein database search to verify predicted frameshifts, even though GeneTack does not use external evidence. GeneTack is freely available at http://topaz.gatech.edu/GeneTack/.

Heart and skeletal muscle inflammation of farmed salmon is associated with infection with a novel reovirus

Atlantic salmon (Salmo salar L.) mariculture has been associated with epidemics of infectious diseases that threaten not only local production, but also wild fish coming into close proximity to marine pens and fish escaping from them. Heart and skeletal muscle inflammation (HSMI) is a frequently fatal disease of farmed Atlantic salmon. First recognized in one farm in Norway in 1999, HSMI was subsequently implicated in outbreaks in other farms in Norway and the United Kingdom. Although pathology and disease transmission studies indicated an infectious basis, efforts to identify an agent were unsuccessful. Here we provide evidence that HSMI is associated with infection with piscine reovirus (PRV). PRV is a novel reovirus identified by unbiased high throughput DNA sequencing and a bioinformatics program focused on nucleotide frequency as well as sequence alignment and motif analyses. Formal implication of PRV in HSMI will require isolation in cell culture and fulfillment of Koch's postulates, or prevention or modification of disease through use of specific drugs or vaccines. Nonetheless, as our data indicate that a causal relationship is plausible, measures must be taken to control PRV not only because it threatens domestic salmon production but also due to the potential for transmission to wild salmon populations.

Segmental duplications in the human genome reveal details of pseudogene formation

Duplicated pseudogenes in the human genome are disabled copies of functioning parent genes. They result from block duplication events occurring throughout evolutionary history. Relatively recent duplications (with sequence similarity ≥90% and length ≥1 kb) are termed segmental duplications (SDs); here, we analyze the interrelationship of SDs and pseudogenes. We present a decision-tree approach to classify pseudogenes based on their (and their parents’) characteristics in relation to SDs. The classification identifies 140 novel pseudogenes and makes possible improved annotation for the 3172 pseudogenes located in SDs. In particular, it reveals that many pseudogenes in SDs likely did not arise directly from parent genes, but are the result of a multi-step process. In these cases, the initial duplication or retrotransposition of a parent gene gives rise to a ‘parent pseudogene’, followed by further duplication creating duplicated–duplicated or duplicated–processed pseudogenes, respectively. Moreover, we can precisely identify these parent pseudogenes by overlap with ancestral SD loci. Finally, a comparison of nucleotide substitutions per site in a pseudogene with its surrounding SD region allows us to estimate the time difference between duplication and disablement events, and this suggests that most duplicated pseudogenes in SDs were likely disabled around the time of the original duplication.

An introduction to kernel-based learning algorithms

This paper provides an introduction to support vector machines, kernel Fisher discriminant analysis, and kernel principal component analysis, as examples for successful kernel-based learning methods. We first give a short background about Vapnik-Chervonenkis theory and kernel feature spaces and then proceed to kernel based learning in supervised and unsupervised scenarios including practical and algorithmic considerations. We illustrate the usefulness of kernel algorithms by discussing applications such as optical character recognition and DNA analysis.

Origin and fate of pseudogenes in Hemiascomycetes: a comparative analysis

Background

Pseudogenes are ubiquitous genetic elements that derive from functional genes after mutational inactivation. Characterization of pseudogenes is important to understand genome dynamics and evolution, and its significance increases when several genomes of related organisms can be compared. Among yeasts, only the genome of the S. cerevisiae reference strain has been analyzed so far for pseudogenes.

Results

We present here the first comparative analysis of pseudogenes within the fully sequenced and annotated genomes of eight yeast species, spanning the entire phylogenetic range of Hemiascomycetes. A total of 871 pseudogenes were found, out of which mutational degradation patterns and consequences on the genetic repertoire of each species could be identified. We found that most pseudogenes in yeasts originate from mutational degradation of gene copies formed after species-specific duplications but duplications of pseudogenes themselves are also encountered. In all yeasts, except in Y. lipolytica, pseudogenes tend to cluster in subtelomeric regions where they can outnumber the number of functional genes from 3 to 16 times. Pseudogenes are generally not conserved between the yeast species studied (except in two cases), consistent with their large evolutionary distances, but tend to be conserved among S. cerevisiae strains. Reiterated pseudogenization of some genes is often observed in different lineages and may affect functions essential in S. cerevisiae, which are, therefore, lost in other species. Although a variety of functions are affected by pseudogenization, there is a bias towards functions involved in the adaptation of the yeasts to their environment, and towards genes of unknown functions.

Conclusions

Our work illustrates for the first time the formation of pseudogenes in different branches of hemiascomycetous yeasts, showing their limited conservation and how they testify for the adaptation of the yeasts functional repertoires.

The Legionella pneumophila LetA/LetS two-component system exhibits rheostat-like behavior

When confronted with metabolic stress, replicative Legionella pneumophila bacteria convert to resilient, infectious cells equipped for transmission. Differentiation is promoted by the LetA/LetS two-component system, which belongs to a family of signal-transducing proteins that employ a four-step phosphorelay to regulate gene expression. Histidine 307 of LetS was essential to switch on the transmission profile, but a threonine substitution at position 311 (T311M) suggested a rheostat-like function. The letS(T311M) bacteria resembled the wild type (WT) for some traits and letS null mutants for others, whereas they displayed intermediate levels of infectivity, cytotoxicity, and lysosome evasion. Although only 30 to 50% of letS(T311M) mutants became motile, flow cytometry determined that every cell eventually activated the flagellin promoter to WT levels, but expression was delayed. Likewise, letS(T311M) mutants exhibited delayed induction of RsmY and RsmZ, regulatory RNAs that relieve CsrA repression of transmission traits. Transcriptional profile analysis revealed that letS(T311M) mutants expressed the flagellar regulon and multiple other transmissive-phase loci at a higher cell density than the WT. Accordingly, we postulate that the letS(T311M) mutant may relay phosphate less efficiently than the WT LetS sensor protein, leading to sluggish gene expression and a variety of phenotypic profiles. Thus, as first described for BvgA/BvgS, rather than acting as on/off switches, this family of two-component systems exhibit rheostat activity that likely confers versatility as microbes adapt to fluctuating environments.

ADAMTSL-6 is a novel extracellular matrix protein that binds to fibrillin-1 and promotes fibrillin-1 fibril formation

ADAMTS (A disintegrin and metalloproteinase with thrombospondin motifs)-like (ADAMTSL) proteins, a subgroup of the ADAMTS superfamily, share several domains with ADAMTS proteinases, including thrombospondin type I repeats, a cysteine-rich domain, and an ADAMTS spacer, but lack a catalytic domain. We identified two new members of ADAMTSL proteins, ADAMTSL-6alpha and -6beta, that differ in their N-terminal amino acid sequences but have common C-terminal regions. When transfected into MG63 osteosarcoma cells, both isoforms were secreted and deposited into pericellular matrices, although ADAMTSL-6alpha, in contrast to -6beta, was barely detectable in the conditioned medium. Immunolabeling at the light and electron microscopic levels showed their close association with fibrillin-1-rich microfibrils in elastic connective tissues. Surface plasmon resonance analyses demonstrated that ADAMTSL-6beta binds to the N-terminal half of fibrillin-1 with a dissociation constant of approximately 80 nm. When MG63 cells were transfected or exogenously supplemented with ADAMTSL-6, fibrillin-1 matrix assembly was promoted in the early but not the late stage of the assembly process. Furthermore, ADAMTSL-6 transgenic mice exhibited excessive fibrillin-1 fibril formation in tissues where ADAMTSL-6 was overexpressed. All together, these results indicated that ADAMTSL-6 is a novel microfibril-associated protein that binds directly to fibrillin-1 and promotes fibrillin-1 matrix assembly.