Gene finding in the chicken genome.

BACKGROUND: Despite the continuous production of genome sequence for a number of organisms, reliable, comprehensive, and cost effective gene prediction remains problematic. This is particularly true for genomes for which there is not a large collection of known gene sequences, such as the recently published chicken genome. We used the chicken sequence to test comparative and homology-based gene-finding methods followed by experimental validation as an effective genome annotation method. RESULTS: We performed experimental evaluation by RT-PCR of three different computational gene finders, Ensembl, SGP2 and TWINSCAN, applied to the chicken genome. A Venn diagram was computed and each component of it was evaluated. The results showed that de novo comparative methods can identify up to about 700 chicken genes with no previous evidence of expression, and can correctly extend about 40% of homology-based predictions at the 5' end. CONCLUSIONS: De novo comparative gene prediction followed by experimental verification is effective at enhancing the annotation of the newly sequenced genomes provided by standard homology-based methods.

Characterization of the rice PHO1 gene family reveals a key role for OsPHO1;2 in phosphate homeostasis and the evolution of a distinct clade in dicotyledons.

Phosphate homeostasis was studied in a monocotyledonous model plant through the characterization of the PHO1 gene family in rice (Oryza sativa). Bioinformatics and phylogenetic analysis showed that the rice genome has three PHO1 homologs, which cluster with the Arabidopsis (Arabidopsis thaliana) AtPHO1 and AtPHO1;H1, the only two genes known to be involved in root-to-shoot transfer of phosphate. In contrast to the Arabidopsis PHO1 gene family, all three rice PHO1 genes have a cis-natural antisense transcript located at the 5 ' end of the genes. Strand-specific quantitative reverse transcription-PCR analyses revealed distinct patterns of expression for sense and antisense transcripts for all three genes, both at the level of tissue expression and in response to nutrient stress. The most abundantly expressed gene was OsPHO1;2 in the roots, for both sense and antisense transcripts. However, while the OsPHO1;2 sense transcript was relatively stable under various nutrient deficiencies, the antisense transcript was highly induced by inorganic phosphate (Pi) deficiency. Characterization of Ospho1;1 and Ospho1;2 insertion mutants revealed that only Ospho1;2 mutants had defects in Pi homeostasis, namely strong reduction in Pi transfer from root to shoot, which was accompanied by low-shoot and high-root Pi. Our data identify OsPHO1;2 as playing a key role in the transfer of Pi from roots to shoots in rice, and indicate that this gene could be regulated by its cis-natural antisense transcripts. Furthermore, phylogenetic analysis of PHO1 homologs in monocotyledons and dicotyledons revealed the emergence of a distinct clade of PHO1 genes in dicotyledons, which include members having roles other than long-distance Pi transport.

Clinical pharmacogenomic testing of KRAS, BRAF and EGFR mutations by high resolution melting analysis and ultra-deep pyrosequencing

BACKGROUND: Epidermal growth factor receptor (EGFR) and its downstream factors KRAS and BRAF are mutated in several types of cancer, affecting the clinical response to EGFR inhibitors. Mutations in the EGFR kinase domain predict sensitivity to the tyrosine kinase inhibitors gefitinib and erlotinib in lung adenocarcinoma, while activating point mutations in KRAS and BRAF confer resistance to the anti-EGFR monoclonal antibody cetuximab in colorectal cancer. The development of new generation methods for systematic mutation screening of these genes will allow more appropriate therapeutic choices. METHODS: We describe a high resolution melting (HRM) assay for mutation detection in EGFR exons 19-21, KRAS codon 12/13 and BRAF V600 using formalin-fixed paraffin-embedded samples. Somatic variation of KRAS exon 2 was also analysed by massively parallel pyrosequencing of amplicons with the GS Junior 454 platform. RESULTS: We tested 120 routine diagnostic specimens from patients with colorectal or lung cancer. Mutations in KRAS, BRAF and EGFR were observed in 41.9%, 13.0% and 11.1% of the overall samples, respectively, being mutually exclusive. For KRAS, six types of substitutions were detected (17 G12D, 9 G13D, 7 G12C, 2 G12A, 2 G12V, 2 G12S), while V600E accounted for all the BRAF activating mutations. Regarding EGFR, two cases showed exon 19 deletions (delE746-A750 and delE746-T751insA) and another two substitutions in exon 21 (one showed L858R with the resistance mutation T590M in exon 20, and the other had P848L mutation). Consistent with earlier reports, our results show that KRAS and BRAF mutation frequencies in colorectal cancer were 44.3% and 13.0%, respectively, while EGFR mutations were detected in 11.1% of the lung cancer specimens. Ultra-deep amplicon pyrosequencing successfully validated the HRM results and allowed detection and quantitation of KRAS somatic mutations. CONCLUSIONS: HRM is a rapid and sensitive method for moderate-throughput cost-effective screening of oncogene mutations in clinical samples. Rather than Sanger sequence validation, next-generation sequencing technology results in more accurate quantitative results in somatic variation and can be achieved at a higher throughput scale.

