Footprints of ancient‐balanced polymorphisms in genetic variation data from closely related species

When long‐lasting, balancing selection can lead to “trans‐species” polymorphisms that are shared by two or more species identical by descent. In such cases, the gene genealogy at the selected site clusters by allele instead of by species, and nearby neutral sites also have unusual genealogies because of linkage. While this scenario is expected to leave discernible footprints in genetic variation data, the specific patterns remain poorly characterized. Motivated by recent findings in primates, we focus on the case of a biallelic polymorphism under ancient balancing selection and derive approximations for summaries of the polymorphism data from two species. Specifically, we characterize the length of the segment that carries most of the footprints, the expected number of shared neutral single nucleotide polymorphisms (SNPs), and the patterns of allelic associations among them. We confirm the accuracy of our approximations by coalescent simulations. We further show that for humans and chimpanzees—more generally, for pairs of species with low genetic diversity levels—these patterns are highly unlikely to be generated by neutral recurrent mutations. We discuss the implications for the design and interpretation of genome scans for ancient balanced polymorphisms in primates and other taxa.     

High-Specificity Targeted Functional Profiling in Microbial Communities with ShortBRED

Profiling microbial community function from metagenomic sequencing data remains a computationally challenging problem. Mapping millions of DNA reads from such samples to reference protein databases requires long run-times, and short read lengths can result in spurious hits to unrelated proteins (loss of specificity). We developed ShortBRED (Short, Better Representative Extract Dataset) to address these challenges, facilitating fast, accurate functional profiling of metagenomic samples. ShortBRED consists of two components: (i) a method that reduces reference proteins of interest to short, highly representative amino acid sequences (“markers”) and (ii) a search step that maps reads to these markers to quantify the relative abundance of their associated proteins. After evaluating ShortBRED on synthetic data, we applied it to profile antibiotic resistance protein families in the gut microbiomes of individuals from the United States, China, Malawi, and Venezuela. Our results support antibiotic resistance as a core function in the human gut microbiome, with tetracycline-resistant ribosomal protection proteins and Class A beta-lactamases being the most widely distributed resistance mechanisms worldwide. ShortBRED markers are applicable to other homology-based search tasks, which we demonstrate here by identifying phylogenetic signatures of antibiotic resistance across more than 3,000 microbial isolate genomes. ShortBRED can be applied to profile a wide variety of protein families of interest; the software, source code, and documentation are available for download at http://huttenhower.sph.harvard.edu/shortbred     

Vertical bars on male Xiphophorus multilineatus: a signal that deters rival males and attracts females

We examined the function of the vertical bar pattern on maleswordtails (Xiphophorus multilinneatus) as a signal in bothmale-male competition and female choice. This pattern had previouslybeen described as an aggressive signal because males intensifiedthe bars during male-male encounters in the laboratory. Ourfield observations supported this observation and also showedthat bars intensified when males courted females. The intensityof bars was correlated with access to females in the field.Within the size range of males that have bars, however, neitherbar number nor male size appeared to influence access to females.We used freeze-branding to remove the bars from males in thelaboratory so that we could control for characters correlatedwith bar intensity, and tested males and females separatelyso that we could separate the influence of these two componentsof sexual selection. We compared the responses of males andfemales to males that had their bars removed and control malesfreeze-branded between the bars. Test males responded more aggressivelyto males without bars as compared to control males. In addition,females showed a preference for control males over males thathad their bars removed. These results suggest that the barsmay function as a signal that deters rival males and attractsfemales.     

Biochemical Screening of Five Protein Kinases from Plasmodium falciparum against 14,000 Cell-Active Compounds

