Origin of minority drug-resistant HIV-1 variants in primary HIV-1 infection.

BACKGROUND: Drug-resistant human immunodeficiency virus type 1 (HIV-1) minority variants (MVs) are present in some antiretroviral therapy (ART)-naive patients. They may result from de novo mutagenesis or transmission. To date, the latter has not been proven. METHODS: MVs were quantified by allele-specific polymerase chain reaction in 204 acute or recent seroconverters from the Zurich Primary HIV Infection study and 382 ART-naive, chronically infected patients. Phylogenetic analyses identified transmission clusters. RESULTS: Three lines of evidence were observed in support of transmission of MVs. First, potential transmitters were identified for 12 of 16 acute or recent seroconverters harboring M184V MVs. These variants were also detected in plasma and/or peripheral blood mononuclear cells at the estimated time of transmission in 3 of 4 potential transmitters who experienced virological failure accompanied by the selection of the M184V mutation before transmission. Second, prevalence between MVs harboring the frequent mutation M184V and the particularly uncommon integrase mutation N155H differed highly significantly in acute or recent seroconverters (8.2% vs 0.5%; P < .001). Third, the prevalence of less-fit M184V MVs is significantly higher in acutely or recently than in chronically HIV-1-infected patients (8.2% vs 2.5%; P = .004). CONCLUSIONS: Drug-resistant HIV-1 MVs can be transmitted. To what extent the origin-transmission vs sporadic appearance-of these variants determines their impact on ART needs to be further explored.

Cryptic diversity among Western Palearctic tree frogs: postglacial range expansion, range limits, and secondary contacts of three European tree frog lineages (Hyla arborea group).

We characterize divergence times, intraspecific diversity and distributions for recently recognized lineages within the Hyla arborea species group, based on mitochondrial and nuclear sequences from 160 localities spanning its whole distribution. Lineages of H. arborea, H. orientalis, H. molleri have at least Pliocene age, supporting species level divergence. The genetically uniform Iberian H. molleri, although largely isolated by the Pyrenees, is parapatric to H. arborea, with evidence for successful hybridization in a small Aquitanian corridor (southwestern France), where the distribution also overlaps with H. meridionalis. The genetically uniform H. arborea, spread from Crete to Brittany, exhibits molecular signatures of a postglacial range expansion. It meets different mtDNA clades of H. orientalis in NE-Greece, along the Carpathians, and in Poland along the Vistula River (there including hybridization). The East-European H. orientalis is strongly structured genetically. Five geographic mitochondrial clades are recognized, with a molecular signature of postglacial range expansions for the clade that reached the most northern latitudes. Hybridization with H. savignyi is suggested in southwestern Turkey. Thus, cryptic diversity in these Pliocene Hyla lineages covers three extremes: a genetically poor, quasi-Iberian endemic (H. molleri), a more uniform species distributed from the Balkans to Western Europe (H. arborea), and a well-structured Asia Minor-Eastern European species (H. orientalis).

The transcriptional repressor REST determines the cell-specific expression of the human MAPK8IP1 gene encoding IB1 (JIP-1).

Islet-brain 1 (IB1) is the human and rat homologue of JIP-1, a scaffold protein interacting with the c-Jun amino-terminal kinase (JNK). IB1 expression is mostly restricted to the endocrine pancreas and to the central nervous system. Herein, we explored the transcriptional mechanism responsible for this preferential islet and neuronal expression of IB1. A 731-bp fragment of the 5' regulatory region of the human MAPK8IP1 gene was isolated from a human BAC library and cloned upstream of a luciferase reporter gene. This construct drove high transcriptional activity in both insulin-secreting and neuron-like cells but not in unrelated cell lines. Sequence analysis of this promoter region revealed the presence of a neuron-restrictive silencer element (NRSE) known to bind repressor zinc finger protein REST. This factor is not expressed in insulin-secreting and neuron-like cells. By mobility shift assay, we confirmed that REST binds to the NRSE present in the IB1 promoter. Once transiently transfected in beta-cell lines, the expression vector encoding REST repressed IB1 transcriptional activity. The introduction of a mutated NRSE in the 5' regulating region of the IB1 gene abolished the repression activity driven by REST in insulin-secreting beta cells and relieved the low transcriptional activity of IB1 observed in unrelated cells. Moreover, transfection in non-beta and nonneuronal cell lines of an expression vector encoding REST lacking its transcriptional repression domain relieved IB1 promoter activity. Last, the REST-mediated repression of IB1 could be abolished by trichostatin A, indicating that deacetylase activity is required to allow REST repression. Taken together, these data establish a critical role for REST in the control of the tissue-specific expression of the human IB1 gene.

