Synonymous co-variation across the E1/E2 gene junction of hepatitis C virus defines virion fitness

Hepatitis C virus is a positive-sense single-stranded RNA virus. The gene junction partitioning the viral glycoproteins E1 and E2 displays concurrent sequence evolution with the 3′-end of E1 highly conserved and the 5′-end of E2 highly heterogeneous. This gene junction is also believed to contain structured RNA elements, with a growing body of evidence suggesting that such structures can act as an additional level of viral replication and transcriptional control. We have previously used ultradeep pyrosequencing to analyze an amplicon library spanning the E1/E2 gene junction from a treatment naïve patient where samples were collected over 10 years of chronic HCV infection. During this timeframe maintenance of an in-frame insertion, recombination and humoral immune targeting of discrete virus sub-populations was reported. In the current study, we present evidence of epistatic evolution across the E1/E2 gene junction and observe the development of co-varying networks of codons set against a background of a complex virome with periodic shifts in population dominance. Overtime, the number of codons actively mutating decreases for all virus groupings. We identify strong synonymous co-variation between codon sites in a group of sequences harbouring a 3 bp in-frame insertion and propose that synonymous mutation acts to stabilize the RNA structural backbone.

Molecular genetics and evolution of Puumala hantavirus

Puumala virus (PUUV) is the causative agent of nephropathia epidemica (NE), a mild form of hemorrhagic fever with renal syndrome. Finland has the highest documented incidence of NE with around 1000 cases diagnosed annually. PUUV is also found in other Scandinavian countries, Central Europe and the European part of Russia. PUUV belongs to the genus Hantavirus in the family Bunyaviridae. Hantaviruses are rodent-borne viruses each carried by a specific host that is persistently and asymptomatically infected by the virus. PUUV is carried by the bank voles (Myodes glareolus, previously known as Clethrionomys glareolus). Hantaviruses have co-evolved with their carrier rodents for millions of years and these host animals are the evolutionary scene of hantaviruses. In this study, PUUV sequences were recovered from bank voles captured in Denmark and Russian Karelia to study the evolution of PUUV in Scandinavia. Phylogenetic analysis of these strains showed a geographical clustering of genetic variants following the presumable migration pattern of bank voles during the recolonization of Scandinavia after the last ice age approximately 10 000 years ago. The currently known PUUV genome sequences were subjected to in-depth phylogenetic analyses and the results showed that genetic drift seems to be the major mechanism of PUUV evolution. In general, PUUV seems to evolve quite slowly following a molecular clock. We also found evidence for recombination in the evolution of some genetic lineages of PUUV. Viral microevolution was studied in controlled virus transmission in colonized bank voles and changes in quasispecies dynamics were recorded as the virus was transmitted from one animal to another. We witnessed PUUV evolution in vivo, as one synonymous mutation became repeatedly fixed in the viral genome during the experiment. The detailed knowledge on the PUUV diversity was used to establish new sensitive and specific detection methods for this virus. Direct viral invasion of the hypophysis was demonstrated for the first time in a lethal case of NE. PUUV detection was done by immunohistochemistry, in situ hybridization and RT-nested-PCR of the autopsy tissue samples.

Sequence evolution of the CCR5 chemokine receptor gene in primates

The chemokine receptor CCR5 can serve as a coreceptor for M-tropic HIV-1 infection and both M-tropic and T-tropic SIV infection. We sequenced the entire CCR5 gene from 10 nonhuman primates: Pongo pygmaeus, Hylobates leucogenys, Trachypithecus francoisi, Trachypithecus phayrei, Pygathrix nemaeus, Rhinopithecus roxellanae, Rhinopithecus bieti, Rhinopithecus avunculus, Macaca assamensis, and Macaca arctoides. When compared with CCR5 sequences from humans and other primates, our results demonstrate that:(1) nucleotide and amino acid sequences of CCR5 among primates are highly homologous, with variations slightly concentrated on the amino and carboxyl termini; and (2) site Asp13, which is critical for CD4-independent binding of SIV gp120 to Macaca mulatta CCR5, was also present in all other nonhuman primates tested here, suggesting that those nonhuman primate CCR5s might also bind SIV gp120 without the presence of CD4. The topologies of CCR5 gene trees constructed here conflict with the putative opinion that the snub-nosed langurs compose a monophyletic group, suggesting that the CCR5 gene may not be a good genetic marker for low-level phylogenetic analysis. The evolutionary rate of CCR5 was calculated, and our results suggest a slowdown in primates after they diverged from rodents. The synonymous mutation rate of CCR5 in primates is constant, about 1.1 x 10(-9) synonymous mutations per site per year. Comparisons of K-a and K-s suggest that the CCR5 genes have undergone negative or purifying selection. K-a/K-s ratios from cercopithecines and colobines are significantly different, implying that selective pressures have played different roles in the two lineages.

