Effective Size of Populations under Selection

Equations to approximate the effective size (N(e)) of populations under continued selection are obtained that include the possibility of partial full-sib mating and other systems such as assortative mating. The general equation for the case of equal number of sexes and constant number of breeding individuals (N) is N(e) = 4N/[2(1 - α(I)) + (S(k)(2) + 4Q(2)C(2)) (1 + α(I) + 2α(O))], where S(k)(2) is the variance of family size due to sampling without selection, C(2) is the variance of selective advantages among families (the squared coefficient of variation of the expected number of offspring per family), α(I) is the deviation from Hardy-Weinberg proportions, α(O) is the correlation between genes of male and female parents, and Q(2) is the term accounting for the cumulative effect of selection on an inherited trait. This is obtained as Q = 2/[2 - G(1 + r)], where G is the remaining proportion of genetic variance in selected individuals and r is the correlation of the expected selective values of male and female parents. The method is also extended to the general case of different numbers of male and female parents. The predictive value of the formulae is tested under a model of truncation selection with the infinitesimal model of gene effects, where C(2) and G are a function of the selection intensity, the heritability and the intraclass correlation of sibs. Under random mating r = α(I) = -1/(N - 1) and α(O) = 0. Under partial full-sib mating with an average proportion β of full-sib matings per generation, r & β and α(O) & α(I) & β/ (4 - 3β). The prediction equation is compared to other approximations based on the long-term contributions of ancestors to descendants. Finally, based on the approach followed, a system of mating (compensatory mating) is proposed to reduce rates of inbreeding without loss of response in selection programs in which selected individuals from the largest families are mated to those from the smallest families.     

A New Interpretation of the Dynamic Changes of the Potassium Conductance in the Squid Giant Axon

The solutions, n(t), of the differential equation dn/dt = α (1 - n) n (4 - 6n + 4n² - n³) - βn² (4 - 6n + 4n² - n³) in which α and β are instantaneous functions of membrane potential, are shown to fit with good accuracy the time courses of the rise of potassium conductance during depolarizing steps in clamp potential, found experimentally by Hodgkin and Huxley and by Cole and Moore. The equation is derived by analysing the dynamic behaviour of a system consisting of a square array of interacting pores. The possible role of Ca⁺⁺ ions in this system is discussed.     

The Effective Size of a Subdivided Population

This paper derives the long-term effective size, N(e), for a general model of population subdivision, allowing for differential deme fitness, variable emigration and immigration rates, extinction, colonization, and correlations across generations in these processes. We show that various long-term measures of N(e) are equivalent. The effective size of a metapopulation can be expressed in a variety of ways. At a demographic equilibrium, N(e) can be derived from the demography by combining information about the ultimate contribution of each deme to the future genetic make-up of the population and Wright's F(ST)'s. The effective size is given by N(e) = 1/(1 + var ( &))<(1 - f(STi))/N(i)n>, where n is the number of demes, &(i) is the eventual contribution of individuals in deme i to the whole population (scaled such that σ(i) &(i) = n), and < > denotes an average weighted by &(i)(2). This formula is applied to a catastrophic extinction model (where sites are either empty or at carrying capacity) and to a metapopulation model with explicit dynamics, where extinction is caused by demographic stochasticity and by chaos. Contrary to the expectation from the standard island model, the usual effect of population subdivision is to decrease the effective size relative to a panmictic population living on the same resource.     

Precedent Fluctuation of Serum hs-CRP to Albumin Ratios and Mortality Risk of Clinically Stable Hemodialysis Patients

Background

A high sensitivity C-reactive protein to albumin ratio (hs-CRP/Alb) predicts mortality risk in patients with acute kidney injury. However, it varies dynamically. This study was conducted to evaluate whether a variation of this marker was associated with long-term outcome in clinically stable hemodialysis (HD) patients.

Methods