Lethal skeletal dysplasia in mice and humans lacking the golgin GMAP-210.

BACKGROUND: Establishing the genetic basis of phenotypes such as skeletal dysplasia in model organisms can provide insights into biologic processes and their role in human disease. METHODS: We screened mutagenized mice and observed a neonatal lethal skeletal dysplasia with an autosomal recessive pattern of inheritance. Through genetic mapping and positional cloning, we identified the causative mutation. RESULTS: Affected mice had a nonsense mutation in the thyroid hormone receptor interactor 11 gene (Trip11), which encodes the Golgi microtubule-associated protein 210 (GMAP-210); the affected mice lacked this protein. Golgi architecture was disturbed in multiple tissues, including cartilage. Skeletal development was severely impaired, with chondrocytes showing swelling and stress in the endoplasmic reticulum, abnormal cellular differentiation, and increased cell death. Golgi-mediated glycosylation events were altered in fibroblasts and chondrocytes lacking GMAP-210, and these chondrocytes had intracellular accumulation of perlecan, an extracellular matrix protein, but not of type II collagen or aggrecan, two other extracellular matrix proteins. The similarities between the skeletal and cellular phenotypes in these mice and those in patients with achondrogenesis type 1A, a neonatal lethal form of skeletal dysplasia in humans, suggested that achondrogenesis type 1A may be caused by GMAP-210 deficiency. Sequence analysis revealed loss-of-function mutations in the 10 unrelated patients with achondrogenesis type 1A whom we studied. CONCLUSIONS: GMAP-210 is required for the efficient glycosylation and cellular transport of multiple proteins. The identification of a mutation affecting GMAP-210 in mice, and then in humans, as the cause of a lethal skeletal dysplasia underscores the value of screening for abnormal phenotypes in model organisms and identifying the causative mutations.

A NOVEL MUTATION IN BCS1L IN A PATIENT WITH AN ISOLATED MITOCHONDRIAL COMPLEX III DEFICIENCY

Background: Isolated complex III deficiencies are caused by mutations in the mitochondrial CytB gene, in the BCS1L gene coding for a CIII assembly factor and in the UQCRQ gene that codes for the ubiquinone binding protein of complex III. Objective: Description of clinical features, mitochondrial function and molecular genetic analysis in a patient with an isolated complex III deficiency. Patient: A 17 year old boy, born to consanguineous parents who presented with hypoglycemia, glycosuria, deafness, growth retardation, Fanconi Syndrome and severe lactic acidosis in the neonatal period. Methods: Activities and assembly of OXPHOS complexes were investigated spectrophotometrically and by BN-PAGE. mt-DNAwas screened for deletions. Cytochrome b (CytB) and the BCS1L gene were sequenced. Results: Isolated complex III deficiency was detected in the patient's skeletal muscle. Using BN-PAGE blotting a complex III of lower molecular weight was detected. Staining the 2D reveals a missing subunit. No mutation was detected in the mitochondrial CytB gene. Sequence analysis of BCS1L revealed a novel homozygous point mutation p.M48V. Conclusion: The patients decreased complex III activity is most likely caused by incomplete assembly of complex III due to the homozygous p. M48V mutation in the BCS1L gene.

Molecular phylogenetics of shrews (Mammalia: Soricidae) reveal timing of transcontinental colonizations.

We sequenced 2167 base pairs (bp) of mitochondrial DNA cytochrome b and 16S, and 1390 bp of nuclear genes BRCA1 and ApoB in shrews taxa (Eulipotyphla, family Soricidae). The aim was to study the relationships at higher taxonomic levels within this family, and in particular the position of difficult clades such as Anourosorex and Myosorex. The data confirmed two monophyletic subfamilies, Soricinae and Crocidurinae. In the former, the tribes Anourosoricini, Blarinini, Nectogalini, Notiosoricini, and Soricini were supported. The latter was formed by the tribes Myosoricini and Crocidurini. The genus Suncus appeared to be paraphyletic and included Sylvisorex. We further suggest a biogeographical hypothesis, which shows that North America was colonized by three independent lineages of Soricinae during middle Miocene. Our hypothesis is congruent with the first fossil records for these taxa. Using molecular dating, the first exchanges between Africa and Eurasia occurred during the middle Miocene. The last one took place in the Late Miocene, with the dispersion of the genus Crocidura through the old world.