不同抗旱性青稞LEA2/LEA3蛋白基因的克隆与表达

作物的抗旱性是一个多基因控制的、极为复杂的数量性状，植物对干旱在分子水平上的差异反应通过植物组织生理和细胞生物学水平，最终表现为植物抗旱性的不同。在我国，旱地农业超过耕地面积的50%，但水资源短缺，因此培育和选育抗旱高产作物是发展节水型农业最有效的途径。 青藏高原气候恶劣、年均降雨量少，也是世界大麦初生起源中心，因而蕴藏了十分丰富的与抗逆相关的种质资源材料，从这些特殊的资源材料克隆抗旱基因，不仅对培育抗旱、优质、高产大麦新品种具有重要理论意义和经济价值，而且对整个作物抗旱基础和育种应用研究都具重大促进作用。 为了筛选青稞（裸大麦，Hordeum vulgare ssp. vulgare）抗旱性材料，本研究选用来自青藏高原不同地区的84份青稞为材料，在叶片失水率(water loss rate, WLR)检测分析的基础上，选择失水率值差异显著的12个品种，通过相对含水量（relative water content, RWC）和反复干旱法评价其抗旱性，并通过植株对干旱胁迫下的丙二醛（MDA）含量和游离脯氨酸（free-proline）含量变化，了解不同抗旱性材料的生理反应特性。选择抗旱性强弱不同的品种各两份进行LEA2蛋白基因(Dhn6基因)、LEA3蛋白基因(HVA1基因)的克隆，比较LEA蛋白结构差异与作物抗旱性之间的关系。同时，对抗旱性不同的青稞品种受到干旱时间不同的失水变化率（dynamics water loss rate, DWLR）进行了检测；对抗旱性不同的青稞对照材料进行2 h、4 h、8 h和12 h的快速干旱处理，通过SYBR Green实时荧光定量RT-PCR技术对Dhn6基因、Dhn11基因、Dhn13基因和HVA1基因在不同抗旱性材料受到不同干旱时间处理后的相对表达水平进行了检测。本研究对LEA蛋白基因在抗旱性不同的青稞材料中的干旱胁迫分子水平上的差异反应进行了研究，也对植物的抗旱机理进行了初步探讨。主要研究结果如下： 1. 青稞苗期进行离体叶片失水率测定结果表明，来自青藏高原的84份青稞材料的WLR在0.086~0.205gh-1g-1DW之间。选择WLR低于0.1gh-1g-1DW和WLR高于0.18gh-1g-1DW的品种各6份，并对苗期分别进行未干旱及干旱12小时的处理。相对含水量检测结果表明，低失水率青稞材料干旱后的具有更高的相对含水量，盆栽缺水试验也显示叶片失水率低的材料耐旱能力强于失水率高的材料。通过水合茚三酮法测定离体叶片游离脯氨酸的含量，结果表明，所有品种未干旱处理时，游离脯氨酸含量差异不大（17.10~25.74 µgg-1FW）；干旱12小时后，低失水率的品种游离脯氨酸含量明显增高（32.99~53.45µgg-1FW），高失水率品种的游离脯氨酸含量与干旱前变化不明显(P<0.05)。硫代巴比妥酸法测定离体叶片丙二醛(MDA)含量，结果显示，12份所选对照品种中，丙二醛的含量在0.97~2.74nmolg-1FW，干旱12小时后丙二醛的含量显著上升（1.46~4.74nmolg-1FW），高失水率的6个品种的丙二醛含量在未干旱和干旱处理时都明显高于低WLR品种。本研究结果表明青稞的低失水率、低丙二醛含量、高相对含水量和高脯氨酸含量具相关性（P<0.05）。综上研究，我们认为作物失水率的测定可以作为快速检测作物抗旱性的指标之一，因此，强抗旱品种喜玛拉10号(TR1)、品比14号(TR2)和弱抗旱品种冬青8号(TS1)、QB24 (TS2)被选作抗旱基因克隆和表达分析的研究材料。 2. 高等植物胚胎发育晚期丰富蛋白(late embryogenesis abundant proteins, LEA proteins)与植物耐脱水性密切相关，为了探讨青稞LEA蛋白结构差异性与植物抗旱性的关系，本研究以强抗旱品种（喜玛拉10号、品比14号）和弱抗旱品种（冬青8号、QB24）为材料，利用同源克隆法，通过RT-PCR，分别克隆了与抗旱性密切相关的Dhn6基因和HVA1基因。Dhn6基因序列分析结果表明，强抗旱品种品比14号和弱抗旱品种冬青8号Dhn6基因所克隆到的序列为1026bp，它们之间只有5个碱基的差异；喜玛拉10号和QB24克隆到的序列长963bp。在强弱不同的抗旱品种中有22个核苷酸易突变位点，相应的脱水素氨基酸序列推导结果表明，22个核苷酸突变位点中，仅有8个位点导致相应的氨基酸残基的改变，其余的位点系同义突变，另外，21个富含甘氨酸序列的缺失并没有联系作物抗旱性特征。推测这些同义突变位点的氨基酸残基对维持青稞DHN6蛋白的正常结构和功能起着非常重要的作用，也可能DHN6蛋白对青稞长期适应逆境胁迫和遗传进化的结果。对HVA1基因的序列分析结果表明，冬青8号、QB24、品比14号和喜玛拉10号的目的基因核苷酸序列全长分别为661bp、697bp、694bp和691bp，它们都包含1个完整的开放阅读框。相应的LEA3蛋白氨基酸序列结果表明，11个高度保守的氨基酸残基组成基元重复序列的拷贝数与青稞抗旱性之间没有必然关系，在强抗旱品种（喜玛拉10号、品比14号）中三个共同的氨基酸突变位点Gln32、Arg33和Ala195可能对抗旱蛋白的结构和功能有影响；另外，强抗旱青稞品种LEA3蛋白质中11-氨基酸保守基元序列拷贝数和极性氨基酸占蛋白的比例更高，推测LEA3蛋白中基元序列拷贝数和极性氨基酸占蛋白的比例对该蛋白的结构和功能影响更大。 3. LEA蛋白基因的表达水平的上调与植物的耐脱水性密切相关，我们对强抗旱性材料（喜玛拉10号、品比14号）和弱抗旱材料（冬青8号、QB24）进行干旱处理2 h、4 h、6 h、8 h和10 h的失水变化率进行测定，结果表明弱抗旱品种在2~4小时之间失水率变化最明显，而四个对照品种的失水率在8小时后和24小时的失水率值变化不大。进一步提取青稞苗期进行2 h、4 h、8 h和12 h的干旱处理后的总RNA，通过SYBR Green实时荧光定量RT-PCR技术对青稞脱水素基因(Dhn6、Dhn11和Dhn13)和LEA3蛋白基因(HVA1)的相对表达水平受干旱时间和作物抗旱性的影响进行了检测。研究发现，抗旱性不同的青稞品种随干旱处理的时间延长，Dhn6、Dhn11、Dhn13和HVA1基因的相对表达水平不同。 Dhn6基因的相对表达水平在强抗旱青稞品种干旱8小时后快速上升，但在弱抗旱青稞品种干旱处理12小时后检测到更高表达量；Dhn11基因在对照青稞抗旱品种的表达累积水平随干旱时间的延长持续下降；整个干旱过程中，Dhn13基因的相对表达水平在弱抗旱品种持续上升，在强抗旱品种中干旱处理8小时快速上升并达到最高，干旱12小时后降低。与脱水素基因相比较，强抗旱青稞品种在干旱2小时后HVA1基因的相对表达水平显著升高，相对表达量随干旱处理的时间持续上升，在干旱12小时后达到最高；与之相比较，在整个干旱过程中，弱抗旱品种的相对表达水平显著低于强抗旱品种，在干旱8小时之前弱抗旱品种的相对表达水平变化不明显；在干旱8~12小时后却显著上升。上述结果表明，不同的LEA蛋白在植物耐脱水过程中的干旱表达累积水平不同；干旱不是诱导高等植物Dhn11基因表达的主要因素；植物的抗旱性不同，不同LEA蛋白基因对干旱的反应有差异。推测某些LEA蛋白基因的干旱胁迫早期表达累积程度与植物的抗旱性直接相关；其中，Dhn11基因和Dhn12基因不同的表达模式可能与干旱调控表达顺式作用成分(dehydration responsive element, DRE)的有无或结构上的差异有关。 本研究结果认为，(1)失水率和相对含水量可作为植物抗旱性检测的指标之一；(2) DHN6同义突变位点的氨基酸残基对维持该蛋白的正常结构和功能起着重要作用；(3) 11-氨基酸保守基元序列拷贝数和极性氨基酸的比例对LEA3蛋白结构和功能有重要影响；(4)LEA蛋白表达随着干旱胁迫程度而增加，但Dhn11基因并不受干旱诱导表达；(5)作物的抗旱性不同，LEA蛋白对干旱的累积反应并不相同，干旱早期LEA蛋白的累积程度可能会影响植物的抗旱性。 Drought resistance was a complex trait which involved multiple physiological and biochemical mechanisms and regulation of numerous genes. Because its complex traits, it is difficult to understand the mechanisms of drought resistance in plants. Plants respond to water stress through multiple physiological mechanisms at the cellular, tissue, and whole-plant levels. Tibetan hulless barley, a pure line, is a selfing annual plant that has predominantly penetrated into the Qinghai-Tibetan Plateau and remains stable populations there. The wide ecological range of Tibetan hulless barley differs in water availability, temperature, soil type and vegetation, which makes it possess a high potential of adaptive diversity to abiotic stresses. This adaptive genetic diversity indicates that the potential of Tibetan hulless barley serves as a good source for drought resistance alleles for breeding purposes. 12 contrasting drought-tolerant genotypes were selected to measure relative water content (RWC), maldondialdehyde (MDA) and proline content, based on values of water loss rate (WLR) and repeated drought methods from Tibetan populations of cultivated hulless barley. As a result of the screening, sensitive and tolerant genotypes were identified to clarify relationships between characteristics of LEA2/LEA3 genes sequences and expression and drought-tolerant genotypes, associated with resistance to water deficit. In addition, dynamics water loss rate (DWLR) was measured to observe the changes on diffrential drought-tolerant genotypes. Real-time quantitative RT-PCR was applied to detect relative expression levels of Dhn6, Dhn11, Dhn13 and HVA1 genes in sensitive and tolerant genotypes with 2 h, 4 h, 8h and 12 h of dehydration. In the present study, differential sequences and expression of LEA2/LEA3 genes were explored in Tibetan hulless barley, associated with phenotypically diverse drought-tolerant genotypes. 1. The assessments of WLR and RWC were considered as an alternative measure of plant water statues reflecting the metabolic activity in plants, and the parameters of MDA and proline contents were usually consistent with the resistance to water stress. The values of detached leaf WLR of the tested genotypes were highly variable among 84 genotypes, ranging from 0.086 to 0.205 g/h.g DW. The 12 most contrasting genotypes (6 genotypes with the lowest values of WLR and 6 genotypes with the highest values of WLR) were further validated by measuring RWC, MDA and free-proline contents, which were well watered and dehydrated for 12 h. Results of RWC indicated that the values of 12 contrasting genotypes RWC ranged from 89.94% to 93.38% under condition of well water, without significant differences, but 6 genotypes with lower WLR had higher RWC suffered from 12 h dehydration. The results indicated that lower MDA contents, lower scores of WLR and higher proline contents were associated with drought-tolerant genotypes in hulless barley. Remarkably, proline amounts were increased more notable in 6 tolerant genotypes than 6 sensitive genotypes after excised leaves were dehydrated for 12 h, with control to slight changes under condition of well water. Results of MDA contents showed that six 6 tolerant genotypes had lower MDA contents than the 6 sensitive genotypes under both stressed and non-stressed conditions. As a result of that screening, drought- resistant genotypes (Ximala 10 and Pinbi 14) and drought-sensitive genotypes (Dongqing 8 and QB 24) were chosen for comparing the differential characteristics of LEA2/LEA3 genes and their expression analysis. It was conclusion that measurements of WLR could be considered an alternative index as screening of drought-tolerant genotypes in crops. 2. Late embryogenesis abundant (LEA) proteins were thought to protect against water stress in plants. To explore the relationships between configuration of LEA proteins and phenotypically diverse drought-tolerant genotypes, sequences of LEA genes and their deduced proteins were compared in Tibetan hulless barley. Results of comparing Dhn6 gene in Ximala 10 and QB24 indicated that absence of 63bp was found, except that only 5 mutant nucleotides were found. While 22 mutant sites were taken place in Dhn6 gene between sensitive and tolerant lines, 14 synonymous mutation sites appeared in the contrasting genotypes. The additional/absent polypeptide of 21 polar amino acid residues was not consistent with phenotypically drought-tolerant genotypes in hulless barley. It was deduced that synonymous mutation sites would play important roles in holding out right configurations and functions on DHN6 protein. The sequencing analysis results indicated that each cloned HVA1 gene from four selected genotypes contained an entire open reading frame. The whole sequence of HVA1 gene from Dongqing 8, QB24, Pinbi 14 and Ximala 10 was respectively 661bp, 697bp, 694bp and 691bp. Results of DNA sequence analyses showed that the differences in nucleotides of HVA1 gene in sensitive genotypes were not consistent with that of tolerant genotypes, except for absence of 33 nucleotides from +154 to +186 (numbering from ATG) in QB24. Database searches using deduced amino acid sequences showed a high homology in LEA3 proteins in the selected genotypes. Multiple sequence alignments revealed that LEA3 protein from Dongqing 8 was composed of 8 repeats of an 11 amino acid motif, less the fourth motif than Pinbi 14, Ximala 10 and QB24. Consistent mutant amino acid residues appeared in contrasting genotypes by aligning and comparing the coding sequence region, including Gln32, Arg33 and Ala195 in tolerant genotypes as compared to Asp32, Glu33 and Thr195 (Thr184 in Dongqing 8) in sensitive lines. It was concluded that consistent appearance of Gln32, Arg33 and Ala195 would contributed to functions of LEA3 protein in crops, as well as higher proportion of 11-amino-repeating motifs and polar amino acid residues. 3. Most of the LEA genes are up-regulated by dehydration, salinity, or low temperature, are also induced by application of exogenous ABA, which increases in concentration in plants under various stress conditions and acts as a mobile stress signal. Higher levels of proteins of LEA group 3 accumulated was correlated well with high level of desiccation tolerance in severely dehydrated plant seedlings. Dehydrins (DHNs), members of LEA2 protein, are an immunologically distinct protein family, and Dhn genes expression is associated with plant response to dehydration. Dynamic water loss rate was measured between sensitive genotypes and tolerant genotypes after they were dehydrated for 2 h, 4 h, 6h and 8 h. Detailed measurements of WLR at the early stage of dehydration (2, 4, 6, and 8 h) showed that WLR was stabilizing after 8 h, and there were no significant changes between these values and WLR after 24 h. Drought stress was applied to 10-day-old seedlings by draining the solution from the container for defined dehydration periods. Leaf tissues of the selected genotypes were harvested from control plants (time 0); and after 2, 4, 8, and 12 h of dehydration. Differential expression trends of Dhn6, Dhn11, Dhn13 and HVA1 genes were detected in phenotypically diverse drought-tolerant hulless barleys, related to different time of dehydration. Results of quantitative real-time PCR indicated that relative level of HVA1 expression was always higher in tolerant genotypes, rapidly increasing at the earlier stages (after 2-4 h of dehydration). However, HVA1 expressions of sensitive genotypes had a fast increase from 8 h to 12 h of stress. Significant differences in expression trends of dehydrin genes between tolerant genotypes and sensitive lines were detected, mainly in Dhn6 and Dhn13 gene, depending on the duration of the dehydration stress. The relative expression levels of Dhn6 gene were significantly higher in tolerant genotypes after 8 h dehydration, by control with notable higher expression levels after 12 h water stress in sensitive ones. The relative expression levels of Dhn13 gene tended to ascend during exposure to dehydration in drought-sensitive genotypes. However, fluctuate trends of Dhn13 expression level were detected in drought-resistant lines, including in lower expression levels of 12 h dehydration as compared to 8 h water stress. It was conclusion that (1) diverse LEA proteins would play variable roles in resisting water stress in plants; (2) expression of Dhn11 gene was not induced by dehydrated signals because of the trends of expression descended in contrasting genotypes suffered from water deficit and (3) variable accumulations on LEA proteins would be appear in diverse drought-tolerant genotypes during dehydrations. It is deduced that higher accumulations of Dhn6 and Dhn13 expression in 8 h dehydration are related to diverse drought-tolerant lines in crops. The present results indicated that different dehydrin genes would play variable functional roles in resisting water stress when plants were suffered from water deficit. The authors suggest physiologically different reactions between resistant and sensitive genotypes may be the results of differential expression of drought-resistant genes and related signal genes in plants. In addition, contrarily induced expression of Dhn11 and Dhn12 was related to dehydration responsive element (DRE) in barleys. The present study indicated that (1) measurements of WLR and RWC could be considered as one index of drought-tolerant screenings; (2) synonymous mutation sites would play important roles in holding out right configurations and functions on DHN6 protein, (3) higher proportion of 11-amino-repeating motifs and polar amino acid residues would contribute to functions on LEA3 protein, (4) the longer drought, the more accumulation on LEA proteins, except for Dhn11 gene in crops and (5) differential responses on expression of LEA protein genes would result in physiological traits of drought tolerance in plants.

