Sources of selection bias in evaluating social programs: An interpretation of conventional measures and evidence on the effectiveness of matching as a program evaluation method

This paper decomposes the conventional measure of selection bias in observational studies into three components. The first two components are due to differences in the distributions of characteristics between participant and nonparticipant (comparison) group members: the first arises from differences in the supports, and the second from differences in densities over the region of common support. The third component arises from selection bias precisely defined. Using data from a recent social experiment, we find that the component due to selection bias, precisely defined, is smaller than the first two components. However, selection bias still represents a substantial fraction of the experimental impact estimate. The empirical performance of matching methods of program evaluation is also examined. We find that matching based on the propensity score eliminates some but not all of the measured selection bias, with the remaining bias still a substantial fraction of the estimated impact. We find that the support of the distribution of propensity scores for the comparison group is typically only a small portion of the support for the participant group. For values outside the common support, it is impossible to reliably estimate the effect of program participation using matching methods. If the impact of participation depends on the propensity score, as we find in our data, the failure of the common support condition severely limits matching compared with random assignment as an evaluation estimator.

Rapid evolution in plant chitinases: Molecular targets of selection in plant-pathogen coevolution

Many pathogen recognition genes, such as plant R-genes, undergo rapid adaptive evolution, providing evidence that these genes play a critical role in plant-pathogen coevolution. Surprisingly, whether rapid adaptive evolution also occurs in genes encoding other kinds of plant defense proteins is unknown. Unlike recognition proteins, plant chitinases attack pathogens directly, conferring disease resistance by degrading chitin, a component of fungal cell walls. Here, we show that nonsynonymous substitution rates in plant class I chitinase often exceed synonymous rates in the plant genus Arabis (Cruciferae) and in other dicots, indicating a succession of adaptively driven amino acid replacements. We identify individual residues that are likely subject to positive selection by using codon substitution models and determine the location of these residues on the three-dimensional structure of class I chitinase. In contrast to primate lysozymes and plant class III chitinases, structural and functional relatives of class I chitinase, the adaptive replacements of class I chitinase occur disproportionately in the active site cleft. This highly unusual pattern of replacements suggests that fungi directly defend against chitinolytic activity through enzymatic inhibition or other forms of chemical resistance and identifies target residues for manipulating chitinolytic activity. These data also provide empirical evidence that plant defense proteins not involved in pathogen recognition also evolve in a manner consistent with rapid coevolutionary interactions.

Possible ecological risks of transgenic organism release when transgenes affect mating success: Sexual selection and the Trojan gene hypothesis

Widespread interest in producing transgenic organisms is balanced by concern over ecological hazards, such as species extinction if such organisms were to be released into nature. An ecological risk associated with the introduction of a transgenic organism is that the transgene, though rare, can spread in a natural population. An increase in transgene frequency is often assumed to be unlikely because transgenic organisms typically have some viability disadvantage. Reduced viability is assumed to be common because transgenic individuals are best viewed as macromutants that lack any history of selection that could reduce negative fitness effects. However, these arguments ignore the potential advantageous effects of transgenes on some aspect of fitness such as mating success. Here, we examine the risk to a natural population after release of a few transgenic individuals when the transgene trait simultaneously increases transgenic male mating success and lowers the viability of transgenic offspring. We obtained relevant life history data by using the small cyprinodont fish, Japanese medaka (Oryzias latipes) as a model. Our deterministic equations predict that a transgene introduced into a natural population by a small number of transgenic fish will spread as a result of enhanced mating advantage, but the reduced viability of offspring will cause eventual local extinction of both populations. Such risks should be evaluated with each new transgenic animal before release.

DNA variation at the Sod locus of Drosophila melanogaster: An unfolding story of natural selection

Patterns of variation at the Sod locus of Drosophila melanogaster suggest that the protein polymorphism at this locus has very recently arisen. In addition, it appears that a previously rare DNA variant has been recently and rapidly driven to intermediate frequency. From the size of the region (>20 kb) that has been swept along with this rare variant, and patterns of linkage disequilibrium in the region, it is inferred that strength of selection was large (s > 0.01) and that the sweep occurred more than 25,000 generations ago. In addition, there are striking similarities to patterns of variation observed at the Est6 and Est-P loci, which are located approximately 1,000 kb from Sod.