Ongoing and future developments at the Universal Protein Resource.

The primary mission of Universal Protein Resource (UniProt) is to support biological research by maintaining a stable, comprehensive, fully classified, richly and accurately annotated protein sequence knowledgebase, with extensive cross-references and querying interfaces freely accessible to the scientific community. UniProt is produced by the UniProt Consortium which consists of groups from the European Bioinformatics Institute (EBI), the Swiss Institute of Bioinformatics (SIB) and the Protein Information Resource (PIR). UniProt is comprised of four major components, each optimized for different uses: the UniProt Archive, the UniProt Knowledgebase, the UniProt Reference Clusters and the UniProt Metagenomic and Environmental Sequence Database. UniProt is updated and distributed every 4 weeks and can be accessed online for searches or download at http://www.uniprot.org.

HLA Class I alleles associated with HCV polymorphisms and response to antiviral therapy in HCV genotype 1-infected patients

Aims: The adaptive immune response against hepatitis C virus (HCV) is significantly shaped by the host's composition of HLA alleles. Thus, the HLA phenotype is a critical determinant of viral evolution during adaptive immune pressure. Potential associations of HLA class I alleles with polymorphisms of HCV immune escape variants are largely unknown. Methods: Direct sequence analysis of the genes encoding the HCV proteins E2, NS3 and NS5B in a cohort of 159 patients with chronic HCV genotype 1 infection who were treated with pegylated interferon-alfa 2b and ribavirin in a prospective controlled trial for 48 weeks was exhibited. HLA class I genotyping was performed by strand-specific reverse hybridization with the INNO-LiPA line probe assays for HLA-A and HLA-B and by strand-specific PCR-SSP. We analyzed each amino acid position of HCV proteins using an extension of Fisher's exact test for associations with HLA alleles. In addition, associations of specific HLA alleles with inflammatory activity, liver fibrosis, HCV RNA viral load and virologic treatment outcome were investigated. Results: Separate analyses of HCV subtype 1a and 1b isolates revealed substantially different patterns of HLA-restricted polymorphisms between subtypes. Only one polymorphism within NS5B (V2758x) was significantly associated with HLA B*15 in HCV genotype 1b infected patients (adjusted p=0,048). However, a number of HLA class I-restricted polymorphisms within novel putative HCV CD8+ T cell epitopes (genotype 1a: HLA-A*11 GTRTIASPK1086-1094 [NS3], HLA-B*07 WPAPQGARSL1111-1120 [NS3]; genotype 1b: HLA-A*24 HYAPRPCGI488-496 [E2], HLA-B*44 GENETDVLL530-538 [E2], HLA-B*15 RVFTEAMTRY2757-2766 [NS5B]) were observed with high predicted epitope binding scores assessed by the web-based software SYFPEITHI (>21). Most of the identified putative epitopes were overlapping with already otherwise published epitopes, indicating a high immunogenicity of the accordant HCV protein region. In addition, certain HLA class I alleles were associated with inflammatory activity, stage of liver fibrosis, and sustained virologic response to antiviral therapy. Conclusions: HLA class I restricted HCV sequence polymorphisms are rare. HCV polymorphisms identified within putative HCV CD8+ T cell epitopes in the present study differ in their genomic distribution between genotype 1a and 1b isolates, implying divergent adaptation to the host's immune pressure on the HCV subtype level.

Genomic organization and gene expression in a chromosomal region of Leishmania major.