Identificação de SNPs (Single Nucleiotide Polymorphisms) no gene colina monooxigenase relacionado ao metabolismo da glicina betaína em Eucalyptus

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Analysis of energetically biased transcripts of viruses and transposable elements

RNA interference (RNAi) is a natural endogenous process by which double-stranded RNA molecules trigger potent and specific gene silencing in eukaryotic cells and is characterized by target RNA cleavage. In mammals, small interfering RNAs (siRNAs) are the trigger molecules of choice and constitute a new class of RNA-based antiviral agents. In an efficient RNAi response, the antisense strand of siRNAs must enter the RNA-induced silencing complex (RISC) in a process mediated by thermodynamic features. In this report, we hypothesize that silent mutations capable of inverting thermodynamic properties can promote resistance to siRNAs. Extensive computational analyses were used to assess whether continuous selective pressure that promotes such mutations could lead to the emergence of viral strains completely resistant to RNAi (i.e., prone to transfer only the sense strands to RISC). Based on our findings, we propose that, although synonymous mutations may produce functional resistance, this strategy cannot be systematically adopted by viruses since the longest RNAi-refractory sequence is only 10 nt long. This finding also suggests that all mRNAs display fluctuating thermodynamic landscapes and that, in terms of thermodynamic features, RNAi is a very efficient antiviral system since there will always be sites susceptible to siRNAs.

