Analysis of Hepatitis C Virus Intrahost Diversity across the Coding Region by Ultradeep Pyrosequencing

Hepatitis C virus (HCV) is the leading cause of liver disease worldwide. In this study, we analyzed four treatment-naive patients infected with subtype 1a and performed Roche/454 pyrosequencing across the coding region. We report the presence of low-level drug resistance mutations that would most likely have been missed using conventional sequencing methods. The approach described here is broadly applicable to studies of viral diversity and could help to improve the efficacy of direct-acting antiviral agents (DAA) in the treatment of HCV-infected patients.

Trypanosoma cruzi mitochondrial maxicircles display species- and strain-specific variation and a conserved element in the non-coding region

Abstract Background The mitochondrial DNA of kinetoplastid flagellates is distinctive in the eukaryotic world due to its massive size, complex form and large sequence content. Comprised of catenated maxicircles that contain rRNA and protein-coding genes and thousands of heterogeneous minicircles encoding small guide RNAs, the kinetoplast network has evolved along with an extreme form of mRNA processing in the form of uridine insertion and deletion RNA editing. Many maxicircle-encoded mRNAs cannot be translated without this post-transcriptional sequence modification. Results We present the complete sequence and annotation of the Trypanosoma cruzi maxicircles for the CL Brener and Esmeraldo strains. Gene order is syntenic with Trypanosoma brucei and Leishmania tarentolae maxicircles. The non-coding components have strain-specific repetitive regions and a variable region that is unique for each strain with the exception of a conserved sequence element that may serve as an origin of replication, but shows no sequence identity with L. tarentolae or T. brucei. Alternative assemblies of the variable region demonstrate intra-strain heterogeneity of the maxicircle population. The extent of mRNA editing required for particular genes approximates that seen in T. brucei. Extensively edited genes were more divergent among the genera than non-edited and rRNA genes. Esmeraldo contains a unique 236-bp deletion that removes the 5'-ends of ND4 and CR4 and the intergenic region. Esmeraldo shows additional insertions and deletions outside of areas edited in other species in ND5, MURF1, and MURF2, while CL Brener has a distinct insertion in MURF2. Conclusion The CL Brener and Esmeraldo maxicircles represent two of three previously defined maxicircle clades and promise utility as taxonomic markers. Restoration of the disrupted reading frames might be accomplished by strain-specific RNA editing. Elements in the non-coding region may be important for replication, transcription, and anchoring of the maxicircle within the kinetoplast network.

Unique Alterations of an Ultraconserved Non-Coding Element in the 3 ' UTR of ZIC2 in Holoprosencephaly

Coding region alterations of ZIC2 are the second most common type of mutation in holoprosencephaly (HPE). Here we use several complementary bioinformatic approaches to identify ultraconserved cis-regulatory sequences potentially driving the expression of human ZIC2. We demonstrate that an 804 bp element in the 3' untranslated region (3'UTR) is highly conserved across the evolutionary history of vertebrates from fish to humans. Furthermore, we show that while genetic variation of this element is unexpectedly common among holoprosencephaly subjects (6/528 or >1%), it is not present in control individuals. Two of six proband-unique variants are de novo, supporting their pathogenic involvement in HPE outcomes. These findings support a general recommendation that the identification and analysis of key ultraconserved elements should be incorporated into the genetic risk assessment of holoprosencephaly cases.

Analysis of plant LTR-retrotransposons at the fine-scale family level reveals individual molecular patterns