In 2010 the identities of thousands of anti-Plasmodium compounds were released publicly to facilitate malaria drug development. Understanding these compounds’ mechanisms of action—i.e., the specific molecular targets by which they kill the parasite—would further facilitate the drug development process. Given that kinases are promising anti-malaria targets, we screened ~14,000 cell-active compounds for activity against five different protein kinases. Collections of cell-active compounds from GlaxoSmithKline (the ~13,000-compound Tres Cantos Antimalarial Set, or TCAMS), St. Jude Children’s Research Hospital (260 compounds), and the Medicines for Malaria Venture (the 400-compound Malaria Box) were screened in biochemical assays of Plasmodium falciparum calcium-dependent protein kinases 1 and 4 (CDPK1 and CDPK4), mitogen-associated protein kinase 2 (MAPK2/MAP2), protein kinase 6 (PK6), and protein kinase 7 (PK7). Novel potent inhibitors (IC₅₀ < 1 μM) were discovered for three of the kinases: CDPK1, CDPK4, and PK6. The PK6 inhibitors are the most potent yet discovered for this enzyme and deserve further scrutiny. Additionally, kinome-wide competition assays revealed a compound that inhibits CDPK4 with few effects on ~150 human kinases, and several related compounds that inhibit CDPK1 and CDPK4 yet have limited cytotoxicity to human (HepG2) cells. Our data suggest that inhibiting multiple Plasmodium kinase targets without harming human cells is challenging but feasible.     

Karyometric image analysis for intraepithelial and invasive cervical lesions

OBJECTIVE: To derive an objective, numeric measure for the progression of intraepithelial and invasive squamous cell cervical lesions. STUDY DESIGN: Thin-layer cervical cytology preparations from colposcopically confirmed normal cervix, low grade squamous intraepithelial lesions, high grade squamous intraepithelial lesions and carcinoma were identified from a cross-sectional study. Fifty-nine cases representing 4 diagnostic categories were selected, and 2,375 nuclei from epithelial cells representative of the diagnostic category were randomly selected for imaging and measurement from these cases. Additionally, 1,378 visually normal appearing intermediate cells from low and high grade squamous intraepithelial lesions, as well as from carcinoma cases, were identified for analysis. The nuclei were quantitatively characterized, and discriminant analyses were performed to derive a progression curve from normal cytology to carcinoma. RESULTS: The lesion signatures show a clear increase in nuclear abnormality with increasing progression. A progression curve was derived based on mean discriminant function scores for each diagnostic category and on the mean nuclear abnormality values for the nuclei in each category, as expressed by their deviation in feature values from normal reference nuclei. CONCLUSION: A numeric assessment of lesion progression for cervical precancerous and cancerous lesions based on karyometric measurements is possible and may provide an objective, precise characterization of each lesion as well as a basis for improved performance in automated cytology-based cervical cancer screening.     

NMR Structural Studies of Domain 1 of Receptor-associated Protein

The 39 kDa receptor-associated protein (RAP) is an endoplasmic reticulum resident protein that binds tightly to the low-density lipoprotein receptor-related protein (LRP) as well as to other members of the low-density lipoprotein receptor superfamily. The association of RAP with LRP prevents this receptor from interacting with ligands. RAP is a three-domain protein that contains two independent LRP binding sites; one located within domains 1 and 2, and one located within domain 3. As the first step toward defining the structure of the full-length protein and understanding the interaction between RAP and this family of receptors, we have determined the 3D structure of domain 1 using constraints derived from heteronuclear multi-dimensional NMR spectra, including NOEs, dihedral angles, J-couplings and chemical shifts, as well as two sets of non-correlated residual dipolar couplings measured from the protein solutions in anisotropic media of Pf1 and 6% polyacrylamide gel. The backbone C(alpha) rmsd between the current structure and a homo-nuclear NOE-based structure is about 2 A. The large rmsd mainly reflects the significant differences in helical orientation and in the structural details of the long helix (helix 2) between the two structures.     

Tackling the plant proteome: practical approaches, hurdles and experimental tools

The study of complex biological questions through comparative proteomics is becoming increasingly attractive to plant biologists as the rapidly expanding plant genomic and expressed sequence tag databases provide improved opportunities for protein identification. This review focuses on practical issues associated with comparative proteomic analysis, including the challenges of effective protein extraction and separation from plant tissues, the pros and cons of two-dimensional gel-based analysis and the problems of identifying proteins from species that are not recognized models for functional genomic studies. Specific points are illustrated using data from an ongoing study of the tomato and pepper fruit proteomes.