Ancient mtDNA sequences in the human nuclear genome: A potential source of errors in identifying pathogenic mutations

Nuclear-localized mtDNA pseudogenes might explain a recent report describing a heteroplasmic mtDNA molecule containing five linked missense mutations dispersed over the contiguous mtDNA CO1 and CO2 genes in Alzheimer’s disease (AD) patients. To test this hypothesis, we have used the PCR primers utilized in the original report to amplify CO1 and CO2 sequences from two independent ρ° (mtDNA-less) cell lines. CO1 and CO2 sequences amplified from both of the ρ° cells, demonstrating that these sequences are also present in the human nuclear DNA. The nuclear pseudogene CO1 and CO2 sequences were then tested for each of the five “AD” missense mutations by restriction endonuclease site variant assays. All five mutations were found in the nuclear CO1 and CO2 PCR products from ρ° cells, but none were found in the PCR products obtained from cells with normal mtDNA. Moreover, when the overlapping nuclear CO1 and CO2 PCR products were cloned and sequenced, all five missense mutations were found, as well as a linked synonymous mutation. Unlike the findings in the original report, an additional 32 base substitutions were found, including two in adjacent tRNAs and a two base pair deletion in the CO2 gene. Phylogenetic analysis of the nuclear CO1 and CO2 sequences revealed that they diverged from modern human mtDNAs early in hominid evolution about 770,000 years before present. These data would be consistent with the interpretation that the missense mutations proposed to cause AD may be the product of ancient mtDNA variants preserved as nuclear pseudogenes.