Using in vitro selection to direct the covalent attachment of human immunodeficiency virus type 1 Rev protein to high-affinity RNA ligands.

We have used an in vitro selection procedure called crosslinking SELEX (SELEX = systematic evolution of ligands by exponential enrichment) to identify RNA sequences that bind with high affinity and crosslink to the Rev protein from human immunodeficiency virus type 1 (HIV-1). A randomized RNA library substituted with the photoreactive chromophore 5-iodouracil was irradiated with monochromatic UV light in the presence of Rev. Those sequences with the ability to photocrosslink to Rev were partitioned from the rest of the RNA pool, amplified, and used for the next round of selection. Rounds of photocrosslinking selection were alternated with rounds of selection for RNA sequences with high affinity to Rev. This iterative, dual-selection method yielded RNA molecules with subnanomolar dissociation constants and high efficiency photocrosslinking to Rev. Some of the RNA molecules isolated by this procedure form a stable complex with Rev that is resistant to denaturing gel electrophoresis in the absence of UV irradiation. In vitro selection of nucleic acids by using modified nucleotides allows the isolation of nucleic acid molecules with potentially limitless chemical capacities to covalently attack a target molecule.

Structure and function of the Bacillus hybrid enzyme GluXyn-1: Native-like jellyroll fold preserved after insertion of autonomous globular domain

The 1,3–1,4-β-glucanase from Bacillus macerans (wtGLU) and the 1,4-β-xylanase from Bacillus subtilis (wtXYN) are both single-domain jellyroll proteins catalyzing similar enzymatic reactions. In the fusion protein GluXyn-1, the two proteins are joined by insertion of the entire XYN domain into a surface loop of cpMAC-57, a circularly permuted variant of wtGLU. GluXyn-1 was generated by protein engineering methods, produced in Escherichia coli and shown to fold spontaneously and have both enzymatic activities at wild-type level. The crystal structure of GluXyn-1 was determined at 2.1 Å resolution and refined to R = 17.7% and R(free) = 22.4%. It shows nearly ideal, native-like folding of both protein domains and a small, but significant hinge bending between the domains. The active sites are independent and accessible explaining the observed enzymatic activity. Because in GluXyn-1 the complete XYN domain is inserted into the compact folding unit of GLU, the wild-type-like activity and tertiary structure of the latter proves that the folding process of GLU does not depend on intramolecular interactions that are short-ranged in the sequence. Insertion fusions of the GluXyn-1 type may prove to be an easy route toward more stable bifunctional proteins in which the two parts are more closely associated than in linear end-to-end protein fusions.

Glu-tRNAGln amidotransferase: A novel heterotrimeric enzyme required for correct decoding of glutamine codons during translation

The three genes, gatC, gatA, and gatB, which constitute the transcriptional unit of the Bacillus subtilis glutamyl-tRNAGln amidotransferase have been cloned. Expression of this transcriptional unit results in the production of a heterotrimeric protein that has been purified to homogeneity. The enzyme furnishes a means for formation of correctly charged Gln-tRNAGln through the transamidation of misacylated Glu-tRNAGln, functionally replacing the lack of glutaminyl-tRNA synthetase activity in Gram-positive eubacteria, cyanobacteria, Archaea, and organelles. Disruption of this operon is lethal. This demonstrates that transamidation is the only pathway to Gln-tRNAGln in B. subtilis and that glutamyl-tRNAGln amidotransferase is a novel and essential component of the translational apparatus.

Recurrent invasion and extinction of a selfish gene

