Improved prediction of critical residues for protein function based on network and phylogenetic analyses

Background  

Phylogenetic approaches are commonly used to predict which amino acid residues are critical to the function of a given protein. However, such approaches display inherent limitations, such as the requirement for identification of multiple homologues of the protein under consideration. Therefore, complementary or alternative approaches for the prediction of critical residues would be desirable. Network analyses have been used in the modelling of many complex biological systems, but only very recently have they been used to predict critical residues from a protein's three-dimensional structure. Here we compare a couple of phylogenetic approaches to several different network-based methods for the prediction of critical residues, and show that a combination of one phylogenetic method and one network-based method is superior to other methods previously employed.     

Phenylketonuria in U.S. blacks: molecular analysis of the phenylalanine hydroxylase gene

We investigated the frequency, origin, and molecular basis of phenylketonuria (PKU) in U.S. blacks. On the basis of 10 years of Maryland newborn-screening data, we found the frequency to be 1/50,000, or one-third that in whites. We performed haplotype analysis of the phenylalanine hydroxylase (PAH) gene of 36 U.S. blacks, 16 from individuals with classical PKU and 20 from controls. In blacks, 20% of wild-type PAH alleles have a common Caucasian haplotype (i.e., haplotype 1), whereas 80% had a variety of haplotypes, all rare in Caucasians and Asians. One of these, haplotype 15, accounted for a large fraction (30%). Among black mutant PAH alleles, 20% have a haplotype (i.e., either haplotype 1 or haplotype 4) common in Caucasians; 40% have a haplotype rare in Caucasians and Asians, and 40% have one of two previously undescribed haplotypes. Both can be derived from known haplotypes by a single event. One of these haplotypes is characterized by a new MspI restriction site, located in intron 8, which was present in five of 16 black mutant alleles but was not present in 60 U.S. black control, 20 U.S. Caucasian control, or 20 Caucasian mutant PAH alleles. Sequence analysis of DNA from a single individual, homozygous for the new MspI associated haplotype, shows homozygosity for a C----T transition at nucleotide 896 in exon 7 of the PAH cDNA, resulting in the conversion of leucine 255 to serine (L255S).     

Increasing calling accuracy,coverage, and read-depth in sequence data by the use of haplotype blocks

High-throughput genotyping of large numbers of lines remains a key challenge in plant genetics, requiring geneticists and breeders to find a balance between data quality and the number of genotyped lines under a variety of different existing genotyping technologies when resources are limited. In this work, we are proposing a new imputation pipeline (“HBimpute”) that can be used to generate high-quality genomic data from low read-depth whole-genome-sequence data. The key idea of the pipeline is the use of haplotype blocks from the software HaploBlocker to identify locally similar lines and subsequently use the reads of all locally similar lines in the variant calling for a specific line. The effectiveness of the pipeline is showcased on a dataset of 321 doubled haploid lines of a European maize landrace, which were sequenced at 0.5X read-depth. The overall imputing error rates are cut in half compared to state-of-the-art software like BEAGLE and STITCH, while the average read-depth is increased to 83X, thus enabling the calling of copy number variation. The usefulness of the obtained imputed data panel is further evaluated by comparing the performance of sequence data in common breeding applications to that of genomic data generated with a genotyping array. For both genome-wide association studies and genomic prediction, results are on par or even slightly better than results obtained with high-density array data (600k). In particular for genomic prediction, we observe slightly higher data quality for the sequence data compared to the 600k array in the form of higher prediction accuracies. This occurred specifically when reducing the data panel to the set of overlapping markers between sequence and array, indicating that sequencing data can benefit from the same marker ascertainment as used in the array process to increase the quality and usability of genomic data.     

Serotypes and Clonal Diversity of Streptococcus pneumoniae Causing Invasive Disease in the Era of PCV13 in Catalonia,Spain

