Assessing Mitochondrial DNA Variation and Copy Number in Lymphocytes of ~2,000 Sardinians Using Tailored Sequencing Analysis Tools

DNA sequencing identifies common and rare genetic variants for association studies, but studies typically focus on variants in nuclear DNA and ignore the mitochondrial genome. In fact, analyzing variants in mitochondrial DNA (mtDNA) sequences presents special problems, which we resolve here with a general solution for the analysis of mtDNA in next-generation sequencing studies. The new program package comprises 1) an algorithm designed to identify mtDNA variants (i.e., homoplasmies and heteroplasmies), incorporating sequencing error rates at each base in a likelihood calculation and allowing allele fractions at a variant site to differ across individuals; and 2) an estimation of mtDNA copy number in a cell directly from whole-genome sequencing data. We also apply the methods to DNA sequence from lymphocytes of ~2,000 SardiNIA Project participants. As expected, mothers and offspring share all homoplasmies but a lesser proportion of heteroplasmies. Both homoplasmies and heteroplasmies show 5-fold higher transition/transversion ratios than variants in nuclear DNA. Also, heteroplasmy increases with age, though on average only ~1 heteroplasmy reaches the 4% level between ages 20 and 90. In addition, we find that mtDNA copy number averages ~110 copies/lymphocyte and is ~54% heritable, implying substantial genetic regulation of the level of mtDNA. Copy numbers also decrease modestly but significantly with age, and females on average have significantly more copies than males. The mtDNA copy numbers are significantly associated with waist circumference (p-value = 0.0031) and waist-hip ratio (p-value = 2.4×10^-5), but not with body mass index, indicating an association with central fat distribution. To our knowledge, this is the largest population analysis to date of mtDNA dynamics, revealing the age-imposed increase in heteroplasmy, the relatively high heritability of copy number, and the association of copy number with metabolic traits.     

Disruption of the ribosomal P complex leads to stress-induced autophagy

The human ribosomal P complex, which consists of the acidic ribosomal P proteins RPLP0, RPLP1, and RPLP2 (RPLP proteins), recruits translational factors, facilitating protein synthesis. Recently, we showed that overexpression of RPLP1 immortalizes primary cells and contributes to transformation. Moreover, RPLP proteins are overexpressed in human cancer, with the highest incidence in breast carcinomas. It is thought that disruption of the P complex would directly affect protein synthesis, causing cell growth arrest and eventually apoptosis. Here, we report a distinct mechanism by which cancer cells undergo cell cycle arrest and induced autophagy when RPLP proteins are downregulated. We found that absence of RPLP0, RPLP1, or RPLP2 resulted in reactive oxygen species (ROS) accumulation and MAPK1/ERK2 signaling pathway activation. Moreover, ROS generation led to endoplasmic reticulum (ER) stress that involved the EIF2AK3/PERK-EIF2S1/eIF2α-EIF2S2-EIF2S3-ATF4/ATF-4- and ATF6/ATF-6-dependent arms of the unfolded protein response (UPR). RPLP protein-deficient cells treated with autophagy inhibitors experienced apoptotic cell death as an alternative to autophagy. Strikingly, antioxidant treatment prevented UPR activation and autophagy while restoring the proliferative capacity of these cells. Our results indicate that ROS are a critical signal generated by disruption of the P complex that causes a cellular response that follows a sequential order: first ROS, then ER stress/UPR activation, and finally autophagy. Importantly, inhibition of the first step alone is able to restore the proliferative capacity of the cells, preventing UPR activation and autophagy. Overall, our results support a role for autophagy as a survival mechanism in response to stress due to RPLP protein deficiency.     

Effects of particulate matter exposure on blood 5-hydroxymethylation: results from the Beijing truck driver air pollution study