Background: Sugarcane is an important crop worldwide for sugar production and increasingly, as a renewable energy source. Modern cultivars have polyploid, large complex genomes, with highly unequal contributions from ancestral genomes. Long Terminal Repeat retrotransposons (LTR-RTs) are the single largest components of most plant genomes and can substantially impact the genome in many ways. It is therefore crucial to understand their contribution to the genome and transcriptome, however a detailed study of LTR-RTs in sugarcane has not been previously carried out. Results: Sixty complete LTR-RT elements were classified into 35 families within four Copia and three Gypsy lineages. Structurally, within lineages elements were similar, between lineages there were large size differences. FISH analysis resulted in the expected pattern of Gypsy/heterochromatin, Copia/euchromatin, but in two lineages there was localized clustering on some chromosomes. Analysis of related ESTs and RT-PCR showed transcriptional variation between tissues and families. Four distinct patterns were observed in sRNA mapping, the most unusual of which was that of Ale1, with very large numbers of 24nt sRNAs in the coding region. The results presented support the conclusion that distinct small RNA-regulated pathways in sugarcane target the lineages of LTR-RT elements. Conclusions: Individual LTR-RT sugarcane families have distinct structures, and transcriptional and regulatory signatures. Our results indicate that in sugarcane individual LTR-RT families have distinct behaviors and can potentially impact the genome in diverse ways. For instance, these transposable elements may affect nearby genes by generating a diverse set of small RNA's that trigger gene silencing mechanisms. There is also some evidence that ancestral genomes contribute significantly different element numbers from particular LTR-RT lineages to the modern sugarcane cultivar genome.

Human telomere disease due to disruption of the CCAAT box of the TERC promoter

Mutations in the coding region of telomerase complex genes can result in accelerated telomere attrition and human disease. Manifestations of telomere disease include the bone marrow failure syndromes dyskeratosis congenita and aplastic anemia, acute myeloid leukemia, liver cirrhosis, and pulmonary fibrosis. Here, we describe a mutation in the CCAAT box (GCAAT) of the TERC gene promoter in a family in which multiple members had typical features of telomeropathy. The genetic alteration in this critical regulatory sequence resulted in reduced reporter gene activity and absent binding of transcription factor NF-Y, likely responsible for reduced TERC levels, decreased telomerase activity, and short telomeres. This is the first description of a pathogenic mutation in the highly con-served CCAAT box and the first instance of a mutation in the promoter region of TERC producing a telomeropathy. We propose that current mutation-screening strategies should include gene promoter regions for the diagnosis of telomere diseases. This clinical trial was registered at www.clinicaltrials.gov as #NCT00071045. (Blood. 2012;119(13):3060-3063)

Association study between the rs165599 catechol-O-methyltransferase genetic polymorphism and schizophrenia in a Brazilian sample

Schizophrenia is a severe psychiatric disorder with frequent recurrent psychotic relapses and progressive functional impairment. It results from a poorly understood gene-environment interaction. The gene encoding catechol-O-methyltransferase (COMT) is a likely candidate for schizophrenia. Its rs165599 (A/G) polymorphism has been shown to be associated with alteration of COMT gene expression. Therefore, the present study aimed to investigate a possible association between schizophrenia and this polymorphism. The distribution of the alleles and genotypes of this polymorphism was investigated in a Brazilian sample of 245 patients and 834 controls. The genotypic frequencies were in Hardy-Weinberg equilibrium and no statistically significant differences were found between cases and controls when analyzed according to gender or schizophrenia subtypes. There was also no difference in homozygosis between cases and controls. Thus, in the sample studied, there was no evidence of any association between schizophrenia and rs165599 (A/G) polymorphism in the non-coding region 3' of the COMT gene.

Non-classical HLA-E gene variability in Brazilians: a nearly invariable locus surrounded by the most variable genes in the human genome

The non-classical human leukocyte antigen (HLA) class I genes present a very low rate of variation. So far, only 10 HLA-E alleles encoding three proteins have been described, but only two are frequently found in worldwide populations. Because of its historical background, Brazilians are very suitable for population genetic studies. Therefore, 104 bone marrow donors from Brazil were evaluated for HLA-E exons 14. Seven variation sites were found, including two known single nucleotide polymorphisms (SNPs) at positions +424 and +756 and five new SNPs at positions +170 (intron 1), +1294 (intron 3), +1625, +1645 and +1857 (exon 4). Haplotyping analysis did show eight haplotypes, three of them known as E*01:01:01, E*01:03:01 and E*01:03:02:01 and five HLA-E new alleles that carry the new variation sites. The HLA-E*01:01:01 allele was the predominant haplotype (62.50%), followed by E*01:03:02:01 (24.52%). Selective neutrality tests have disclosed an interesting pattern of selective pressures in which balancing selection is probably shaping allele frequency distributions at an SNP at exon 3 (codon 107), sequence diversity at exon 4 and the non-coding regions is facing significant purifying pressure. Even in an admixed population such as the Brazilian one, the HLA-E locus is very conserved, presenting few polymorphic SNPs in the coding region.