The S. pombe cdc15 gene is a key element in the reorganization of F-actin at mitosis.

The S. pombe cdc15 gene is essential for cell division. cdc15ts mutants do not form a septum, but growth and nuclear division continue, leading to formation of multinucleate cells. The earliest step in septum formation and cytokinesis, rearrangement of actin to the center of the cell, is associated with appearance of hypophosphorylated cdc15p and formation of a cdc15p ring, which colocalizes with actin. Loss of cdc15p function impairs formation of the actin ring. The abundance of cdc15 mRNA varies through the cell division cycle, peaking in early mitosis before septation. Expression of cdc15 in G2-arrested cells induces actin rearrangement to the center of the cell. These data implicate cdc15p as a key element in mediating the cytoskeletal rearrangements required for cytokinesis.

Evaluation of two long synthetic merozoite surface protein 2 peptides as malaria vaccine candidates.

Merozoite surface protein 2 (MSP2) is a promising vaccine candidate against Plasmodium falciparum blood stages. A recombinant 3D7 form of MSP2 was a subunit of Combination B, a blood stage vaccine tested in the field in Papua New Guinea. A selective effect in favour of the allelic family not represented by the vaccine argued for a MSP2 vaccine consisting of both dimorphic variants. An alternative approach to recombinant manufacture of vaccines is the production of long synthetic peptides (LSP). LSP exceeding a length of well over 100 amino acids can now be routinely synthesized. Synthetic production of vaccine antigens cuts the often time-consuming steps of protein expression and purification short. This considerably reduces the time for a candidate to reach the phase of clinical trials. Here we present the evaluation of two long synthetic peptides representing both allelic families of MSP2 as potential vaccine candidates. The constructs were well recognized by human immune sera from different locations and different age groups. Furthermore, peptide-specific antibodies in human immune sera were associated with protection from clinical malaria. The synthetic fragments share major antigenic properties with native MSP2. Immunization of mice with these antigens yielded high titre antibody responses and monoclonal antibodies recognized parasite-derived MSP2. Our results justify taking these candidate poly-peptides into further vaccine development.

Characterization and binding affinities of SmLANP: a new Schistosoma mansoni member of the ANP32 family of regulatory proteins.

Members of the leucine-rich repeat protein family are involved in diverse functions including protein phosphatase 2-inhibition, cell cycle regulation, gene regulation and signalling pathways. A novel Schistosoma mansoni gene, called SmLANP, presenting homology to various genes coding for proteins that belong to the super family of leucine-rich repeat proteins, was characterized here. SmLANP was 1184bp in length as determined from cDNA and genomic sequences and encoded a 296 amino acid open reading frame that spanning from 6 to 894bp. The predicted amino acid sequence had a calculated molecular weight of 32kDa. Analysis of the predicted sequence indicated the presence of 3 leucine-rich domains (LRR) located in the N-terminal region and an aspartic acid rich region in the C-terminal end. SmLANP transcript is expressed in all stages of the S. mansoni life cycle analyzed, exhibiting the highest expression level in males. The SmLANP protein was expressed in a GST expression system and antibodies raised in mice against the recombinant protein. By immunolocalization assay, using adult worms, it was shown that the protein is mainly present in the cell nucleus through the whole body and strongly expressed along the tegument cell body nuclei of adult worms. As members of this family are usually involved in protein-protein interaction, a yeast two hybrid assay was conducted to identify putative binding partners for SmLANP. Thirty-six possible partners were identified, and a protein ATP synthase subunit alpha was confirmed by pull down assays, as a binding partner of the SmLANP protein.