Homing endonuclease genes show super-Mendelian inheritance, which allows them to spread in populations even when they are of no benefit to the host organism. To test the idea that regular horizontal transmission is necessary for the long-term persistence of these genes, we surveyed 20 species of yeasts for the ω-homing endonuclease gene and associated group I intron. The status of ω could be categorized into three states (functional, nonfunctional, or absent), and status was not clustered on the host phylogeny. Moreover, the phylogeny of ω differed significantly from that of the host, strong evidence of horizontal transmission. Further analyses indicate that horizontal transmission is more common than transposition, and that it occurs preferentially between closely related species. Parsimony analysis and coalescent theory suggest that there have been 15 horizontal transmission events in the ancestry of our yeast species, through simulations indicate that this value is probably an underestimate. Overall, the data support a cyclical model of invasion, degeneration, and loss, followed by reinvasion, and each of these transitions is estimated to occur about once every 2 million years. The data are thus consistent with the idea that frequent horizontal transmission is necessary for the long-term persistence of homing endonuclease genes, and further, that this requirement limits these genes to organisms with easily accessible germ lines. The data also show that mitochondrial DNA sequences are transferred intact between yeast species; if other genes do not show such high levels of horizontal transmission, it would be due to lack of selection, rather than lack of opportunity.

Somatic hypermutation of the new antigen receptor gene (NAR) in the nurse shark does not generate the repertoire: Possible role in antigen-driven reactions in the absence of germinal centers

The new antigen receptor (NAR) gene in the nurse shark diversifies extensively by somatic hypermutation. It is not known, however, whether NAR somatic hypermutation generates the primary repertoire (like in the sheep) or rather is used in antigen-driven immune responses. To address this issue, the sequences of NAR transmembrane (Tm) and secretory (Sec) forms, presumed to represent the primary and secondary repertoires, respectively, were examined from the peripheral blood lymphocytes of three adult nurse sharks. More than 40% of the Sec clones but fewer than 11% of Tm clones contained five mutations or more. Furthermore, more than 75% of the Tm clones had few or no mutations. Mutations in the Sec clones occurred mostly in the complementarity-determining regions (CDR) with a significant bias toward replacement substitutions in CDR1; in Tm clones there was no significant bias toward replacements and only a low level of targeting to the CDRs. Unlike the Tm clones where the replacement mutational pattern was similar to that seen for synonymous changes, Sec replacements displayed a distinct pattern of mutations. The types of mutations in NAR were similar to those found in mouse Ig genes rather than to the unusual pattern reported for shark and Xenopus Ig. Finally, an oligoclonal family of Sec clones revealed a striking trend toward acquisition of glutamic/aspartic acid, suggesting some degree of selection. These data strongly suggest that hypermutation of NAR does not generate the repertoire, but instead is involved in antigen-driven immune responses.

In vitro selection for altered divalent metal specificity in the RNase P RNA

The ribozyme RNase P absolutely requires divalent metal ions for catalytic function. Multiple Mg2+ ions contribute to the optimal catalytic efficiency of RNase P, and it is likely that the tertiary structure of the ribozyme forms a specific metal-binding pocket for these ions within the active-site. To identify base moieties that contribute to catalytic metal-binding sites, we have used in vitro selection to isolate variants of the Escherichia coli RNase P RNA with altered specificities for divalent metal. RNase P RNA variants with increased activity in Ca2+ were enriched over 18 generations of selection for catalysis in the presence of Ca2+, which is normally disfavored relative to Mg2+. Although a wide spectrum of mutations was found in the generation-18 clones, only a single point mutation was common to all clones: a cytosine-to-uracil transition at position 70 (E. coli numbering) of RNase P. Analysis of the C70U point mutant in a wild-type background confirmed that the identity of the base at position 70 is the sole determinant of Ca2+ selectivity. It is noteworthy that C70 lies within the phylogenetically well conserved J3/4-P4-J2/4 region, previously implicated in Mg2+ binding. Our finding that a single base change is sufficient to alter the metal preference of RNase P is further evidence that the J3/4-P4-J2/4 domain forms a portion of the ribozyme’s active site.

Strong balancing selection at HLA loci: Evidence from segregation in South Amerindian families