Molecular evolution of a malaria resistance gene (DARC) in primates

Genes involved in host-pathogen interactions are often strongly affected by positive natural selection. The Duffy antigen, coded by the Duffy antigen receptor for chemokines (DARC) gene, serves as a receptor for Plasmodium vivax in humans and for Plasmodium knowlesi in some nonhuman primates. In the majority of sub-Saharan Africans, a nucleic acid variant in GATA-1 of the gene promoter is responsible for the nonexpression of the Duffy antigen on red blood cells and consequently resistance to invasion by P. vivax. The Duffy antigen also acts as a receptor for chemokines and is expressed in red blood cells and many other tissues of the body. Because of this dual role, we sequenced a 3,000-bp region encompassing the entire DARC gene as well as part of its 5' and 3' flanking regions in a phylogenetic sample of primates and used statistical methods to evaluate the nature of selection pressures acting on the gene during its evolution. We analyzed both coding and regulatory regions of the DARC gene. The regulatory analysis showed accelerated rates of substitution at several sites near known motifs. Our tests of positive selection in the coding region using maximum likelihood by branch sites and maximum likelihood by codon sites did not yield statistically significant evidence for the action of positive selection. However, the maximum likelihood test in which the gene was subdivided into different structural regions showed that the known binding region for P. vivax/P. knowlesi is under very different selective pressures than the remainder of the gene. In fact, most of the gene appears to be under strong purifying selection, but this is not evident in the binding region. We suggest that the binding region is under the influence of two opposing selective pressures, positive selection possibly exerted by the parasite and purifying selection exerted by chemokines.

Análise de mutações nos genes TAC3 e TACR3 em pacientes com distúrbios puberais centrais idiopáticos

OBJETIVO: Investigar a presença de variantes nos genes TAC3 e TACR3, os quais codificam a NKB e seu receptor (NK3R), respectivamente, em uma coorte de pacientes com distúrbios puberais centrais idiopáticos. SUJEITOS E MÉTODOS: Duzentos e trinta e sete pacientes foram estudados: 114 com puberdade precoce central (PPC), 73 com hipogonadismo hipogonadotrófico isolado normósmico (HHI) e 50 com retardo constitucional do crescimento e desenvolvimento (RCCD). O grupo controle consistiu de 150 indivíduos brasileiros que apresentaram desenvolvimento puberal normal. O DNA genômico foi extraído de sangue periférico, e as regiões codificadoras dos genes TAC3 e TACR3 foram amplificadas e sequenciadas automaticamente. RESULTADOS: Uma variante (p.A63P) foi identificada na NKB, e quatro variantes, p.G18D, p.L58L (c.172C>T), p.W275X e p.A449S, foram identificadas no NK3R, as quais foram ausentes no grupo controle. A variante p.A63P foi identificada em uma menina com PPC, e a variante p.A449S, em uma menina com RCCD. As variantes previamente descritas, p.G18D, p.L58L e p.W275X, foram identificadas em três indivíduos com HHI normósmico do sexo masculino não relacionados. CONCLUSÃO: Variantes raras nos genes TAC3 e TACR3 foram identificadas em pacientes com distúrbios puberais centrais idiopáticos. Mutações de perda de função no gene TACR3 foram associadas com o fenótipo de HHI normósmico. Arq Bras Endocrinol Metab. 2012;56(9):646-52

Germline mutation analysis of Tpit in Poodle dogs with ACTH-dependent hyperadrenocorticism