Characterization of a protein of the plastid inner envelope having homology to animal inorganic phosphate, chloride and organic-anion transporters.

A protein from Arabidopsis thaliana (L.) Heynh. showing homology to animal proteins of the NaPi-1 family, involved in the transport of inorganic phosphate, chloride, glutamate and sialic acid, has been characterized. This protein, named ANTR2 (for anion transporters) was shown by chloroplast subfractionation to be localized to the plastid inner envelope in both A. thaliana and Spinacia oleracea (L.). Immunolocalization revealed that ANTR2 was expressed in the leaf mesophyll cells as well as in the developing embryo at the upturned-U stage. Five additional homologues of ANTR2 are found in the Arabidopsis genome, of which one was shown by green fluorescent protein (GFP) fusion to be also located in the chloroplast. All ANTR proteins share homology to the animal NaPi-1 family, as well as to other organic-anion transporters that are members of the Anion:Cation Symporter (ACS) family, and share the main features of transporters from this family, including the presence of 12 putative transmembrane domains and of a 7-amino acid motif in the fourth putative transmembrane domain. ANTR2 thus represent a novel protein of the plastid inner envelope that is likely to be involved in anion transport.

Molecular characterisation of Bartonella species in cats from São Luís, state of Maranhão, north-eastern Brazil

Bartonella species are fastidious bacteria that predominantly infect mammalian erythrocytes and endothelial cells and cause long-lasting bacteraemia in their reservoir hosts. Reports that describe the epidemiology of bartonellosis in Brazil are limited. This study aimed to detect and characterise Bartonella spp DNA from cat blood samples in São Luís, Maranhão, north-eastern Brazil. Among 200 cats tested for multiple genes, nine (4.5%) were positive for Bartonella spp: six cats for Bartonella henselae and three for Bartonella clarridgeiae. Based on the phylogenetic analysis of four genes, the B. henselae strain matched strains previously observed in Brazil and was positioned in the same clade as B. henselae isolates from the United States of America. Moreover, sequence alignment demonstrated that the B. clarridgeiae strain detected in the present study was the same as the one recently detected in cats from southern Brazil.

Evolutionary simulations to detect functional lineage-specific genes.

MOTIVATION: Supporting the functionality of recent duplicate gene copies is usually difficult, owing to high sequence similarity between duplicate counterparts and shallow phylogenies, which hamper both the statistical and experimental inference. RESULTS: We developed an integrated evolutionary approach to identify functional duplicate gene copies and other lineage-specific genes. By repeatedly simulating neutral evolution, our method estimates the probability that an ORF was selectively conserved and is therefore likely to represent a bona fide coding region. In parallel, our method tests whether the accumulation of non-synonymous substitutions reveals signatures of selective constraint. We show that our approach has high power to identify functional lineage-specific genes using simulated and real data. For example, a coding region of average length (approximately 1400 bp), restricted to hominoids, can be predicted to be functional in approximately 94-100% of cases. Notably, the method may support functionality for instances where classical selection tests based on the ratio of non-synonymous to synonymous substitutions fail to reveal signatures of selection. Our method is available as an automated tool, ReEVOLVER, which will also be useful to systematically detect functional lineage-specific genes of closely related species on a large scale. AVAILABILITY: ReEVOLVER is available at http://www.unil.ch/cig/page7858.html.

A histone-fold complex and FANCM form a conserved DNA-remodeling complex to maintain genome stability.