Detection of unusual mutation within the VP1 region of different re-isolates of poliovirus Sabin vaccine

In the present study, a genomic analysis of full VP1 sequence region of 15 clinical re-isolates (14 healthy vaccinees and one bone marrow tumor patient) was conducted, aiming to the identification of mutations and to the assessment of their impact on virus fitness, providing also insights relevant with the natural evolution of Sabin strains. Clinical re-isolates were analyzed by RT-PCR, sequencing and computational analysis. Some re-isolates were characterized by an unusual mutational pattern in which non-synonymous mutations outnumbered the synonymous ones. Furthermore, the majority of amino-acid substitutions were located in the capsid exterior, specifically in N-Ags, near N-Ags and in the north rim of the canyon. Also mutations, which are well-known determinants of attenuation, were identified. The results of this study propose that some re-isolates are characterized by an evolutionary pattern in which non-synonymous mutations with a direct phenotypic impact on viral fitness are fixed in viral genomes, in spite of synonymous ones with no phenotypic impact on viral fitness. Results of the present retrospective characterization of Sabin clinical re-isolates, based on the full VP1 sequence, suggest that vaccine-derived viruses may make their way through narrow breaches and may evolve into transmissible pathogens even in adequately immunized populations. For this reason increased poliovirus laboratory surveillance should be permanent and full VP1 sequence analysis should be conducted even in isolates originating from healthy vaccinees.