Little is known about the relation between the genome organization and gene expression in Leishmania. Bioinformatic analysis can be used to predict genes and find homologies with known proteins. A model was proposed, in which genes are organized into large clusters and transcribed from only one strand, in the form of large polycistronic primary transcripts. To verify the validity of this model, we studied gene expression at the transcriptional, post-transcriptional and translational levels in a unique locus of 34kb located on chr27 and represented by cosmid L979. Sequence analysis revealed 115 ORFs on either DNA strand. Using computer programs developed for Leishmania genes, only nine of these ORFs, localized on the same strand, were predicted to code for proteins, some of which show homologies with known proteins. Additionally, one pseudogene, was identified. We verified the biological relevance of these predictions. mRNAs from nine predicted genes and proteins from seven were detected. Nuclear run-on analyses confirmed that the top strand is transcribed by RNA polymerase II and suggested that there is no polymerase entry site. Low levels of transcription were detected in regions of the bottom strand and stable transcripts were identified for four ORFs on this strand not predicted to be protein-coding. In conclusion, the transcriptional organization of the Leishmania genome is complex, raising the possibility that computer predictions may not be comprehensive.

Evolutionary patterns of MHC class II B in owls and their implications for the understanding of avian MHC evolution

Owing to its special mode of evolution and central role in the adaptive immune system, the major histocompatibility complex (MHC) has become the focus of diverse disciplines such as immunology, evolutionary ecology, and molecular evolution. MHC evolution has been studied extensively in diverse vertebrate lineages over the last few decades, and it has been suggested that birds differ from the established mammalian norm. Mammalian MHC genes evolve independently, and duplication history (i.e., orthology) can usually be traced back within lineages. In birds, this has been observed in only 3 pairs of closely related species. Here we report strong evidence for the persistence of orthology of MHC genes throughout an entire avian order. Phylogenetic reconstructions of MHC class II B genes in 14 species of owls trace back orthology over tens of thousands of years in exon 3. Moreover, exon 2 sequences from several species show closer relationships than sequences within species, resembling transspecies evolution typically observed in mammals. Thus, although previous studies suggested that long-term evolutionary dynamics of the avian MHC was characterized by high rates of concerted evolution, resulting in rapid masking of orthology, our results question the generality of this conclusion. The owl MHC thus opens new perspectives for a more comprehensive understanding of avian MHC evolution.

Differential T cell receptor photoaffinity labeling among H-2Kd restricted cytotoxic T lymphocyte clones specific for a photoreactive peptide derivative. Labeling of the alpha-chain correlates with J alpha segment usage.

Using a direct binding assay based on photoaffinity labeling, we studied the interaction of T cell receptor (TCR) with a Kd-bound photoreactive peptide derivative on living cells. The Kd-restricted Plasmodium berghei circumsporozoite (PbCS) peptide 253-260 (YIPSAEKI) was reacted NH2-terminally with biotin and at the TCR contact residue Lys259 with photoreactive iodo, 4-azido salicylic acid (IASA) to make biotin-YIPSAEK(IASA)I. Cytotoxic T lymphocyte (CTL) clones derived from mice immunized with this derivative recognized this conjugate, but not a related one lacking the IASA group nor the parental PbCS peptide. The clones were Kd restricted. Recognition experiments with variant conjugates, lacking substituents from IASA, revealed a diverse fine specificity pattern and indicated that this group interacted directly with the TCR. The TCR of four clones could be photoaffinity labeled by biotin-YIPSAEK(125IASA)I. This labeling was dependent on the conjugates binding to the Kd molecule and was selective for the TCR alpha (2 clones) or beta chain (1 clone), or was common for both chains (1 clone). TCR sequence analysis showed a preferential usage of J alpha TA28 containing alpha chains that were paired with V beta 1 expressing beta chains. The TCR that were photoaffinity labeled at the alpha chain expressed these J alpha and V beta segments. The tryptophan encoded by the J alpha TA28 segment is rarely found in other J alpha segments. Moreover, we show that the IASA group interacts preferentially with tryptophan in aqueous solution. We thus propose that for these CTL clones, labeling of the alpha chain occurs via the J alpha-encoded tryptophan residue.

Pancreatic-specific expression of the glucose transporter type 2 gene: identification of cis-elements and islet-specific trans-acting factors.

