Toward a unified theory of sexual mixing and pair formation.

Sexually transmitted diseases such as gonorrhea, syphilis, herpes, and AIDS are driven and maintained in populations by epidemiological and sociological factors that are not completely understood. One such factor is the way in which people mix sexually. In this paper, we outline a unified approach to modeling sexual mixing structures, where such structures are defined in terms of a set of axioms for a finite number of distinct groups of people. Theorems for homosexual, heterosexual, and arbitrary group mixing are presented, leading to a representation of all mixing structures defined by the axioms. The representation and its parameters are interpreted in terms of intergroup affinities for sexual mixing. The use of the approach in sexually transmitted disease modeling is discussed.     

The dynamics of population models with distributed maturation periods

An integro-differential equation for the dynamics of a subpopulation of adults in a closed system where only the adults compete and where there is a distribution of maturation periods is described. We show how the careful choice of a general weighting function based on the gamma distribution with a shift in origin enables us to characterize adequately some observed maturation-period distributions, and also makes local stability and numerical analyses straightforward. Using these results we examine the progression in the behavior of the distributed-delay model as the distribution is narrowed toward the limit of a discrete delay. We conclude that while local stability properties approach those of the limiting equation very rapidly, the persistent fluctuation behavior converges more slowly, with the dominant period and maximum amplitude being least affected by the details of the distribution, and the fine structure of solutions being most sensitive. Finally, we examine the consequences for population modeling, and using several examples of insect populations, conclude that although quite often a full maturation-period distribution should be incorporated in a given model, in many cases a discrete-delay approximation will suffice.     

Arrangement of Integrated Avian Sarcoma Virus DNA Sequences Within the Cellular Genomes of Transformed and Revertant Mammalian Cells

We have examined the arrangement of integrated avian sarcoma virus (ASV) DNA sequences in several different avian sarcoma virus transformed mammalian cell lines, in independently isolated clones of avian sarcoma virus transformed rat liver cells, and in morphologically normal revertants of avian sarcoma virus transformed rat embryo cells. By using restriction endonuclease digestion, agarose gel electrophoresis, Southern blotting, and hybridization with labeled avian sarcoma virus complementary DNA probes, we have compared the restriction enzyme cleavage maps of integrated viral DNA and adjacent cellular DNA sequences in four different mouse and rat cell lines transformed with either Bratislava 77 or Schmidt-Ruppin strains of avian sarcoma virus. The results of these experiments indicated that the integrated viral DNA resided at a different site within the host cell genome in each transformed cell line. A similar analysis of several independently derived clones of Schmidt-Ruppin transformed rat liver cells also revealed that each clone contained a unique cellular site for the integration of proviral DNA. Examination of several morphologically normal revertants and spontaneous retransformants of Schmidt-Ruppin transformed rat embryo cells revealed that the internal arrangement and cellular integration site of viral DNA sequences was identical with that of the transformed parent cell line. The loss of the transformed phenotype in these revertant cell lines, therefore, does not appear to be the result of rearrangement or deletions either within the viral genome or in adjacent cellular DNA sequences. The data presented support a model for ASV proviral DNA integration in which recombination can occur at multiple sites within the mammalian cell genome. The integration and maintenance of at least one complete copy of the viral genome appear to be required for continuous expression of the transformed phenotype in mammalian cells.     

Interleukin-15 rescues tolerant CD8+ T cells for use in adoptive immunotherapy of established tumors

CD8+ T cells can mediate eradication of established tumors, and strategies to amplify tumor-reactive T-cell numbers by immunization or ex vivo expansion followed by adoptive transfer are currently being explored in individuals with cancer. Generating effective CD8+ T cell-mediated responses to tumors is often impeded by T-cell tolerance to relevant tumor antigens, as most of these antigens are also expressed in normal tissues. We examined whether such tolerant T cells could be rescued and functionally restored for use in therapy of established tumors. We used a transgenic T-cell receptor (TCR) mouse model in which peripheral CD8+ T cells specific for a candidate tumor antigen also expressed in liver are tolerant, failing to proliferate or secrete interleukin (IL)-2 in response to antigen. Molecular and cellular analysis showed that these tolerant T cells expressed the IL-15 receptor alpha chain, and could be induced to proliferate in vitro in response to exogenous IL-15. Such proliferation abrogated tolerance and the rescued cells became effective in treating leukemia. Therefore, high-affinity CD8+ T cells are not necessarily deleted by encounter with self-antigen in the periphery, and can potentially be rescued and expanded for use in tumor immunotherapy.     

Exome sequencing in multiplex families with left-sided cardiac defects has high yield for disease gene discovery