A Nonsense Mutation in the IKBKG Gene in Mares with Incontinentia Pigmenti.

Ectodermal dysplasias (EDs) are a large and heterogeneous group of hereditary disorders characterized by abnormalities in structures of ectodermal origin. Incontinentia pigmenti (IP) is an ED characterized by skin lesions evolving over time, as well as dental, nail, and ocular abnormalities. Due to X-linked dominant inheritance IP symptoms can only be seen in female individuals while affected males die during development in utero. We observed a family of horses, in which several mares developed signs of a skin disorder reminiscent of human IP. Cutaneous manifestations in affected horses included the development of pruritic, exudative lesions soon after birth. These developed into wart-like lesions and areas of alopecia with occasional wooly hair re-growth. Affected horses also had streaks of darker and lighter coat coloration from birth. The observation that only females were affected together with a high number of spontaneous abortions suggested an X-linked dominant mechanism of transmission. Using next generation sequencing we sequenced the whole genome of one affected mare. We analyzed the sequence data for non-synonymous variants in candidate genes and found a heterozygous nonsense variant in the X-chromosomal IKBKG gene (c.184C>T; p.Arg62*). Mutations in IKBKG were previously reported to cause IP in humans and the homologous p.Arg62* variant has already been observed in a human IP patient. The comparative data thus strongly suggest that this is also the causative variant for the observed IP in horses. To our knowledge this is the first large animal model for IP.