A defect in glucose sensing of the pancreatic beta-cells has been observed in several animal models of type II diabetes and has been correlated with a reduced gene expression of the glucose transporter type 2 (Glut2). In a transgenic mouse model, expression of Glut2 antisense RNA in pancreatic beta-cells has recently been shown to be associated with an impaired glucose-induced insulin secretion and the development of diabetes. To identify factors that may be involved in the specific decrease of Glut2 in the beta-cells of the diabetic animal, an attempt was made to localize the cis-elements and trans-acting factors involved in the control of Glut2 expression in the endocrine pancreas. It was demonstrated by transient transfection studies that only 338 base pairs (bp) of the murine Glut2 proximal promoter are needed for reporter gene expression in pancreatic islet-derived cell lines, whereas no activity was detected in nonpancreatic cells. Three cis-elements, GTI, GTII, and GTIII, have been identified by DNAse I footprinting and gel retardation experiments within these 338 bp. GTI and GTIII bind distinct but ubiquitously expressed trans-acting factors. On the other hand, nuclear proteins specifically expressed in pancreatic cell lines interact with GTII, and their relative abundance correlates with endogenous Glut2 expression. These GTII-binding factors correspond to nuclear proteins of 180 and 90 kilodaltons as defined by Southwestern analysis. The 180-kilodalton factor is present in pancreatic beta-cell lines but not in an alpha-cell line. Mutation of the GTI or GTIII cis-elements decreases transcriptional activity directed by the 338-bp promoter, whereas mutation of GTII increases gene transcription. Thus negative and positive regulatory sequences are identified within the proximal 338 bp of the GLUT2 promoter and may participate in the islet-specific expression of the gene by binding beta-cell specific trans-acting factors.

Characterization and evolution of the novel gene family FAM90A in primates originated by multiple duplication and rearrangement events

Genomic plasticity of human chromosome 8p23.1 region is highly influenced by two groups of complex segmental duplications (SDs), termed REPD and REPP, that mediate different kinds of rearrangements. Part of the difficulty to explain the wide range of phenotypes associated with 8p23.1 rearrangements is that REPP and REPD are not yet well characterized, probably due to their polymorphic status. Here, we describe a novel primate-specific gene family, named FAM90A (family with sequence similarity 90), found within these SDs. According to the current human reference sequence assembly, the FAM90A family includes 24 members along 8p23.1 region plus a single member on chromosome 12p13.31, showing copy number variation (CNV) between individuals. These genes can be classified into subfamilies I and II, which differ in their upstream and 5′-untranslated region sequences, but both share the same open reading frame and are ubiquitously expressed. Sequence analysis and comparative fluorescence in situ hybridization studies showed that FAM90A subfamily II suffered a big expansion in the hominoid lineage, whereas subfamily I members were likely generated sometime around the divergence of orangutan and African great apes by a fusion process. In addition, the analysis of the Ka/Ks ratios provides evidence of functional constraint of some FAM90A genes in all species. The characterization of the FAM90A gene family contributes to a better understanding of the structural polymorphism of the human 8p23.1 region and constitutes a good example of how SDs, CNVs and rearrangements within themselves can promote the formation of new gene sequences with potential functional consequences.

A general definition and nomenclature for alternative splicing events

Understanding the molecular mechanisms responsible for the regulation of the transcriptome present in eukaryotic cells isone of the most challenging tasks in the postgenomic era. In this regard, alternative splicing (AS) is a key phenomenoncontributing to the production of different mature transcripts from the same primary RNA sequence. As a plethora ofdifferent transcript forms is available in databases, a first step to uncover the biology that drives AS is to identify thedifferent types of reflected splicing variation. In this work, we present a general definition of the AS event along with anotation system that involves the relative positions of the splice sites. This nomenclature univocally and dynamically assignsa specific ‘‘AS code’’ to every possible pattern of splicing variation. On the basis of this definition and the correspondingcodes, we have developed a computational tool (AStalavista) that automatically characterizes the complete landscape of ASevents in a given transcript annotation of a genome, thus providing a platform to investigate the transcriptome diversityacross genes, chromosomes, and species. Our analysis reveals that a substantial part—in human more than a quarter—ofthe observed splicing variations are ignored in common classification pipelines. We have used AStalavista to investigate andto compare the AS landscape of different reference annotation sets in human and in other metazoan species and found thatproportions of AS events change substantially depending on the annotation protocol, species-specific attributes, andcoding constraints acting on the transcripts. The AStalavista system therefore provides a general framework to conductspecific studies investigating the occurrence, impact, and regulation of AS.

In silico meets in vivo

A report of the 6th Georgia Tech-Oak Ridge National Lab International Conference on Bioinformatics 'In silico Biology: Gene Discovery and Systems Genomics', Atlanta, USA, 15-17 November, 2007.