hs-CRP/Alb was checked bimonthly in 284 clinically stable HD outpatients throughout all of 2008. Based on the “slope” of trend equation derived from 5–6 hs-CRP/alb ratios for each patient, the total number of patients was divided into quartiles—Group 1: β≦ −0.13, n = 71; group 2: β>-0.13≦0.003; n = 71, group 3: β>0.003≦0.20; and group 4: β>0.20, n = 71. The observation period was from January 1, 2009 to August 31, 2012.

Results

Group 1+4 showed a worse long-term survival (p = 0.04) and a longer 5-year hospitalization stay than Group 2+3 (38.7±44.4 vs. 16.7±22.4 days, p<0.001). Group 1+4 were associated with older age (OR = 1.03, 95% CI = 1.01–1.05) and a high prevalence of congestive heart failure (OR = 2.02, 95% CI = 1.00–4.11). Standard deviation (SD) of hs-CRP/Alb was associated with male sex (β = 0.17, p = 0.003), higher Davies co-morbidity score (β = 0.16, p = 0.03), and baseline hs-CRP (β = 0.39, p<0.001). Patients with lower baseline and stable trend of hs-CRP/Alb had a better prognosis. By multivariate Cox proportional methods, SD of hs-CRP/alb (HR: 1.05, 95% CI: 1.01–1.08) rather than baseline hs-CRP/Alb was an independent predictive factor for long-term mortality after adjusting for sex and HD vintage.

Conclusion

Clinically stable HD patients with a fluctuating variation of hs-CRP/Alb are characterized by old age, and more co-morbidity, and they tend to have longer subsequent hospitalization stay and higher mortality risk.     

Selection Response in Finite Populations

Formulae were derived to predict genetic response under various selection schemes assuming an infinitesimal model. Account was taken of genetic drift, gametic (linkage) disequilibrium (Bulmer effect), inbreeding depression, common environmental variance, and both initial segregating variance within families (σ(AW0)(2)) and mutational (σ(M)(2)) variance. The cumulative response to selection until generation t(CR(t)) can be approximated as & where N(e) is the effective population size, σ(AW &)(2) = N(e)σ(M)(2) is the genetic variance within families at the steady state (or one-half the genic variance, which is unaffected by selection), and D is the inbreeding depression per unit of inbreeding. R(0) is the selection response at generation 0 assuming preselection so that the linkage disequilibrium effect has stabilized. β is the derivative of the logarithm of the asymptotic response with respect to the logarithm of the within-family genetic variance, i.e., their relative rate of change. R(0) is the major determinant of the short term selection response, but σ(M)(2), N(e) and β are also important for the long term. A selection method of high accuracy using family information gives a small N(e) and will lead to a larger response in the short term and a smaller response in the long term, utilizing mutation less efficiently.     

Maximum Likelihood Estimation of Population Parameters

One of the most important parameters in population genetics is θ = 4N(e)μ where N(e) is the effective population size and μ is the rate of mutation per gene per generation. We study two related problems, using the maximum likelihood method and the theory of coalescence. One problem is the potential improvement of accuracy in estimating the parameter θ over existing methods and the other is the estimation of parameter λ which is the ratio of two θ's. The minimum variances of estimates of the parameter θ are derived under two idealized situations. These minimum variances serve as the lower bounds of the variances of all possible estimates of θ in practice. We then show that Watterson's estimate of θ based on the number of segregating sites is asymptotically an optimal estimate of θ. However, for a finite sample of sequences, substantial improvement over Watterson's estimate is possible when θ is large. The maximum likelihood estimate of λ = θ(1)/θ(2) is obtained and the properties of the estimate are discussed.     

Effective Size of Nonrandom Mating Populations