A mutation in the SUV39H2 gene in Labrador Retrievers with hereditary nasal parakeratosis (HNPK) provides insights into the epigenetics of keratinocyte differentiation.

Hereditary nasal parakeratosis (HNPK), an inherited monogenic autosomal recessive skin disorder, leads to crusts and fissures on the nasal planum of Labrador Retrievers. We performed a genome-wide association study (GWAS) using 13 HNPK cases and 23 controls. We obtained a single strong association signal on chromosome 2 (p(raw) = 4.4×10⁻¹⁴). The analysis of shared haplotypes among the 13 cases defined a critical interval of 1.6 Mb with 25 predicted genes. We re-sequenced the genome of one case at 38× coverage and detected 3 non-synonymous variants in the critical interval with respect to the reference genome assembly. We genotyped these variants in larger cohorts of dogs and only one was perfectly associated with the HNPK phenotype in a cohort of more than 500 dogs. This candidate causative variant is a missense variant in the SUV39H2 gene encoding a histone 3 lysine 9 (H3K9) methyltransferase, which mediates chromatin silencing. The variant c.972T>G is predicted to change an evolutionary conserved asparagine into a lysine in the catalytically active domain of the enzyme (p.N324K). We further studied the histopathological alterations in the epidermis in vivo. Our data suggest that the HNPK phenotype is not caused by hyperproliferation, but rather delayed terminal differentiation of keratinocytes. Thus, our data provide evidence that SUV39H2 is involved in the epigenetic regulation of keratinocyte differentiation ensuring proper stratification and tight sealing of the mammalian epidermis.

A frameshift mutation in the cubilin gene (CUBN) in Border Collies with Imerslund-Gräsbeck syndrome (selective cobalamin malabsorption).

Imerslund-Gräsbeck syndrome (IGS) or selective cobalamin malabsorption has been described in humans and dogs. IGS occurs in Border Collies and is inherited as a monogenic autosomal recessive trait in this breed. Using 7 IGS cases and 7 non-affected controls we mapped the causative mutation by genome-wide association and homozygosity mapping to a 3.53 Mb interval on chromosome 2. We re-sequenced the genome of one affected dog at ∼10× coverage and detected 17 non-synonymous variants in the critical interval. Two of these non-synonymous variants were in the cubilin gene (CUBN), which is known to play an essential role in cobalamin uptake from the ileum. We tested these two CUBN variants for association with IGS in larger cohorts of dogs and found that only one of them was perfectly associated with the phenotype. This variant, a single base pair deletion (c.8392delC), is predicted to cause a frameshift and premature stop codon in the CUBN gene. The resulting mutant open reading frame is 821 codons shorter than the wildtype open reading frame (p.Q2798Rfs*3). Interestingly, we observed an additional nonsense mutation in the MRC1 gene encoding the mannose receptor, C type 1, which was in perfect linkage disequilibrium with the CUBN frameshift mutation. Based on our genetic data and the known role of CUBN for cobalamin uptake we conclude that the identified CUBN frameshift mutation is most likely causative for IGS in Border Collies.

A COL11A2 mutation in Labrador retrievers with mild disproportionate dwarfism.

We describe a mild form of disproportionate dwarfism in Labrador Retrievers, which is not associated with any obvious health problems such as secondary arthrosis. We designate this phenotype as skeletal dysplasia 2 (SD2). It is inherited as a monogenic autosomal recessive trait with incomplete penetrance primarily in working lines of the Labrador Retriever breed. Using 23 cases and 37 controls we mapped the causative mutation by genome-wide association and homozygosity mapping to a 4.44 Mb interval on chromosome 12. We re-sequenced the genome of one affected dog at 30x coverage and detected 92 non-synonymous variants in the critical interval. Only two of these variants, located in the lymphotoxin A (LTA) and collagen alpha-2(XI) chain gene (COL11A2), respectively, were perfectly associated with the trait. Previously described COL11A2 variants in humans or mice lead to skeletal dysplasias and/or deafness. The dog variant associated with disproportionate dwarfism, COL11A2:c.143G>C or p.R48P, probably has only a minor effect on collagen XI function, which might explain the comparatively mild phenotype seen in our study. The identification of this candidate causative mutation thus widens the known phenotypic spectrum of COL11A2 mutations. We speculate that non-pathogenic COL11A2 variants might even contribute to the heritable variation in height.

