Reconstruction of the Evolutionary Dynamics of A(H3N2) Influenza Viruses Circulating in Italy from 2004 to 2012

Background

Influenza A viruses are characterised by their rapid evolution, and the appearance of point mutations in the viral hemagglutinin (HA) domain causes seasonal epidemics. The A(H3N2) virus has higher mutation rate than the A(H1N1) virus. The aim of this study was to reconstruct the evolutionary dynamics of the A(H3N2) viruses circulating in Italy between 2004 and 2012 in the light of the forces driving viral evolution.

Methods

Phylodinamic analyses were made using a Bayesian method, and codon-specific positive selection acting on the HA coding sequence was evaluated.

Results

Global and local phylogenetic analyses showed that the Italian strains collected between 2004 and 2012 grouped into five significant Italian clades that included viral sequences circulating in different epidemic seasons. The time of the most recent common ancestor (tMRCA) of the tree root was between May and December 2003. The tMRCA estimates of the major clades suggest that the origin of a new viral strain precedes the effective circulation of the strain in the Italian population by 6–31 months, thus supporting a central role of global migration in seeding the epidemics in Italy. The study of selection pressure showed that four codons were under positive selection, three of which were located in antigenic sites. Analysis of population dynamics showed the alternation of periods of exponential growth followed by a decrease in the effective number of infections corresponding to epidemic and inter-epidemic seasons.

Conclusions

Our analyses suggest that a complex interaction between the immune status of the population, migrations, and a few selective sweeps drive the influenza A(H3N2) virus evolution. Our findings suggest the possibility of the year-round survival of local strains even in temperate zones, a hypothesis that warrants further investigation.     

Discovery of Candidate Disease Genes in ENU–Induced Mouse Mutants by Large-Scale Sequencing,Including a Splice-Site Mutation in Nucleoredoxin

An accurate and precisely annotated genome assembly is a fundamental requirement for functional genomic analysis. Here, the complete DNA sequence and gene annotation of mouse Chromosome 11 was used to test the efficacy of large-scale sequencing for mutation identification. We re-sequenced the 14,000 annotated exons and boundaries from over 900 genes in 41 recessive mutant mouse lines that were isolated in an N-ethyl-N-nitrosourea (ENU) mutation screen targeted to mouse Chromosome 11. Fifty-nine sequence variants were identified in 55 genes from 31 mutant lines. 39% of the lesions lie in coding sequences and create primarily missense mutations. The other 61% lie in noncoding regions, many of them in highly conserved sequences. A lesion in the perinatal lethal line l11Jus13 alters a consensus splice site of nucleoredoxin (Nxn), inserting 10 amino acids into the resulting protein. We conclude that point mutations can be accurately and sensitively recovered by large-scale sequencing, and that conserved noncoding regions should be included for disease mutation identification. Only seven of the candidate genes we report have been previously targeted by mutation in mice or rats, showing that despite ongoing efforts to functionally annotate genes in the mammalian genome, an enormous gap remains between phenotype and function. Our data show that the classical positional mapping approach of disease mutation identification can be extended to large target regions using high-throughput sequencing.     

Anisotropic elastic network modeling of entire microtubules

Microtubules are supramolecular structures that make up the cytoskeleton and strongly affect the mechanical properties of the cell. Within the cytoskeleton filaments, the microtubule (MT) exhibits by far the highest bending stiffness. Bending stiffness depends on the mechanical properties and intermolecular interactions of the tubulin dimers (the MT building blocks). Computational molecular modeling has the potential for obtaining quantitative insights into this area. However, to our knowledge, standard molecular modeling techniques, such as molecular dynamics (MD) and normal mode analysis (NMA), are not yet able to simulate large molecular structures like the MTs; in fact, their possibilities are normally limited to much smaller protein complexes. In this work, we developed a multiscale approach by merging the modeling contribution from MD and NMA. In particular, MD simulations were used to refine the molecular conformation and arrangement of the tubulin dimers inside the MT lattice. Subsequently, NMA was used to investigate the vibrational properties of MTs modeled as an elastic network. The coarse-grain model here developed can describe systems of hundreds of interacting tubulin monomers (corresponding to up to 1,000,000 atoms). In particular, we were able to simulate coarse-grain models of entire MTs, with lengths up to 350 nm. A quantitative mechanical investigation was performed; from the bending and stretching modes, we estimated MT macroscopic properties such as bending stiffness, Young modulus, and persistence length, thus allowing a direct comparison with experimental data.     

A novel haplotype within C-reactive protein gene influences CRP levels and coronary heart disease risk in Northwest Indians

According to several epidemiological and clinical studies, the concentration of C-reactive protein (CRP) in blood is associated with the risk of coronary heart disease (CHD). However, these studies are limited in high incidence and prevalence area of North-West India. The present case control study investigated the contribution of three relevant CRP single nucleotide polymorphisms: ?717A>G located in the promoter region (rs2794521), +1059G>C on exon2 (rs1800947) and +1444C>T in the 3′ UTR (rs1130864) in 180 angiographically verified CHD cases and 175 control subjects. Minor allele frequencies (G, C and T) of rs2794521, rs1800947 and rs1130864 are observed to be 21.1, 11.7, 29.4 and 11.4, 10.0, 19.7 % in CHD cases and controls respectively. AA genotype of ?717A>G and TT genotype of +1444C>T were significantly associated (P = 0.02 & 0.03 respectively) with the risk of CHD whereas, +1059G and +1444T were found to be strongly related (P = 0.023 & P = 0.008 respectively) with multivariable adjusted CRP levels. AGT Haplotype was significantly associated with the adjusted CRP levels (P < 0.05). Disease association analysis revealed that haplotype AGT influences CHD risk (OR 2.4, 95 % CI 1.23–4.84, P = 0.006) which exacerbates after correcting the confounding effects of risk variables (OR 2.5, 95 % CI 1.27–4.99, P = 0.004). With the global index of Akaike information criterion, it has been observed that the carrying each single unit of this susceptibility haplotype increases CHD risk by a value of 2.41 ± 0.439 (β ± SE) in the recessive mode.     