Nonrandom mating whereby parents are related is expected to cause a reduction in effective population size because their gene frequencies are correlated and this will increase the genetic drift. The published equation for the variance effective size, Ne, which includes the possibility of nonrandom mating, does not take into account such a correlation, however. Further, previous equations to predict effective sizes in populations with partial sib mating are shown to be different, but also incorrect. In this paper, a corrected form of these equations is derived and checked by stochastic simulation. For the case of stable census number, N, and equal progeny distributions for each sex, the equation is [formula: see text], where Sk2 is the variance of family size and alpha is the departure from Hardy-Weinberg proportions. For a Poisson distribution of family size (Sk2 = 2), it reduces to Ne = N/(1 + alpha), as when inbreeding is due to selfing. When nonrandom mating occurs because there is a specified system of partial inbreeding every generation, alpha can be substituted by Wright's FIS statistic, to give the effective size as a function of the proportion of inbred mates.     

The Contribution of Non-Conventional T Cells and NK Cells in the Mycobacterial-Specific IFNγ Response in Bacille Calmette-Guérin (BCG)-Immunized Infants

Background

The Mycobacterium bovis Bacille Calmette-Guérin (BCG) vaccine is given to >120 million infants each year worldwide. Most studies investigating the immune response to BCG have focused on adaptive immunity. However the importance of TCR-gamma/delta (γδ) T cells and NK cells in the mycobacterial-specific immune response is of increasing interest.

Methods

Participants in four age-groups were BCG-immunized. Ten weeks later, in vitro BCG-stimulated blood was analyzed for NK and T cell markers, and intracellular IFNgamma (IFNγ) by flow cytometry. Total functional IFNγ response was calculated using integrated median fluorescence intensity (iMFI).

Results

In infants and children, CD4 and CD4-CD8- (double-negative (DN)) T cells were the main IFNγ-expressing cells representing 43-56% and 27-37% of total CD3+ IFNγ+ T cells respectively. The iMFI was higher in DN T cells compared to CD4 T cells in all age groups, with the greatest differences seen in infants immunized at birth (p=0.002) or 2 months of age (p<0.0001). When NK cells were included in the analysis, they accounted for the majority of total IFNγ-expressing cells and, together with DN Vδ2 γδ T cells, had the highest iMFI in infants immunized at birth or 2 months of age.

Conclusion

In addition to CD4 T cells, NK cells and DN T cells, including Vδ2 γδ T cells, are the key populations producing IFNγ in response to BCG immunization in infants and children. This suggests that innate immunity and unconventional T cells play a greater role in the mycobacterial immune response than previously recognized and should be considered in the design and assessment of novel tuberculosis vaccines.     

Mechanistic Studies of Semicarbazone Triapine Targeting Human Ribonucleotide Reductase in Vitro and in Mammalian Cells: TYROSYL RADICAL QUENCHING NOT INVOLVING REACTIVE OXYGEN SPECIES*

Triapine® (3-aminopyridine-2-carboxaldehyde thiosemicarbazone (3-AP)) is a drug in Phase II trials. One of its established cellular targets is the β₂ subunit of ribonucleotide reductase that requires a diferric-tyrosyl-radical [(Fe^III₂-Y·)(Fe^III₂)] cofactor for de novo DNA biosynthesis. Several mechanisms for 3-AP inhibition of β₂ have been proposed; one involves direct iron chelation from β₂, whereas a second involves Y· destruction by reactive oxygen species formed in situ in the presence of O₂ and reductant by Fe(II)-(3-AP). Inactivation of β₂ can thus arise from cofactor destruction by loss of iron or Y·. In vitro kinetic data on the rates of ⁵⁵Fe and Y· loss from [(⁵⁵Fe^III₂-Y·)(⁵⁵Fe^III₂)]-β₂ under aerobic and anaerobic conditions reveal that Y· loss alone is sufficient for rapid β₂ inactivation. Oxyblot^TM and mass spectrometric analyses of trypsin-digested inhibited β₂, and lack of Y· loss from H₂O₂ and O₂^˙̄ treatment together preclude reactive oxygen species involvement in Y· loss. Three mammalian cell lines treated with 5 μm 3-AP reveal Y· loss and β₂ inactivation within 30-min of 3-AP-exposure, analyzed by whole-cell EPR and lysate assays, respectively. Selective degradation of apo- over [(Fe^III₂-Y·)(Fe^III₂)]-β₂ in lysates, similar iron-content in β₂ immunoprecipitated from 3-AP-treated and untreated [⁵⁵Fe]-prelabeled cells, and prolonged (12 h) stability of the inhibited β₂ are most consistent with Y· loss being the predominant mode of inhibition, with β₂ remaining iron-loaded and stable. A model consistent with in vitro and cell-based biochemical studies is presented in which Fe(II)-(3-AP), which can be cycled with reductant, directly reduces Y· of the [(Fe^III₂-Y·)(Fe^III₂)] cofactor of β₂.     