The aim of this study was to study the serotypes and clonal diversity of pneumococci causing invasive pneumococcal disease in Catalonia, Spain, in the era of 13-valent pneumococcal conjugate vaccine (PCV13). In our region, this vaccine is only available in the private market and it is estimated a PCV13 vaccine coverage around 55% in children. A total of 1551 pneumococcal invasive isolates received between 2010 and 2013 in the Molecular Microbiology Department at Hospital Sant Joan de Déu, Barcelona, were included. Fifty-two serotypes and 249 clonal types—defined by MLST—were identified. The most common serotypes were serotype 1 (n = 182; 11.7%), 3 (n = 145; 9.3%), 19A (n = 137; 8.8%) and 7F (n = 122; 7.9%). Serotype 14 was the third most frequent serotype in children < 2 years (15 of 159 isolates). PCV7 serotypes maintained their proportion along the period of study, 16.6% in 2010 to 13.4% in 2013, whereas there was a significant proportional decrease in PCV13 serotypes, 65.3% in 2010 to 48.9% in 2013 (p<0.01). This decrease was mainly attributable to serotypes 19A and 7F. Serotype 12F achieved the third position in 2013 (n = 22, 6.4%). The most frequent clonal types found were ST306 (n = 154, 9.9%), ST191 (n = 111, 7.2%), ST989 (n = 85, 5.5%) and ST180 (n = 80, 5.2%). Despite their decrease, PCV13 serotypes continue to be a major cause of disease in Spain. These results emphasize the need for complete PCV13 vaccination.     

Improvement of culture conditions to overproduce β-galactosidase from Escherichia coli in Bacillus subtilis

The effect of some culture variables in the production of β-galactosidase from Escherichia coli in Bacillus subtilis was evaluated. The lacZ gene was expressed in B. subtilis using the regulatory region of the subtilisin gene aprE. The host contained also the hpr2 and degU32 mutations, which are known to overexpress the aprE gene. We found that, when this overproducing B. subtilis strain was grown in mineral medium supplemented with glucose (MMG), β-galactosidase production was partially growth-associated, as 40%–60% of the maximum enzyme activity was produced before the onset of the stationary phase. In contrast, when a complex medium was used, β-galactosidase was produced only at low levels during vegetative growth, whereas it accumulated to high levels during early stationary phase. Compared with the results obtained in complex media, a 20% increase in specific β-galactosidase activity in MMG supplemented with 11.6 g/l glucose was obtained. On the 1-l fermenter scale, a threefold increase in volumetric β-galactosidase activity was obtained when the glucose concentration was varied from 11 g/l to 26 g/l. In addition, glucose feeding during the stationary phase resulted in a twofold increase in volumetric enzyme activity as cellular lysis was prevented. Finally, we showed that oxygen uptake and carbon dioxide evolution rates can be used for on-line determination of the onset of stationary phase, glucose depletion and biomass concentration. Received: 18 April 1996 / Received revision: 27 August 1996 / Accepted: 6 September 1996     

Operation of an aerobic granular pilot scale SBR plant to treat swine slurry

A pilot-scale Sequencing Batch Reactor was operated during 307 days in order to treat swine slurry characterized by its high variable composition: organic and nitrogen applied loading rates and C/N ratio were 1.4–6.3 kg COD_s/(m³ d), 0.5–2.5 kg N/(m³ d) and 1.9–9.4 g COD_s/(g N), respectively. Aerobic granules successfully developed in the reactor and their physical properties remained rather stable despite the feeding composition variability. Organic and ammonia removal efficiency reached 61–73% and 56–77%, respectively, however ammonia was mainly oxidized to nitrite. The reactor had a good biomass retention capacity to select for granular biomass. However, its efficiency to retain the solids present in the feeding was low. Aerobic granulation in SBR systems appears as an interesting alternative to treat slurry in small livestock facilities where the implementation of anaerobic digestion systems is not a feasible option or the removal of nitrogenous compounds is required.     

Virological Efficacy in Cerebrospinal Fluid and Neurocognitive Status in Patients with Long-Term Monotherapy Based on Lopinavir/Ritonavir: An Exploratory Study

Background

Data on suppression of HIV replication in the CNS and on the subsequent risk of neurocognitive impairment using monotherapy with boosted protease inhibitors are limited.