Previous studies have reported epigenetic changes induced by environmental exposures. However, previous investigations did not distinguish 5-methylcytosine (5mC) from a similar oxidative form with opposite functions, 5-hydroxymethylcytosine (5hmC). Here, we measured blood DNA global 5mC and 5hmC by ELISA and used adjusted mixed-effects regression models to evaluate the effects of ambient PM₁₀ and personal PM_2.5 and its elemental components—black carbon (BC), aluminum (Al), calcium (Ca), potassium (K), iron (Fe), sulfur (S), silicon (Si), titanium (Ti), and zinc (Zn)—on blood global 5mC and 5hmC levels. The study was conducted in 60 truck drivers and 60 office workers in Beijing, China from The Beijing Truck Driver Air Pollution Study at 2 exams separated by one to 2 weeks. Blood 5hmC level (0.08%) was ∼83-fold lower than 5mC (6.61%). An inter-quartile range (IQR) increase in same-day PM₁₀ was associated with increases in 5hmC of 26.1% in office workers (P = 0.004), 20.2% in truck drivers (P = 0.014), and 21.9% in all participants combined (P < 0.001). PM₁₀ effects on 5hmC were increasingly stronger when averaged over 4, 7, and 14 d preceding assessment (up to 132.6% for the 14-d average in all participants, P < 0.001). PM₁₀ effects were also significant after controlling for multiple testing (family-wise error rate; FWER < 0.05). 5hmC was not correlated with personal measures of PM_2.5 and elemental components (FWER > 0.05). 5mC showed no correlations with PM₁₀, PM_2.5, and elemental components measures (FWER > 0.05). Our study suggests that exposure to ambient PM₁₀ affects 5hmC over time, but not 5mC. This finding demonstrates the need to differentiate 5hmC and 5mC in environmental studies of DNA methylation.     

Post‐translational modifications of transthyretin affect the triiodonine‐binding potential

Transthyretin (TTR) is a visceral protein, which facilitates the transport of thyroid hormones in blood and cerebrospinal fluid. The homotetrameric structure of TTR enables the simultaneous binding of two thyroid hormones per molecule. Each TTR subunit provides a single cysteine residue (Cys₁₀), which is frequently affected by oxidative post‐translational modifications. As Cys₁₀ is part of the thyroid hormone‐binding channel within the TTR molecule, PTM of Cys₁₀ may influence the binding of thyroid hormones. Therefore, we analysed the effects of Cys₁₀ modification with sulphonic acid, cysteine, cysteinylglycine and glutathione on binding of triiodothyronine (T3) by molecular modelling. Furthermore, we determined the PTM pattern of TTR in serum of patients with thyroid disease by immunoprecipitation and mass spectrometry to evaluate this association in vivo. The in silico assays demonstrated that oxidative PTM of TTR resulted in substantial reorganization of the intramolecular interactions and also affected the binding of T3 in a chemotype‐ and site‐specific manner with S‐glutathionylation as the most potent modulator of T3 binding. These findings were supported by the in vivo results, which indicated thyroid function‐specific patterns of TTR with a substantial decrease in S‐sulphonated, S‐cysteinylglycinated and S‐glutathionylated TTR in hypothyroid patients. In conclusion, this study provides evidence that oxidative modifications of Cys₁₀ seem to affect binding of T3 to TTR probably because of the introduction of a sterical hindrance and induction of conformational changes. As oxidative modifications can be dynamically regulated, this may represent a sensitive mechanism to adjust thyroid hormone availability.     

The Structure of the T190M Mutant of Murine α-Dystroglycan at High Resolution: Insight into the Molecular Basis of a Primary Dystroglycanopathy

The severe dystroglycanopathy known as a form of limb-girdle muscular dystrophy (LGMD2P) is an autosomal recessive disease caused by the point mutation T192M in α-dystroglycan. Functional expression analysis in vitro and in vivo indicated that the mutation was responsible for a decrease in posttranslational glycosylation of dystroglycan, eventually interfering with its extracellular-matrix receptor function and laminin binding in skeletal muscle and brain. The X-ray crystal structure of the missense variant T190M of the murine N-terminal domain of α-dystroglycan (50-313) has been determined, and showed an overall topology (Ig-like domain followed by a basket-shaped domain reminiscent of the small subunit ribosomal protein S6) very similar to that of the wild-type structure. The crystallographic analysis revealed a change of the conformation assumed by the highly flexible loop encompassing residues 159–180. Moreover, a solvent shell reorganization around Met190 affects the interaction between the B1–B5 anti-parallel strands forming part of the floor of the basket-shaped domain, with likely repercussions on the folding stability of the protein domain(s) and on the overall molecular flexibility. Chemical denaturation and limited proteolysis experiments point to a decreased stability of the T190M variant with respect to its wild-type counterpart. This mutation may render the entire L-shaped protein architecture less flexible. The overall reduced flexibility and stability may affect the functional properties of α-dystroglycan via negatively influencing its binding behavior to factors needed for dystroglycan maturation, and may lay the molecular basis of the T190M-driven primary dystroglycanopathy.     