Factors Controlling the Stable Nitrogen Isotopic Composition (δ15N) of Lipids in Marine Animals

Lipid extraction of biomass prior to stable isotope analysis is known to cause variable changes in the stable nitrogen isotopic composition (δ¹⁵N) of residual biomass. However, the underlying factors causing these changes are not yet clear. Here we address this issue by comparing the δ¹⁵N of bulk and residual biomass of several marine animal tissues (fish, crab, cockle, oyster, and polychaete), as well as the δ¹⁵N of the extracted lipids. As observed previously, lipid extraction led to a variable offset in δ¹⁵N of biomass (differences ranging from -2.3 to +1.8 ‰). Importantly, the total lipid extract (TLE) was highly depleted in ¹⁵N compared to bulk biomass, and also highly variable (differences ranging from -14 to +0.7 ‰). The TLE consisted mainly of phosphatidylcholines, a group of lipids with one nitrogen atom in the headgroup. To elucidate the cause for the ¹⁵N-depletion in the TLE, the δ¹⁵N of amino acids was determined, including serine because it is one of the main sources of nitrogen to N-containing lipids. Serine δ¹⁵N values differed by -7 to +2 ‰ from bulk biomass δ¹⁵N, and correlated well with the ¹⁵N depletion in TLEs. On average, serine was less depleted (-3‰) than the TLE (-7 ‰), possibly due to fractionation during biosynthesis of N-containing headgroups, or that other nitrogen-containing compounds, such as urea and choline, or recycled nitrogen contribute to the nitrogen isotopic composition of the TLE. The depletion in ¹⁵N of the TLE relative to biomass increased with the trophic level of the organisms.     

Requirement for Calcium in Adhesion of a Fouling Diatom to Glass

An assay for the adhesion of Amphora coffeaeformis (Bacillariophyceae) to glass is described. Ca²⁺ and Sr²⁺ promoted adhesion, but Ba²⁺ and Mn²⁺ were ineffective. Adhesion was inhibited by cycloheximide, carbonyl cyanide 3-chlorophenyl hydrazone, and α-isopropyl-α-[(N-methyl-N-homoveratryl)-γ-amino- propyl]-3,4,5-trimethyoxy phenyl acetonitrile (D-600), but not by 3-(3,4-dichlorophenyl)-1, 1-dimethyl urea or darkness. The assay has potential use in comparing chemical antifoulants.     

Structure and Location of the Regulatory β Subunits in the (αβγδ)4 Phosphorylase Kinase Complex

Phosphorylase kinase (PhK) is a hexadecameric (αβγδ)₄ complex that regulates glycogenolysis in skeletal muscle. Activity of the catalytic γ subunit is regulated by allosteric activators targeting the regulatory α, β, and δ subunits. Three-dimensional EM reconstructions of PhK show it to be two large (αβγδ)₂ lobes joined with D₂ symmetry through interconnecting bridges. The subunit composition of these bridges was unknown, although indirect evidence suggested the β subunits may be involved in their formation. We have used biochemical, biophysical, and computational approaches to not only address the quaternary structure of the β subunits within the PhK complex, i.e. whether they compose the bridges, but also their secondary and tertiary structures. The secondary structure of β was determined to be predominantly helical by comparing the CD spectrum of an αγδ subcomplex with that of the native (αβγδ)₄ complex. An atomic model displaying tertiary structure for the entire β subunit was constructed using chemical cross-linking, MS, threading, and ab initio approaches. Nearly all this model is covered by two templates corresponding to glycosyl hydrolase 15 family members and the A subunit of protein phosphatase 2A. Regarding the quaternary structure of the β subunits, they were directly determined to compose the four interconnecting bridges in the (αβγδ)₄ kinase core, because a β₄ subcomplex was observed through both chemical cross-linking and top-down MS of PhK. The predicted model of the β subunit was docked within the bridges of a cryoelectron microscopic density envelope of PhK utilizing known surface features of the subunit.     