Congenital heart disease (CHD) is a common group of birth defects with a strong genetic contribution to their etiology, but historically the diagnostic yield from exome studies of isolated CHD has been low. Pleiotropy, variable expressivity, and the difficulty of accurately phenotyping newborns contribute to this problem. We hypothesized that performing exome sequencing on selected individuals in families with multiple members affected by left-sided CHD, then filtering variants by population frequency, in silico predictive algorithms, and phenotypic annotations from publicly available databases would increase this yield and generate a list of candidate disease-causing variants that would show a high validation rate. In eight of the nineteen families in our study (42%), we established a well-known gene/phenotype link for a candidate variant or performed confirmation of a candidate variant’s effect on protein function, including variants in genes not previously described or firmly established as disease genes in the body of CHD literature: BMP10, CASZ1, ROCK1 and SMYD1. Two plausible variants in different genes were found to segregate in the same family in two instances suggesting oligogenic inheritance. These results highlight the need for functional validation and demonstrate that in the era of next-generation sequencing, multiplex families with isolated CHD can still bring high yield to the discovery of novel disease genes.     

Intratumor heterogeneity of biomarker expression in breast carcinomas

Small biopsy samples are used increasingly to assess the biomarker expression for prognostic information and for monitoring therapeutic responses prior to and during neoadjuvant therapy. The issue of intratumor heterogeneity of expression of biomarkers, however, has raised questions about the validity of the assessment of biomarker expression based on limited tissue samples. We examined immunohistochemically the expression of HER-2neu (p185^erbB-2), epidermal growth factor receptor (EGFR), Bcl-2, p53, and proliferating cell nuclear antigen (PCNA) in 30 breast carcinomas using archived, paraffin embedded tissue and determined the extent of intratumor heterogeneity. Each section was divided into four randomly oriented discrete regions, each containing a portion of the infiltrating carcinoma. For each tumor, the entire lesion and four regions were analyzed for the expression of these markers. Scores of both membrane and cytoplasmic staining of HER-2neu and EGFR, scores of cytoplasmic staining of Bcl-2, and scores of nuclear staining of both p53 and PCNA were recorded. The intensity of staining and the proportion of immunostained cells were determined. A semiquantitative immunoscore was calculated by determining the sum of the products of the intensity and corresponding proportion of stained tumor cells. We analyzed both invasive (IDC) and in situ (DCIS) carcinomas. The Wilcoxon signed-rank test was used for paired comparisons between overall and regional immunoscores and between overall and regional percentages of stained cells. Spearman's correlation coefficients were used to assess the level of agreement of overall biomarker expression with each of the regions. Generalized linear models were used to assess overall and pair-wise differences in the absolute values of percent changes between overall and regional expression of biomarkers. For IDCs, there were no statistically significant differences in the expression of the biomarkers in terms of either the percentage of cells staining or the immunoscores when comparing the entire tumor with each region except for the lower EGFR expression of arbitrarily selected region 1 and lower p53 expression of region 1 compared to that of the entire tumor section. For DCIS, there were no statistically significant differences in the expression of the biomarkers between the entire tumor and each region except in PCNA of region 2 compared to that of entire tumor section. Positive correlation of immunoscores was observed between the entire tumor and each region as well as across all four regions for IDC. Similar observations were noted with DCIS except for HER-2neu and PCNA. No statistically significant differences were observed in the absolute values of percent changes of biomarker expression between overall and the four regions for both DCIS and IDC. Therefore, no significant intratumor heterogeneity in the expression of HER-2neu, Bcl-2, and PCNA was observed in IDC. Minor regional variations were observed for EGFR and p53 in IDC. Similarly, no significant regional variation in the expression of markers was observed in DCIS except for PCNA.     

Thymineless death is inhibited by CsrA in Escherichia coli lacking the SOS response

Thymineless death (TLD) is the rapid loss of colony-forming ability in bacterial, yeast and human cells starved for thymine, and is the mechanism of action of common chemotherapeutic drugs. In Escherichia coli, significant loss of viability during TLD requires the SOS replication-stress/DNA-damage response, specifically its role in inducing the inhibitor of cell division, SulA. An independent RecQ- and RecJ-dependent TLD pathway accounts for a similarly large additional component of TLD, and a third SOS- and RecQ/J-independent TLD pathway has also been observed. Although two groups have implicated the SOS-response in TLD, an SOS-deficient mutant strain from an earlier study was found to be sensitive to thymine deprivation. We performed whole-genome resequencing on that SOS-deficient strain and find that, compared with the SOS-proficient control strain, it contains five mutations in addition to the SOS-blocking lexA(Ind⁻) mutation. One of the additional mutations, csrA, confers TLD sensitivity specifically in SOS-defective strains. We find that CsrA, a carbon storage regulator, reduces TLD in SOS- or SulA-defective cells, and that the increased TLD that occurs in csrA⁻ SOS-defective cells is dependent on RecQ. We consider a hypothesis in which the modulation of nucleotide pools by CsrA might inhibit TLD specifically in SOS-deficient (SulA-deficient) cells.     