FANCM remodels branched DNA structures and plays essential roles in the cellular response to DNA replication stress. Here, we show that FANCM forms a conserved DNA-remodeling complex with a histone-fold heterodimer, MHF. We find that MHF stimulates DNA binding and replication fork remodeling by FANCM. In the cell, FANCM and MHF are rapidly recruited to forks stalled by DNA interstrand crosslinks, and both are required for cellular resistance to such lesions. In vertebrates, FANCM-MHF associates with the Fanconi anemia (FA) core complex, promotes FANCD2 monoubiquitination in response to DNA damage, and suppresses sister-chromatid exchanges. Yeast orthologs of these proteins function together to resist MMS-induced DNA damage and promote gene conversion at blocked replication forks. Thus, FANCM-MHF is an essential DNA-remodeling complex that protects replication forks from yeast to human.

Molecular characterisation of the NSP4 gene of group A human rotavirus G2P[4] strains circulating in São Paulo, Brazil, from 1994 and 2006 to 2010

Group A human rotaviruses (HuRVA) are causative agents of acute gastroenteritis. Six viral structural proteins (VPs) and six nonstructural proteins (NSPs) are produced in RV-infected cells. NSP4 is a diarrhoea-inducing viral enterotoxin and NSP4 gene analysis revealed at least 15 (E1-E15) genotypes. This study analysed the NSP4 genetic diversity of HuRVA G2P[4] strains collected in the state of São Paulo (SP) from 1994 and 2006-2010 using reverse transcription-polymerase chain reaction, sequencing and phylogenetic analysis. Forty (97.6%) G2P[4] strains displayed genotype E2; one strain (2.4%) displayed genotype E1. These results are consistent with the proposed linkage between VP4/VP7 (G2P[4]) and the NSP4 (E2) genotype of HuRVA. NSP4 phylogenetic analysis showed distinct clusters, with grouping of most strains by their genotype and collection year, and most strains from SP were clustered together with strains from other Brazilian states. A deduced amino acid sequence alignment for E2 showed many variations in the C-terminal region, including the VP4-binding domain. Considering the ability of NSP4 to generate host immunity, monitoring NSP4 variations, along with those in the VP4 or VP7 protein, is important for evaluating the circulation and pathogenesis of RV. Finally, the presence of one G2P[4]E1 strain reinforces the idea that new genotype combinations emerge through reassortment and independent segregation.

Describing ancient horizontal gene transfers at the nucleotide and gene levels by comparative pathogenicity island genometrics.

MOTIVATION: Lateral gene transfer is a major mechanism contributing to bacterial genome dynamics and pathovar emergence via pathogenicity island (PAI) spreading. However, since few of these genomic exchanges are experimentally reproducible, it is difficult to establish evolutionary scenarios for the successive PAI transmissions between bacterial genera. Methods initially developed at the gene and/or nucleotide level for genomics, i.e. comparisons of concatenated sequences, ortholog frequency, gene order or dinucleotide usage, were combined and applied here to homologous PAIs: we call this approach comparative PAI genometrics. RESULTS: YAPI, a Yersinia PAI, and related islands were compared with measure evolutionary relationships between related modules. Through use of our genometric approach designed for tracking codon usage adaptation and gene phylogeny, an ancient inter-genus PAI transfer was oriented for the first time by characterizing the genomic environment in which the ancestral island emerged and its subsequent transfers to other bacterial genera.

Identification of key amino acids of the mouse mammary tumor virus superantigen involved in the specific interaction with T-cell receptor V(beta) domains.