Role of Sulfate Reduction Versus Methanogenesis in Terminal Carbon Flow in Polluted Intertidal Sediment of Waimea Inlet, Nelson, New Zealand

An investigation of the terminal anaerobic processes occurring in polluted intertidal sediments indicated that terminal carbon flow was mainly mediated by sulfate-reducing organisms in sediments with high sulfate concentrations (>10 mM in the interstitial water) exposed to low loadings of nutrient (equivalent to <10² kg of N · day⁻¹) and biochemical oxygen demand (<0.7 × 10³ kg · day⁻¹) in effluents from different pollution sources. However, in sediments exposed to high loadings of nutrient (>10² kg of N · day⁻¹) and biochemical oxygen demand (>0.7 × 10³ kg · day⁻¹), methanogenesis was the major process in the mediation of terminal carbon flow, and sulfate concentrations were low (≤2 mM). The respiratory index [¹⁴CO₂/(¹⁴CO₂ + ¹⁴CH₄)] for [2-¹⁴C]acetate catabolism, a measure of terminal carbon flow, was ≥0.96 for sediment with high sulfate, but in sediments with sulfate as little as 10 μM in the interstitial water, respiratory index values of ≤0.22 were obtained. In the latter sediment, methane production rates as high as 3 μmol · g⁻¹ (dry weight) · h⁻¹ were obtained, and there was a potential for active sulfate reduction.     

Dynamic phospholipid interaction of β2e subunit regulates the gating of voltage-gated Ca2+ channels

High voltage-activated Ca²⁺ (Ca_V) channels are protein complexes containing pore-forming α1 and auxiliary β and α2δ subunits. The subcellular localization and membrane interactions of the β subunits play a crucial role in regulating Ca_V channel inactivation and its lipid sensitivity. Here, we investigated the effects of membrane phosphoinositide (PI) turnover on Ca_V2.2 channel function. The β2 isoform β2e associates with the membrane through electrostatic and hydrophobic interactions. Using chimeric β subunits and liposome-binding assays, we determined that interaction between the N-terminal 23 amino acids of β2e and anionic phospholipids was sufficient for β2e membrane targeting. Binding of the β2e subunit N terminus to liposomes was significantly increased by inclusion of 1% phosphatidylinositol 4,5-bisphosphate (PIP₂) in the liposomes, suggesting that, in addition to phosphatidylserine, PIs are responsible for β2e targeting to the plasma membrane. Membrane binding of the β2e subunit slowed Ca_V2.2 current inactivation. When membrane phosphatidylinositol 4-phosphate and PIP₂ were depleted by rapamycin-induced translocation of pseudojanin to the membrane, however, channel opening was decreased and fast inactivation of Ca_V2.2(β2e) currents was enhanced. Activation of the M₁ muscarinic receptor elicited transient and reversible translocation of β2e subunits from membrane to cytosol, but not that of β2a or β3, resulting in fast inactivation of Ca_V2.2 channels with β2e. These results suggest that membrane targeting of the β2e subunit, which is mediated by nonspecific electrostatic insertion, is dynamically regulated by receptor stimulation, and that the reversible association of β2e with membrane PIs results in functional changes in Ca_V channel gating. The phospholipid–protein interaction observed here provides structural insight into mechanisms of membrane–protein association and the role of phospholipids in ion channel regulation.     