Enhanced Amphiphilic Profile of a Short β-Stranded Peptide Improves Its Antimicrobial Activity

SB056 is a novel semi-synthetic antimicrobial peptide with a dimeric dendrimer scaffold. Active against both Gram-negative and -positive bacteria, its mechanism has been attributed to a disruption of bacterial membranes. The branched peptide was shown to assume a β-stranded conformation in a lipidic environment. Here, we report on a rational modification of the original, empirically derived linear peptide sequence [WKKIRVRLSA-NH₂, SB056-lin]. We interchanged the first two residues [KWKIRVRLSA-NH₂, β-SB056-lin] to enhance the amphipathic profile, in the hope that a more regular β-strand would lead to a better antimicrobial performance. MIC values confirmed that an enhanced amphiphilic profile indeed significantly increases activity against both Gram-positive and -negative strains. The membrane binding affinity of both peptides, measured by tryptophan fluorescence, increased with an increasing ratio of negatively charged/zwitterionic lipids. Remarkably, β-SB056-lin showed considerable binding even to purely zwitterionic membranes, unlike the original sequence, indicating that besides electrostatic attraction also the amphipathicity of the peptide structure plays a fundamental role in binding, by stabilizing the bound state. Synchrotron radiation circular dichroism and solid-state ¹⁹F-NMR were used to characterize and compare the conformation and mobility of the membrane bound peptides. Both SB056-lin and β-SB056-lin adopt a β-stranded conformation upon binding POPC vesicles, but the former maintains an intrinsic structural disorder that also affects its aggregation tendency. Upon introducing some anionic POPG into the POPC matrix, the sequence-optimized β-SB056-lin forms well-ordered β-strands once electro-neutrality is approached, and it aggregates into more extended β-sheets as the concentration of anionic lipids in the bilayer is raised. The enhanced antimicrobial activity of the analogue correlates with the formation of these extended β-sheets, which also leads to a dramatic alteration of membrane integrity as shown by ³¹P-NMR. These findings are generally relevant for the design and optimization of other membrane-active antimicrobial peptides that can fold into amphipathic β-strands.     

Spectral Properties of Single Gold Nanoparticles in Close Proximity to Biological Fluorophores Excited by 2-Photon Excitation

Metallic nanoparticles (NPs) are able to modify the excitation and emission rates (plasmonic enhancement) of fluorescent molecules in their close proximity. In this work, we measured the emission spectra of 20 nm Gold Nanoparticles (AuNPs) fixed on a glass surface submerged in a solution of different fluorophores using a spectral camera and 2-photon excitation. While on the glass surface, we observed the presence in the emission at least 3 components: i) second harmonic signal (SHG), ii) a broad emission from AuNPS and iii) fluorescence arising from fluorophores nearby. When on the glass surface, we found that the 3 spectral components have different relative intensities when the incident direction of linear polarization was changed indicating different physical origins for these components. Then we measured by fluctuation correlation spectroscopy (FCS) the scattering and fluorescence signal of the particles alone and in a solution of 100 nM EGFP using the spectral camera or measuring the scattering and fluorescence from the particles. We observed occasional fluorescence bursts when in the suspension we added fluorescent proteins. The spectrum of these burst was devoid of the SHG and of the broad emission in contrast to the signal collected from the gold nanoparticles on the glass surface. Instead we found that the spectrum during the burst corresponded closely to the spectrum of the fluorescent protein. An additional control was obtained by measuring the cross-correlation between the reflection from the particles and the fluorescence arising from EGFP both excited at 488 nm. We found a very weak cross-correlation between the AuNPs and the fluorescence confirming that the burst originate from a few particles with a fluorescence signal.