The genotypic proportions for major histocompatibility complex loci, HLA-A and HLA-B, of progeny in families in 23 South Amerindian tribes in which segregation for homozygotes and heterozygotes could occur are examined. Overall, there is a large deficiency of homozygotes compared with Mendelian expectations (for HLA-A, 114 observed and 155.50 expected and for HLA-B 110 observed and 144.75 expected), consistent with strong balancing selection favoring heterozygotes. There is no evidence that these deficiencies were associated with particular alleles or with the age of the individuals sampled. When these families were divided into four mating types, there was strong selection against homozygotes, averaging 0.462 for three of the mating types over the two loci. For the other mating type in which the female parent is homozygous and shares one allele with the heterozygous male parent, there was no evidence of selection against homozygotes. A theoretical model incorporating these findings surprisingly does not result in a stable polymorphism for two alleles but does result in an excess of heterozygotes and a minimum fitness at intermediate allele frequencies. However, for more than two alleles, balancing selection does occur and the model approaches the qualities of the symmetrical heterozygote advantage model as the number of alleles increases.

Role of left inferior prefrontal cortex in retrieval of semantic knowledge: A reevaluation

A number of neuroimaging findings have been interpreted as evidence that the left inferior frontal gyrus (IFG) subserves retrieval of semantic knowledge. We provide a fundamentally different interpretation, that it is not retrieval of semantic knowledge per se that is associated with left IFG activity but rather selection of information among competing alternatives from semantic memory. Selection demands were varied across three semantic tasks in a single group of subjects. Functional magnetic resonance imaging signal in overlapping regions of left IFG was dependent on selection demands in all three tasks. In addition, the degree of semantic processing was varied independently of selection demands in one of the tasks. The absence of left IFG activity for this comparison counters the argument that the effects of selection can be attributed solely to variations in degree of semantic retrieval. Our findings suggest that it is selection, not retrieval, of semantic knowledge that drives activity in the left IFG.

Spectral karyotyping suggests additional subsets of colorectal cancers characterized by pattern of chromosome rearrangement

The abundant chromosome abnormalities in most carcinomas are probably a reflection of genomic instability present in the tumor, so the pattern and variability of chromosome abnormalities will reflect the mechanism of instability combined with the effects of selection. Chromosome rearrangement was investigated in 17 colorectal carcinoma-derived cell lines. Comparative genomic hybridization showed that the chromosome changes were representative of those found in primary tumors. Spectral karyotyping (SKY) showed that translocations were very varied and mostly unbalanced, with no translocation occurring in more than three lines. At least three karyotype patterns could be distinguished. Some lines had few chromosome abnormalities: they all showed microsatellite instability, the replication error (RER)+ phenotype. Most lines had many chromosome abnormalities: at least seven showed a surprisingly consistent pattern, characterized by multiple unbalanced translocations and intermetaphase variation, with chromosome numbers around triploid, 6–16 structural aberrations, and similarities in gains and losses. Almost all of these were RER−, but one, LS411, was RER+. The line HCA7 showed a novel pattern, suggesting a third kind of genomic instability: multiple reciprocal translocations, with little numerical change or variability. This line was also RER+. The coexistence in one tumor of two kinds of genomic instability is to be expected if the underlying defects are selected for in tumor evolution.

Does recombination improve selection on codon usage? Lessons from nematode and fly complete genomes

Understanding the factors responsible for variations in mutation patterns and selection efficacy along chromosomes is a prerequisite for deciphering genome sequences. Population genetics models predict a positive correlation between the efficacy of selection at a given locus and the local rate of recombination because of Hill–Robertson effects. Codon usage is considered one of the most striking examples that support this prediction at the molecular level. In a wide range of species including Caenorhabditis elegans and Drosophila melanogaster, codon usage is essentially shaped by selection acting for translational efficiency. Codon usage bias correlates positively with recombination rate in Drosophila, apparently supporting the hypothesis that selection on codon usage is improved by recombination. Here we present an exhaustive analysis of codon usage in C. elegans and D. melanogaster complete genomes. We show that in both genomes there is a positive correlation between recombination rate and the frequency of optimal codons. However, we demonstrate that in both species, this effect is due to a mutational bias toward G and C bases in regions of high recombination rate, possibly as a direct consequence of the recombination process. The correlation between codon usage bias and recombination rate in these species appears to be essentially determined by recombination-dependent mutational patterns, rather than selective effects. This result highlights that it is necessary to take into account the mutagenic effect of recombination to understand the evolutionary role and impact of recombination.