Ultra-deep sequencing reveals the microRNA expression pattern of the human stomach

Background

While microRNAs (miRNAs) play important roles in tissue differentiation and in maintaining basal physiology, little is known about the miRNA expression levels in stomach tissue. Alterations in the miRNA profile can lead to cell deregulation, which can induce neoplasia.

Methodology/Principal Findings

A small RNA library of stomach tissue was sequenced using high-throughput SOLiD sequencing technology. We obtained 261,274 quality reads with perfect matches to the human miRnome, and 42% of known miRNAs were identified. Digital Gene Expression profiling (DGE) was performed based on read abundance and showed that fifteen miRNAs were highly expressed in gastric tissue. Subsequently, the expression of these miRNAs was validated in 10 healthy individuals by RT-PCR showed a significant correlation of 83.97% (P<0.05). Six miRNAs showed a low variable pattern of expression (miR-29b, miR-29c, miR-19b, miR-31, miR-148a, miR-451) and could be considered part of the expression pattern of the healthy gastric tissue.

Conclusions/Significance

This study aimed to validate normal miRNA profiles of human gastric tissue to establish a reference profile for healthy individuals. Determining the regulatory processes acting in the stomach will be important in the fight against gastric cancer, which is the second-leading cause of cancer mortality worldwide.     

eNOS Activation by HDL Is Impaired in Genetic CETP Deficiency

Mutations in the CETP gene resulting in defective CETP activity have been shown to cause remarkable elevations of plasma HDL-C levels, with the accumulation in plasma of large, buoyant HDL particles enriched in apolipoprotein E. Genetic CETP deficiency thus represents a unique tool to evaluate how structural alterations of HDL impact on HDL atheroprotective functions. Aim of the present study was to assess the ability of HDL obtained from CETP-deficient subjects to protect endothelial cells from the development of endothelial dysfunction. HDL isolated from one homozygous and seven heterozygous carriers of CETP null mutations were evaluated for their ability to down-regulate cytokine-induced cell adhesion molecule expression and to promote NO production in cultured endothelial cells. When compared at the same protein concentration, HDL and HDL3 from carriers proved to be as effective as control HDL and HDL3 in down-regulating cytokine-induced VCAM-1, while carrier HDL2 were more effective than control HDL2 in inhibiting VCAM-1 expression. On the other hand, HDL and HDL fractions from carriers of CETP deficiency were significantly less effective than control HDL and HDL fractions in stimulating NO production, due to a reduced eNOS activating capacity, likely because of a reduced S1P content. In conclusion, the present findings support the notion that genetic CETP deficiency, by affecting HDL particle structure, impacts on HDL vasculoprotective functions. Understanding of these effects might be important for predicting the outcomes of pharmacological CETP inhibition.     

Pallial origin of basal forebrain cholinergic neurons in the nucleus basalis of Meynert and horizontal limb of the diagonal band nucleus

The majority of the cortical cholinergic innervation implicated in attention and memory originates in the nucleus basalis of Meynert and in the horizontal limb of the diagonal band nucleus of the basal prosencephalon. Functional alterations in this system give rise to neuropsychiatric disorders as well as to the cognitive alterations described in Parkinson and Alzheimer's diseases. Despite the functional importance of these basal forebrain cholinergic neurons very little is known about their origin and development. Previous studies suggest that they originate in the medial ganglionic eminence of the telencephalic subpallium; however, our results identified Tbr1-expressing, reelin-positive neurons migrating from the ventral pallium to the subpallium that differentiate into cholinergic neurons in the basal forebrain nuclei projecting to the cortex. Experiments with Tbr1 knockout mice, which lack ventropallial structures, confirmed the pallial origin of cholinergic neurons in Meynert and horizontal diagonal band nuclei. Also, we demonstrate that Fgf8 signaling in the telencephalic midline attracts these neurons from the pallium to follow a tangential migratory route towards the basal forebrain.     

Purification and characterization of the enantioselective esterase from Kluyveromyces marxianus CBS 1553

An intracellular esterase from the yeast Kluyveromyces marxianus CBS 1553 with interesting enantioselective hydrolytic activity towards racemic esters of 1,2-O-isopropylidene glycerol (IPG) was purified and characterized. Optimal culture conditions for the obtainment of the enantioselective esterase on a 5 l-fermentation scale were investigated. Two esterase activities (EST1 and EST2) in the crude cell extract were identified by native PAGE with specific activity staining and separated from each other by anion-exchange chromatography. EST1 showed higher activity and enantioselectivity than EST2 in the resolution of racemic IPG acetate and was further purified by hydrophobic interaction chromatography and preparative electrophoresis (final specific activity approximately = 300 U mg(-1), showing a main protein band with a molecular mass of 29 kDa. EST1 showed optimal activity between pH 8.0 and 10.0 and was stable in the pH range 7.0-10.0. Moreover, it was rather thermostable and active up to 80 degrees C, and retained most of its activity in the presence of 15% (v/v) of various organic solvents. The enzyme showed similar Vmax in the hydrolysis of the acetate esters of IPG, whereas the Km value towards (S)-IPG acetate was significantly lower than the one towards the (R)-enantiomer (5.3 and 70 microM, respectively). Finally, comparison of EST1 activity in the presence of different glycerol esters and synthetic substrates with different chain lengths showed a strong preference of this biocatalyst for short-chain substrates.