A Gly98Val mutation in the N-Myc downstream regulated gene 1 (NDRG1) in Alaskan Malamutes with polyneuropathy.

The first cases of early-onset progressive polyneuropathy appeared in the Alaskan Malamute population in Norway in the late 1970s. Affected dogs were of both sexes and were ambulatory paraparetic, progressing to non-ambulatory tetraparesis. On neurologic examination, affected dogs displayed predominantly laryngeal paresis, decreased postural reactions, decreased spinal reflexes and muscle atrophy. The disease was considered eradicated through breeding programmes but recently new cases have occurred in the Nordic countries and the USA. The N-myc downstream-regulated gene (NDRG1) is implicated in neuropathies with comparable symptoms or clinical signs both in humans and in Greyhound dogs. This gene was therefore considered a candidate gene for the polyneuropathy in Alaskan Malamutes. The coding sequence of the NDRG1 gene derived from one healthy and one affected Alaskan Malamute revealed a non-synonymous G>T mutation in exon 4 in the affected dog that causes a Gly98Val amino acid substitution. This substitution was categorized to be "probably damaging" to the protein function by PolyPhen2 (score: 1.000). Subsequently, 102 Alaskan Malamutes from the Nordic countries and the USA known to be either affected (n = 22), obligate carriers (n = 7) or healthy (n = 73) were genotyped for the SNP using TaqMan. All affected dogs had the T/T genotype, the obligate carriers had the G/T genotype and the healthy dogs had the G/G genotype except for 13 who had the G/T genotype. A protein alignment showed that residue 98 is conserved in mammals and also that the entire NDRG1 protein is highly conserved (94.7%) in mammals. We conclude that the G>T substitution is most likely the mutation that causes polyneuropathy in Alaskan Malamutes. Our characterization of a novel candidate causative mutation for polyneuropathy offers a new canine model that can provide further insight into pathobiology and therapy of human polyneuropathy. Furthermore, selection against this mutation can now be used to eliminate the disease in Alaskan Malamutes.

Mutation and association analysis of the PVR and PVRL2 genes in patients with non-syndromic cleft lip and palate.

Orofacial clefts (OFC; MIM 119530) are among the most common major birth defects. Here, we carried out mutation screening of the PVR and PVRL2 genes, which are both located at an OFC linkage region at 19q13 (OFC3) and are closely related to PVRL1, which has been associated with both syndromic and non-syndromic cleft lip and palate (nsCLP). We screened a total of 73 nsCLP patients and 105 non-cleft controls from the USA for variants in PVR and PVRL2, including all exons and encompassing all isoforms. We identified four variants in PVR and five in PVRL2. One non-synonymous PVR variant, A67T, was more frequent among nsCLP patients than among normal controls, but this difference did not achieve statistical significance.

A mutation in the FAM83G gene in dogs with hereditary footpad hyperkeratosis (HFH)

Hereditary footpad hyperkeratosis (HFH) represents a palmoplantar hyperkeratosis, which is inherited as a monogenic autosomal recessive trait in several dog breeds, such as e.g. Kromfohrländer and Irish Terriers. We performed genome-wide association studies (GWAS) in both breeds. In Kromfohrländer we obtained a single strong association signal on chromosome 5 (p(raw) = 1.0×10(-13)) using 13 HFH cases and 29 controls. The association signal replicated in an independent cohort of Irish Terriers with 10 cases and 21 controls (p(raw) = 6.9×10(-10)). The analysis of shared haplotypes among the combined Kromfohrländer and Irish Terrier cases defined a critical interval of 611 kb with 13 predicted genes. We re-sequenced the genome of one affected Kromfohrländer at 23.5× coverage. The comparison of the sequence data with 46 genomes of non-affected dogs from other breeds revealed a single private non-synonymous variant in the critical interval with respect to the reference genome assembly. The variant is a missense variant (c.155G>C) in the FAM83G gene encoding a protein with largely unknown function. It is predicted to change an evolutionary conserved arginine into a proline residue (p.R52P). We genotyped this variant in a larger cohort of dogs and found perfect association with the HFH phenotype. We further studied the clinical and histopathological alterations in the epidermis in vivo. Affected dogs show a moderate to severe orthokeratotic hyperplasia of the palmoplantar epidermis. Thus, our data provide the first evidence that FAM83G has an essential role for maintaining the integrity of the palmoplantar epidermis.