Estimation of transformation parameters for microarray data

MOTIVATION AND RESULTS: Durbin et al. (2002), Huber et al. (2002) and Munson (2001) independently introduced a family of transformations (the generalized-log family) which stabilizes the variance of microarray data up to the first order. We introduce a method for estimating the transformation parameter in tandem with a linear model based on the procedure outlined in Box and Cox (1964). We also discuss means of finding transformations within the generalized-log family which are optimal under other criteria, such as minimum residual skewness and minimum mean-variance dependency. AVAILABILITY: R and Matlab code and test data are available from the authors on request.     

The cortisol awakening response (CAR) in male children with autism spectrum disorder

Our ability to adapt to change is fundamental. The cortisol awakening response (CAR) is a sharp rise in cortisol 30 min after waking to help prepare an individual for ensuing stress. Children with autism spectrum disorder (ASD) often have difficulty adapting to change. Exploration of the CAR is warranted; yet, the few studies investigating it are inconclusive. The CAR was investigated in 94 pre-pubertal male children 8-to-12 years of age with ASD (46) and typical development (TD, 48). Salivary samples were collected over three diurnal cycles involving two morning samples: M1: Immediately upon Waking and M2: 30-min Post Waking (M2 − M1 = CAR). The magnitude of the CAR was measured by independent two sample t-tests, variability was measured using Levene's Test, the sequence of the CAR was analyzed by a linear mixed-effects model and proportion of children exhibiting a CAR by chi-square test of independence. There were no significant differences on the CAR between the groups based on magnitude (t(92) = − 0.14, p = 0.89, d = 0.04), variability (F(45,47) = 1.11, p = 0.72, η² = 0.11) or the sequence over three days (F(2,88) = 0.26, p = 0.77, η² = 0.01). No significant differences were shown in the proportion of children exhibiting a CAR across the groups based on child (χ²(1) = 0.02, p = 0.89) or adult criterion (χ²(1) = 1.82, p = 0.18). Despite group differences in the regulation and responsivity of cortisol, the CAR is indistinguishable between children with and without ASD. Inconsistencies across studies may be due to age, criterion used, and diagnostic distinctions.     

Adenoviral expression of 15-lipoxygenase-1 in rabbit aortic endothelium: role in arachidonic acid-induced relaxation

Endothelium-dependent vasorelaxation of the rabbit aorta is mediated by either nitric oxide (NO) or arachidonic acid (AA) metabolites from cyclooxygenase (COX) and 15-lipoxygenase (15-LO) pathways. 15-LO-1 metabolites of AA, 11,12,15-trihydroxyeicosatrienoic acid (THETA), and 15-hydroxy-11,12-epoxyeicosatrienoic acid (HEETA) cause concentration-dependent relaxation. We tested the hypothesis that in the 15-LO pathway of AA metabolism, 15-LO-1 is sufficient and is the rate-limiting step in inducing relaxations in rabbit aorta. Aorta and rabbit aortic endothelial cells were treated with adenoviruses containing human 15-LO-1 cDNA (Ad-15-LO-1) or beta-galactosidase (Ad-beta-Gal). Ad-15-LO-1-transduction increased the expression of a 75-kDa protein corresponding to 15-LO-1, detected by immunoblotting with an anti-human15-LO-1 antibody, and increased the production of HEETA and THETA from [(14)C]AA. Immunohistochemical studies on Ad-15-LO-1-transduced rabbit aorta showed the presence of 15-LO-1 in endothelial cells. Ad-15-LO-1-treated aortic rings showed enhanced relaxation to AA (max 31.7 +/- 3.2%) compared with Ad-beta-Gal-treated (max 12.7 +/- 3.2%) or control nontreated rings (max 13.1 +/- 1.6%) (P < 0.01). The relaxations in Ad-15-LO-1-treated aorta were blocked by the 15-LO inhibitor cinnamyl-3,4-dihydroxy-a-cyanocinnamate. Overexpression of 15-LO-1 in the rabbit aortic endothelium is sufficient to increase the production of the vasodilatory HEETA and THETA and enhance the relaxations to AA. This confirms the role of HEETA and THETA as endothelium-derived relaxing factors.