There is a high incidence of pituitary-dependent hyperadrenocorticism (PDH) in Poodle dogs, with family members being affected by the disease, suggesting a genetic involvement. Tpit is an obligate transcription factor for the expression of pro-opiomelanocortingene and for corticotroph terminal differentiation. The aim of the present study was to screen the Tpit gene for germline mutations in Poodles with PDH. Fifty Poodle dogs (33 female, 8.71 +/- 2.8 years) with PDH and 50 healthy Poodle dogs (32 females, 9.4241 2.8 years) were studied. Genomic DNA was isolated from peripheral blood, amplified by PCR and submitted to automatic sequence. No mutation in the coding region of Tpit was found, whereas the new single nucleotide polymorphism p.S343G, in heterozygous state, was found in the same frequency in both PDH and control groups. We concluded that Tpit gain-of-function mutations are not involved in the etiology of PDH in Poodle dogs.

GH-Releasing Hormone Receptor Gene: A Novel Splice-Disrupting Mutation and Study of Founder Effects

Background: Mutations in GH-releasing hormone receptor gene (GHRHR) are emerging as the most common cause of autosomal recessive isolated GH deficiency (IGHD). Objective: To search for GHRHR mutations in patients with familial or sporadic IGHD and to investigate founder effects in recurring mutations. Methods: The coding region of GHRHR was entirely amplified and sequenced from DNA of 18 patients with IGHD (16 unrelated) with topic posterior pituitary lobe on MRI. Haplotypes containing promoter SNPs and microsatellites flanking GHRHR were analyzed in patients with c.57+1G>A (IVS1+1G>A) mutation of our previously published kindred and also a Brazilian patient and 2 previously reported Japanese sisters with c. 1146G>A (p.E382E) mutation. Results: A novel homozygous intronic GHRHR c.752-1G>A (IVS7-1G>A) mutation, predicting loss of the constitutive splice acceptor site, was identified in two siblings with IGHD. A compound heterozygous c.[57+1G>A];[1146G>A] and a heterozygous c.527C>T (p.A176V) were found in two sporadic cases. Haplotype analysis provided evidence for a founder effect for the c.57+1G>A mutation and independent recurrence for the c.1146G>A mutation. Conclusion: We report a novel splice-disrupting mutation in GHRHR in 2 siblings and provide evidence that all c.57+1G>A (IVS1+1G>A) mutant chromosomes have the same haplotype ancestor, indicating the occurrence of a founder effect in Brazilian patients with IGHD. Copyright (C) 2012 S. Karger AG, Basel

Structural studies of Helicase NS3 variants from Hepatitis C virus genotype 3 in virological sustained responder and non-responder patients

Abstract Background About 130 million people are infected with the hepatitis C virus (HCV) worldwide, but effective treatment options are not yet available. One of the most promising targets for antiviral therapy is nonstructural protein 3 (NS3). To identify possible changes in the structure of NS3 associated with virological sustained response or non-response of patients, a model was constructed for each helicase NS3 protein coding sequence. From this, the goal was to verify the interaction between helicases variants and their ligands. Findings Evidence was found that the NS3 helicase portion of non-responder patients contained substitutions in its ATP and RNA binding sites. K210E substitution can cause an imbalance in the distribution of loads, leading to a decrease in the number of ligations between the essential amino acids required for the hydrolysis of ATP. W501R substitution causes an imbalance in the distribution of loads, leading and forcing the RNA to interact with the amino acid Thr269, but not preventing binding of ribavirin inhibitor. Conclusions Useful information is provided on the genetic profiling of the HCV genotype 3, specifically the coding region of the NS3 protein, improving our understanding of the viral genome and the regions of its protein catalytic site.

Characterization and optimization of ArtinM lectin expression in Escherichia coli