Methods

Ours was an exploratory cross-sectional study in patients on lopinavir/ritonavir-based monotherapy (LPV/r-MT) or standard triple therapy (LPV/r-ART) for at least 96 weeks who maintained a plasma viral load <50 copies/mL. HIV-1 RNA in CSF was determined by HIV-1 SuperLow assay (lower limit of detection, 1 copy/mL). Neurocognitive functioning was assessed using a recommended battery of neuropsychological tests covering 7 areas. Neurocognitive impairment (NCI) was determined and also a global deficit score (GDS) for study comparisons.

Results

Seventeen patients on LPV/r-MT and 17 on LPV/r-ART were included. Fourteen (82.4%) patients on LPV/r-MT and 16 (94.1%) on LPV/r-ART had HIV-1 RNA <1 copy/mL in CSF (p = 0.601). NCI was observed in 7 patients on LPV/r-MT and in 10 on LPV/r-ART (41% vs 59%; p = 0.494). Mean (SD) GDS was 0.22 (0.20) in patients on LPV/r-MT and 0.47 (0.34) in those on LPV/r-ART (p = 0.012).

Conclusions

Suppression of HIV in CSF is similar in individuals with durable plasma HIV-1 RNA suppression who are receiving LPV/r-MT or LPV/r-ART for at least 96 weeks. Findings for HIV-1 replication in CSF and neurocognitive status indicate that this strategy seems to be safe for CNS functioning.     

Reciprocal Regulation of Epileptiform Neuronal Oscillations and Electrical Synapses in the Rat Hippocampus

Gap junction (GJ) channels have been recognized as an important mechanism for synchronizing neuronal networks. Herein, we investigated the participation of GJ channels in the pilocarpine-induced status epilepticus (SE) by analyzing electrophysiological activity following the blockade of connexins (Cx)-mediated communication. In addition, we examined the regulation of gene expression, protein levels, phosphorylation profile and distribution of neuronal Cx36, Cx45 and glial Cx43 in the rat hippocampus during the acute and latent periods. Electrophysiological recordings revealed that the GJ blockade anticipates the occurrence of low voltage oscillations and promotes a marked reduction of power in all analyzed frequencies.Cx36 gene expression and protein levels remained stable in acute and latent periods, whereas upregulation of Cx45 gene expression and protein redistribution were detected in the latent period. We also observed upregulation of Cx43 mRNA levels followed by changes in the phosphorylation profile and protein accumulation. Taken together, our results indisputably revealed that GJ communication participates in the epileptiform activity induced by pilocarpine. Moreover, considering that specific Cxs undergo alterations through acute and latent periods, this study indicates that the control of GJ communication may represent a focus in reliable anti-epileptogenic strategies.     

Variation in estimated recombination rates across human populations

Recently it has been reported that recombination hotspots appear to be highly variable between humans and chimpanzees, and there is evidence for between-person variability in hotspots, and evolutionary transience. To understand the nature of variation in human recombination rates, it is important to describe patterns of variability across populations. Direct measurement of recombination rates remains infeasible on a large scale, and population-genetic approaches can be imprecise, and are affected by demographic history. Reports to date have suggested broad similarity in recombination rates at large genomic scales and across human populations. Here, we examine recombination rate estimates at a finer population and genomic scale: 28 worldwide populations and 107 SNPs in a 1 Mb stretch of chromosome 22q. We employ analysis of variance of recombination rate estimates, corrected for differences in effective population size using genome-wide microsatellite mutation rate estimates. We find substantial variation in fine-scale rates between populations, but reduced variation within continental groups. All effects examined (SNP-pair, region, population and interactions) were highly significant. Adjustment for effective population size made little difference to the conclusions. Observed hotspots tended to be conserved across populations, albeit at varying intensities. This holds particularly for populations from the same region, and also to a considerable degree across geographical regions. However, some hotspots appear to be population-specific. Several results from studies on the population history of humans are in accordance with our analysis. Our results suggest that between-population variation in DNA sequences may underly recombination rate variation.