Mouse mammary tumor virus (MMTV) is a retrovirus encoding a superantigen that is recognized in association with major histocompatibility complex class II by the variable region of the beta chain (V(beta)) of the T-cell receptor. The C-terminal 30 to 40 amino acids of the superantigen of different MMTVs display high sequence variability that correlates with the recognition of particular T-cell receptor V(beta) chains. Interestingly, MMTV(SIM) and mtv-8 superantigens are highly homologous but have nonoverlapping T-cell receptor V(beta) specificities. To determine the importance of these few differences for specific V(beta) interaction, we studied superantigen responses in mice to chimeric and mutant MMTV(SIM) and mtv-8 superantigens expressed by recombinant vaccinia viruses. We show that only a few changes (two to six residues) within the C terminus are necessary to modify superantigen recognition by specific V(beta)s. Thus, the introduction of the MMTV(SIM) residues 314-315 into the mtv-8 superantigen greatly decreased its V(beta)12 reactivity without gain of MMTV(SIM)-specific function. The introduction of MMTV(SIM)-specific residues 289 to 295, however, induced a recognition pattern that was a mixture of MMTV(SIM)- and mtv-8-specific V(beta) reactivities: both weak MMTV(SIM)-specific V(beta)4 and full mtv-8-specific V(beta)11 recognition were observed while V(beta)12 interaction was lost. The combination of the two MMTV(SIM)-specific regions in the mtv-8 superantigen established normal MMTV(SIM)-specific V(beta)4 reactivity and completely abolished mtv-8-specific V(beta)5, -11, and -12 interactions. These new functional superantigens with mixed V(beta) recognition patterns allowed us to precisely delineate sites relevant for molecular interactions between the SIM or mtv-8 superantigen and the T-cell receptor V(beta) domain within the 30 C-terminal residues of the viral superantigen.

The polymorphic nature of HIV type 1 env V4 affects the patterns of potential N-glycosylation sites in proviral DNA at the intrahost level.

We have previously shown that env V4 from HIV-1 plasma RNA is highly heterogeneous within a single patient, due to indel-associated polymorphism. In this study, we have analyzed the variability of V4 in proviral DNA from unfractionated PBMC and sorted T and non-T cell populations within individual patients. Our data show that the degree of sequence variability and length polymorphism in V4 from HIV provirus is even higher than we previously reported in plasma. The data also show that the sequence of V4 depends largely on the experimental approach chosen. We could observe no clear trend for compartmentalization of V4 variants in specific cell types. Of interest is the fact that some variants that had been found to be predominant in plasma were not detected in any of the cell subsets analyzed. Consistently with our observations in plasma, V3 was found to be relatively conserved at both interpatient and intrapatient level. Our data show that V4 polymorphism involving insertions and deletions in addition to point mutations results in changes in the patterns of sequons in HIV-1 proviral DNA as well as in plasma RNA. These rearrangements may result in the coexistence, within the same individual, of a swarm of different V4 regions, each characterized by a different carbohydrate surface shield. Further studies are needed to investigate the mechanism responsible for the variability observed in V4 and its role in HIV pathogenesis.

Structure-function analyses point to a polynucleotide-accommodating groove essential for APOBEC3A restriction activities.

Members of the human APOBEC3 family of editing enzymes can inhibit various mobile genetic elements. APOBEC3A (A3A) can block the retrotransposon LINE-1 and the parvovirus adeno-associated virus type 2 (AAV-2) but does not inhibit retroviruses. In contrast, APOBEC3G (A3G) can block retroviruses but has only limited effects on AAV-2 or LINE-1. What dictates this differential target specificity remains largely undefined. Here, we modeled the structure of A3A based on its homology with the C-terminal domain of A3G and further compared the sequence of human A3A to those of 11 nonhuman primate orthologues. We then used these data to perform a mutational analysis of A3A, examining its ability to restrict LINE-1, AAV-2, and foreign plasmid DNA and to edit a single-stranded DNA substrate. The results revealed an essential functional role for the predicted single-stranded DNA-docking groove located around the A3A catalytic site. Within this region, amino acid differences between A3A and A3G are predicted to affect the shape of the polynucleotide-binding groove. Correspondingly, transferring some of these A3A residues to A3G endows the latter protein with the ability to block LINE-1 and AAV-2. These results suggest that the target specificity of APOBEC3 family members is partly defined by structural features influencing their interaction with polynucleotide substrates.