Background ArtinM is a D-mannose-specific lectin from Artocarpus integrifolia seeds that induces neutrophil migration and activation, degranulation of mast cells, acceleration of wound healing, induction of interleukin-12 production by macrophages and dendritic cells, and protective T helper 1 immune response against Leishmania major, Leishmania amazonensis and Paracoccidioides brasiliensis infections. Considering the important biological properties of ArtinM and its therapeutic applicability, this study was designed to produce high-level expression of active recombinant ArtinM (rArtinM) in Escherichia coli system. Results The ArtinM coding region was inserted in pET29a(+) vector and expressed in E. coli BL21(DE3)-Codon Plus-RP. The conditions for overexpression of soluble ArtinM were optimized testing different parameters: temperatures (20, 25, 30 or 37°C) and shaking speeds (130, 200 or 220 rpm) during induction, concentrations of the induction agent IPTG (0.01-4 mM) and periods of induction (1-19 h). BL21-CodonPlus(DE3)-RP cells induced under the optimized conditions (incubation at 20°C, at a shaking speed of 130 rpm, induction with 0.4 mM IPTG for 19 h) resulted in the accumulation of large amounts of soluble rArtinM. The culture provided 22.4 mg/L of rArtinM, which activity was determined by its one-step purification through affinity chromatography on immobilized D-mannose and glycoarray analysis. Gel filtration showed that rArtinM is monomeric, contrasting with the tetrameric form of the plant native protein (jArtinM). The analysis of intact rArtinM by mass spectrometry revealed a 16,099.5 Da molecular mass, and the peptide mass fingerprint and esi-cid-ms/ms of amino acid sequences of peptides from a tryptic digest covered 41% of the total ArtinM amino acid sequence. In addition, circular dichroism and fluorescence spectroscopy of rArtinM indicated that its global fold comprises β-sheet structure. Conclusions Overall, the optimized process to express rArtinM in E. coli provided high amounts of soluble, correctly folded and active recombinant protein, compatible with large scale production of the lectin.

Dengue virus type 3 adaptive changes during epidemics in São Jose de Rio Preto, Brazil, 2006–2007

Global dengue virus spread in tropical and sub-tropical regions has become a major international public health concern. It is evident that DENV genetic diversity plays a significant role in the immunopathology of the disease and that the identification of polymorphisms associated with adaptive responses is important for vaccine development. The investigation of naturally occurring genomic variants may play an important role in the comprehension of different adaptive strategies used by these mutants to evade the human immune system. In order to elucidate this role we sequenced the complete polyprotein-coding region of thirty-three DENV-3 isolates to characterize variants circulating under high endemicity in the city of São José de Rio Preto, Brazil, during the onset of the 2006-07 epidemic. By inferring the evolutionary history on a local-scale and estimating rates of synonymous (dS) and nonsynonimous (dN) substitutions, we have documented at least two different introductions of DENV-3 into the city and detected 10 polymorphic codon sites under significant positive selection (dN/dS > 1) and 8 under significant purifying selection (dN/dS < 1). We found several polymorphic amino acid coding sites in the envelope (15), NS1 (17), NS2A (11), and NS5 (24) genes, which suggests that these genes may be experiencing relatively recent adaptive changes. Furthermore, some polymorphisms correlated with changes in the immunogenicity of several epitopes. Our study highlights the existence of significant and informative DENV variability at the spatio-temporal scale of an urban outbreak.

Stable water isotopes of precipitation and firn cores from the northern Antarctic Peninsula region as a proxy for climate reconstruction

In order to investigate the climate variability in the northern Antarctic Peninsula region, this paper focuses on the relationship between stable isotope content of precipitation and firn, and main meteorological variables (air temperature, relative humidity, sea surface temperature, and sea ice extent). Between 2008 and 2010, we collected precipitation samples and retrieved firn cores from several key sites in this region. We conclude that the deuterium excess oscillation represents a robust indicator of the meteorological variability on a seasonal to sub-seasonal scale. Low absolute deuterium excess values and the synchronous variation of both deuterium excess and air temperature imply that the evaporation of moisture occurs in the adjacent Southern Ocean. The delta O-18-air temperature relationship is complicated and significant only at a (multi)seasonal scale. Backward trajectory calculations show that air-parcels arriving at the region during precipitation events predominantly originate at the South Pacific Ocean and Bellingshausen Sea. These investigations will be used as a calibration for ongoing and future research in the area, suggesting that appropriate locations for future ice core research are located above 600 m a.s.l. We selected the Plateau Laclavere, Antarctic Peninsula as the most promising site for a deeper drilling campaign.