The N-terminal region of DNA polymerase delta catalytic subunit is necessary for holoenzyme function

Genetic and biochemical studies have shown that DNA polymerase δ (Polδ) is the major replicative Pol in the eukaryotic cell. Its functional form is the holoenzyme composed of Polδ, proliferating cell nuclear antigen (PCNA) and replication factor C (RF-C). In this paper, we describe an N-terminal truncated form of DNA polymerase δ (ΔN Polδ) from calf thymus. The ΔN Polδ was stimulated as the full-length Polδ by PCNA in a RF-C-independent Polδ assay. However, when tested for holoenzyme function in a RF-C-dependent Polδ assay in the presence of RF-C, ATP and replication protein A (RP-A), the ΔN Polδ behaved differently. First, the ΔN Polδ lacked holoenzyme functions to a great extent. Second, product size analysis and kinetic experiments showed that the holoenzyme containing ΔN Polδ was much less efficient and synthesized DNA at a much slower rate than the holoenzyme containing full-length Polδ. The present study provides the first evidence that the N-terminal part of the large subunit of Polδ is involved in holoenzyme function.     

Alanine Scan of α-Conotoxin RegIIA Reveals a Selective α3β4 Nicotinic Acetylcholine Receptor Antagonist

Activation of the α3β4 nicotinic acetylcholine receptor (nAChR) subtype has recently been implicated in the pathophysiology of various conditions, including development and progression of lung cancer and in nicotine addiction. As selective α3β4 nAChR antagonists, α-conotoxins are valuable tools to evaluate the functional roles of this receptor subtype. We previously reported the discovery of a new α4/7-conotoxin, RegIIA. RegIIA was isolated from Conus regius and inhibits acetylcholine (ACh)-evoked currents mediated by α3β4, α3β2, and α7 nAChR subtypes. The current study used alanine scanning mutagenesis to understand the selectivity profile of RegIIA at the α3β4 nAChR subtype. [N11A] and [N12A] RegIIA analogs exhibited 3-fold more selectivity for the α3β4 than the α3β2 nAChR subtype. We also report synthesis of [N11A,N12A]RegIIA, a selective α3β4 nAChR antagonist (IC₅₀ of 370 nm) that could potentially be used in the treatment of lung cancer and nicotine addiction. Molecular dynamics simulations of RegIIA and [N11A,N12A]RegIIA bound to α3β4 and α3β2 suggest that destabilization of toxin contacts with residues at the principal and complementary faces of α3β2 (α3-Tyr⁹², Ser¹⁴⁹, Tyr¹⁸⁹, Cys¹⁹², and Tyr¹⁹⁶; β2-Trp⁵⁷, Arg⁸¹, and Phe¹¹⁹) may form the molecular basis for the selectivity shift.     

Microbial Transformation of β-Ionone and β-Methylionone

Aspergillus niger JTS 191 was selected from many microorganisms tested as capable of converting ionones to other compounds having aromas. The individual transformation products from β-ionone were isolated and identified by comparison with synthetically derived compounds. The major products were (R)-4-hydroxy-β-ionone and (S)-2-hydroxy-β-ionone. 2-Oxo-, 4-oxo-, 3,4-dehydro-, 2,3-dehydro-4-oxo-, 3,4-dehydro-2-oxo-, (S)-2-acetoxy-, (R)-4-acetoxy-, and 5,6-epoxy-β-ionone and 4-(2,3,6-trimethylphenyl)-but-3-en-2-one were also identified. Analogous transformation products of β-methylionone also were identified. Based on gas-liquid chromatographic analysis during the fermentation, we propose two main oxidative pathways of β-ionone. The results of this study suggest that these transformations of β-ionones may be useful as tobacco-flavoring compounds.     

Calcium is Necessary for Motility in the Diatom Amphora coffeaeformis

The marine diatom Amphora coffeaeformis required Ca²⁺ and bicarbonate for motility. Movement was inhibited by the Ca²⁺-blocking agents ruthenium red and α-isopropyl-α-[(N-methyl-N-homoveratryl)-α- aminopropyl]-3,4,5-trimethoxy phenyl acetonitrile and the metabolic energy uncoupler, carbonyl cyanide 3-chlorophenylhydrazone. 3-(3′,4-Dichlorophenyl)-1,1-Dimethyl urea was without effect on cells at a concentration that prevented O₂ production in the light. Although Sr²⁺ could replace Ca²⁺ in the attachment of cells to glass, it did not substitute for Ca²⁺ in motility.     

Determination of Fluorescence Polarization of Double Chromophore Complexes

Polarized fluorescence of rigid double-chromophore complexes with intracomplex energy exchange between chromophores was analyzed, and the formula for the degree of polarization derived for the case of steady-state excitation: P = (3 cos²θ - 1 + 2A)/(3 + cos²θ + 4A). In this formula θ is the angle between the transition dipole moments of chromophores in complexes, and A is the parameter dependent on the spectroscopic features of chromophores and energy migration rates. The case of excitation by a δ-pulse was also analyzed, and a formula for fluorescence polarization kinetics was derived.As an example of the application of the derived formulae, the polarized fluorescence spectra and their picosecond kinetics were calculated for the β-subunits of the blue-green algae Agmenellum quadruplicatum. The results obtained were compared with experimental measurements of Mimuro et al. (1986, Biochim. Biophys. Acta848, 155-166) and found to match these data well.     

Cerebrospinal Fluid (CSF) CD8+ T-Cells That Express Interferon-Gamma Contribute to HIV Associated Neurocognitive Disorders (HAND)

BackgroundHIV associated neurocognitive disorders (HAND) continue to affect cognition and everyday functioning despite anti-retroviral treatment (ART). Previous studies focused on mechanisms related to monocyte/macrophage mediated inflammation. However, in the ART era, there is increasing evidence for the involvement of CD8+ T-cells in CNS pathogenesis.MethodsTo investigate the relationship between T-cell responses and neurocognitive impairment (NCI), cerebrospinal fluid (CSF) and peripheral blood CD4+ and CD8+ T-cell intracellular cytokine (IFNγ, IL-2, TNFα) and lytic marker (CD107a) expression were assessed in HIV infected subjects who underwent comprehensive neurocognitive (NC) evaluation and either initiated or changed ART.ResultsData were collected from 31 participants at 70 visits. The frequency of cytokine expressing T-cells in CSF was significantly higher than in peripheral blood for CD4+T-cells: TNFα, IL-2, IFNγ and CD8+T-cells: IL-2 and IFNγ. Analysis of T-cell activity and NCI as a function of CSF HIV RNA levels suggested a general association between NCI, high CSF CD8+ (but not CD4+T-cell) cytokine expression and CSF HIV RNA <10³ copies/ml (p<0.0001). Specifically, CSF CD8+ T-cell IFNγ expression correlated with severity of NCI (r = 0.57, p = 0.004). Multivariable analyses indicated that CSF CD8+T-cell IFNγ and myeloid activation (CD163) contributed equally and independently to cognitive status and a composite variable produced the strongest correlation with NCI (r = 0.83, p = 0.0001). In contrast, CD8+ cytolytic activity (CD107a expression) was negatively correlated with NCI (p = 0.05) but was dependent on CD4 levels >400/μl and low CSF HIV RNA levels (<10³ copies/ml). In our longitudinal analysis of 16 subjects, higher CSF CD8+IFNγ expression at baseline predicted NC decline at follow-up (p = 0.02). Severity of NCI at follow-up correlated with level of residual HIV RNA in CSF.ConclusionsPresence of IFNγ expressing CD8+ T-cells, absence of cytolytic CD8+ T-cells, high myeloid activation, and failure of ART to suppress HIV replication in CSF contribute to increased risk of HAND.