In vivo selection of basic region-leucine zipper proteins with altered DNA-binding specificities.

A transcription interference assay was used to generate mutant basic region-leucine zipper proteins with altered DNA-binding specificities. A library of mutants of a CCAAT/enhancer binding protein was constructed by randomizing five DNA-contacting amino acids in the basic region Asn-18, Ala-15, Val-14, Ser-11, and Arg-10. These mutants were then selected for their ability to bind mutant recognition sequences containing substitutions at the 2 and 3 positions of the wild-type sequence 5'-A5T4T3G2C1G1'C2'A3A4'T5'-3'. Mutants containing the sequence Leu-18Tyr-15Xaa-14Tyr-11Arg-10, in which four of the five contact residues are altered, were identified that recognize the palindromic sequence 5'-ATCYCGY'GAT-3' (Xaa = asparagine when Y = G; Xaa = methionine when Y = A). Moreover, in a selection against the sequence 5'-ATTACGTAAT-3', mutants were obtained containing substitutions not only in the basic region but also in the hinge region between the basic and leucine zipper regions. The mutant proteins showed high specificity in a functional transcription interference assay. A model for the interaction of these mutants with the target DNA sequences is discussed.

Distribution of haplotypes from a chromosome 21 region distinguishes multiple prehistoric human migrations

Despite mounting genetic evidence implicating a recent origin of modern humans, the elucidation of early migratory gene-flow episodes remains incomplete. Geographic distribution of haplotypes may show traces of ancestral migrations. However, such evolutionary signatures can be erased easily by recombination and mutational perturbations. A 565-bp chromosome 21 region near the MX1 gene, which contains nine sites frequently polymorphic in human populations, has been found. It is unaffected by recombination and recurrent mutation and thus reflects only migratory history, genetic drift, and possibly selection. Geographic distribution of contemporary haplotypes implies distinctive prehistoric human migrations: one to Oceania, one to Asia and subsequently to America, and a third one predominantly to Europe. The findings with chromosome 21 are confirmed by independent evidence from a Y chromosome phylogeny. Loci of this type will help to decipher the evolutionary history of modern humans.

Evidence for balancing selection operating at the het-c heterokaryon incompatibility locus in a group of filamentous fungi

In filamentous fungi, het loci (for heterokaryon incompatibility) are believed to regulate self/nonself-recognition during vegetative growth. As filamentous fungi grow, hyphal fusion occurs within an individual colony to form a network. Hyphal fusion can occur also between different individuals to form a heterokaryon, in which genetically distinct nuclei occupy a common cytoplasm. However, heterokaryotic cells are viable only if the individuals involved have identical alleles at all het loci. One het locus, het-c, has been characterized at the molecular level in Neurospora crassa and encodes a glycine-rich protein. In an effort to understand the role of this locus in filamentous fungi, we chose to study its evolution by analyzing het-c sequence variability in species within Neurospora and related genera. We determined that the het-c locus was polymorphic in a field population of N. crassa with close to equal frequency of each of the three allelic types. Different species and even genera within the Sordariaceae shared het-c polymorphisms, indicating that these polymorphisms originated in an ancestral species. Finally, an analysis of the het-c specificity region shows a high occurrence of nonsynonymous substitution. The persistence of allelic lineages, the nearly equal allelic distribution within populations, and the high frequency of nonsynonymous substitutions in the het-c specificity region suggest that balancing selection has operated to maintain allelic diversity at het-c. Het-c shares this particular evolutionary characteristic of departing from neutrality with other self/nonself-recognition systems such as major histocompatibility complex loci in mammals and the S (self-incompatibility) locus in angiosperms.

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.

Emergence of FY*Anull in a Plasmodium vivax-endemic region of Papua New Guinea

In Papua New Guinea (PNG), numerous blood group polymorphisms and hemoglobinopathies characterize the human population. Human genetic polymorphisms of this nature are common in malarious regions, and all four human malaria parasites are holoendemic below 1500 meters in PNG. At this elevation, a prominent condition characterizing Melanesians is α+-thalassemia. Interestingly, recent epidemiological surveys have demonstrated that α+-thalassemia is associated with increased susceptibility to uncomplicated malaria among young children. It is further proposed that α+-thalassemia may facilitate so-called “benign” Plasmodium vivax infection to act later in life as a “natural vaccine” against severe Plasmodium falciparum malaria. Here, in a P. vivax-endemic region of PNG where the resident Abelam-speaking population is characterized by a frequency of α+-thalassemia ≥0.98, we have discovered the mutation responsible for erythrocyte Duffy antigen-negativity (Fy[a−b−]) on the FY*A allele. In this study population there were 23 heterozygous and no homozygous individuals bearing this new allele (allele frequency, 23/1062 = 0.022). Flow cytometric analysis illustrated a 2-fold difference in erythroid-specific Fy-antigen expression between heterozygous (FY*A/FY*Anull) and homozygous (FY*A/FY*A) individuals, suggesting a gene-dosage effect. In further comparisons, we observed a higher prevalence of P. vivax infection in FY*A/FY*A (83/508 = 0.163) compared with FY*A/FY*Anull (2/23 = 0.087) individuals (odds ratio = 2.05, 95% confidence interval = 0.47–8.91). Emergence of FY*Anull in this population suggests that P. vivax is involved in selection of this erythroid polymorphism. This mutation would ultimately compromise α+-thalassemia/P. vivax-mediated protection against severe P. falciparum malaria.

Spectral bifurcations in dispersive wave turbulence

Dispersive wave turbulence is studied numerically for a class of one-dimensional nonlinear wave equations. Both deterministic and random (white noise in time) forcings are studied. Four distinct stable spectra are observed—the direct and inverse cascades of weak turbulence (WT) theory, thermal equilibrium, and a fourth spectrum (MMT; Majda, McLaughlin, Tabak). Each spectrum can describe long-time behavior, and each can be only metastable (with quite diverse lifetimes)—depending on details of nonlinearity, forcing, and dissipation. Cases of a long-live MMT transient state dcaying to a state with WT spectra, and vice-versa, are displayed. In the case of freely decaying turbulence, without forcing, both cascades of weak turbulence are observed. These WT states constitute the clearest and most striking numerical observations of WT spectra to date—over four decades of energy, and three decades of spatial, scales. Numerical experiments that study details of the composition, coexistence, and transition between spectra are then discussed, including: (i) for deterministic forcing, sharp distinctions between focusing and defocusing nonlinearities, including the role of long wavelength instabilities, localized coherent structures, and chaotic behavior; (ii) the role of energy growth in time to monitor the selection of MMT or WT spectra; (iii) a second manifestation of the MMT spectrum as it describes a self-similar evolution of the wave, without temporal averaging; (iv) coherent structures and the evolution of the direct and inverse cascades; and (v) nonlocality (in k-space) in the transferral process.

RNA polymerase sigma factor determines start-site selection but is not required for upstream promoter element activation on heteroduplex (bubble) templates

Sequence-selective transcription by bacterial RNA polymerase (RNAP) requires σ factor that participates in both promoter recognition and DNA melting. RNAP lacking σ (core enzyme) will initiate RNA synthesis from duplex ends, nicks, gaps, and single-stranded regions. We have used DNA templates containing short regions of heteroduplex (bubbles) to compare initiation in the presence and absence of various σ factors. Using bubble templates containing the σD-dependent flagellin promoter, with or without its associated upstream promoter (UP) element, we demonstrate that UP element stimulation occurs efficiently even in the absence of σ. This supports a model in which the UP element acts primarily through the α subunit of core enzyme to increase the initial association of RNAP with the promoter. Core and holoenzyme do differ substantially in the template positions chosen for initiation: σD restricts initiation to sites 8–9 nucleotides downstream of the conserved −10 element. Remarkably, σA also has a dramatic effect on start-site selection even though the σA holoenzyme is inactive on the corresponding homoduplexes. The start sites chosen by the σA holoenzyme are located 8 nucleotides downstream of sequences on the nontemplate strand that resemble the conserved −10 hexamer recognized by σA. Thus, σA appears to recognize the −10 region even in a single-stranded state. We propose that in addition to its described roles in promoter recognition and start-site melting, σ also localizes the transcription start site.

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.

Novel gene conversion between X-Y homologues located in the nonrecombining region of the Y chromosome in Felidae (Mammalia)

Genes located on the mammalian Y chromosome outside of the pseudoautosomal region do not recombine with those on the X and are predicted to either undergo selection for male function or gradually degenerate because of an accumulation of deleterious mutations. Here, phylogenetic analyses of X-Y homologues, Zfx and Zfy, among 26 felid species indicate two ancestral episodes of directed genetic exchange (ectopic gene conversion) from X to Y: once during the evolution of pallas cat and once in a common predecessor of ocelot lineage species. Replacement of the more rapidly evolving Y homologue with the evolutionarily constrained X copy may represent a mechanism for adaptive editing of functional genes on the nonrecombining region of the mammalian Y chromosome.

A first trial of retrospective collaboration for positional cloning in complex inheritance: Assay of the cytokine region on chromosome 5 by the Consortium on Asthma Genetics (COAG)

The central problem of complex inheritance is to map oligogenes for disease susceptibility, integrating linkage and association over samples that differ in several ways. Combination of evidence over multiple samples with 1,037 families supports loci contributing to asthma susceptibility in the cytokine region on 5q [maximum logarithm of odds (lod) = 2.61 near IL-4], but no evidence for atopy. The principal problems with retrospective collaboration on linkage appear to have been solved, providing far more information than a single study. A multipoint lod table evaluated at commonly agreed reference loci is required for both collaboration and metaanalysis, but variations in ascertainment, pedigree structure, phenotype definition, and marker selection are tolerated. These methods are invariant with statistical methods that increase the power of lods and are applicable to all diseases, motivating collaboration rather than competition. In contrast to linkage, positional cloning by allelic association has yet to be extended to multiple samples, a prerequisite for efficient combination with linkage and the greatest current challenge to genetic epidemiology.

Normal cardiovascular development in mice deficient for 16 genes in 550 kb of the velocardiofacial/ DiGeorge syndrome region

Hemizygous interstitial deletions in human chromosome 22q11 are associated with velocardiofacial syndrome and DiGeorge syndrome and lead to multiple congenital abnormalities, including cardiovascular defects. The gene(s) responsible for these disorders is thought to reside in a 1.5-Mb region of 22q11 in which 27 genes have been identified. We have used Cre-mediated recombination of LoxP sites in embryonic stem cells and mice to generate a 550-kb deletion encompassing 16 of these genes in the corresponding region on mouse chromosome 16. Mice heterozygous for this deletion are normal and do not exhibit cardiovascular abnormalities. Because mice with a larger deletion on mouse chromosome 16 do have heart defects, the results allow us to exclude these 16 genes as being solely, or in combination among themselves, responsible for the cardiovascular abnormalities in velocardiofacial/DiGeorge syndrome. We also generated mice with a duplication of the 16 genes that may help dissect the genetic basis of “cat eye” and derivative 22 syndromes that are characterized by extra copies of portions of 22q11, including these 16 genes. We also describe a strategy for selecting cell lines with defined chromosomal rearrangements. The method is based on reconstitution of a dominant selection marker after Cre-mediated recombination of LoxP sites. Therefore it should be widely applicable to many cell lines.

How Escherichia coli can bias the results of molecular cloning: Preferential selection of defective genomes of hepatitis C virus during the cloning procedure

Cloned PCR products containing hepatitis C virus (HCV) genomic fragments have been used for analyses of HCV genomic heterogeneity and protein expression. These studies assume that the clones derived are representative of the entire virus population and that subsets are not inadvertently selected. The aim of the present study was to express HCV structural proteins. However, we found that there was a strong cloning selection for defective genomes and that most clones generated initially were incapable of expressing the HCV proteins. The HCV structural region (C-E1-E2-p7) was directly amplified by long reverse transcription–PCR from the plasma of an HCV-infected patient or from a control plasmid containing a viable full-length cDNA of HCV derived from the same patient but cloned in a different vector. The PCR products were cloned into a mammalian expression vector, amplified in Escherichia coli, and tested for their ability to produce HCV structural proteins. Twenty randomly picked clones derived from the HCV-infected patient all contained nucleotide mutations leading to absence or truncation of the expected HCV products. Of 25 clones derived from the control plasmid, only 8% were fully functional for polyprotein synthesis. The insertion of extra nucleotides in the region just upstream of the start codon of the HCV insert led to a statistically significant increase in the number of fully functional clones derived from the patient (42%) and from the control plasmid (72–92%). Nonrandom selection of clones during the cloning procedure has enormous implications for the study of viral heterogeneity, because it can produce a false spectrum of genomic diversity. It can also be an impediment to the construction of infectious viral clones.

Positive Darwinian selection drives the evolution of several female reproductive proteins in mammals

Rapid evolution driven by positive Darwinian selection is a recurrent theme in male reproductive protein evolution. In contrast, positive selection has never been demonstrated for female reproductive proteins. Here, we perform phylogeny-based tests on three female mammalian fertilization proteins and demonstrate positive selection promoting their divergence. Two of these female fertilization proteins, the zona pellucida glycoproteins ZP2 and ZP3, are part of the mammalian egg coat. Several sites identified in ZP3 as likely to be under positive selection are located in a region previously demonstrated to be involved in species-specific sperm-egg interaction, suggesting the selective pressure is related to male-female interaction. The results provide long-sought evidence for two evolutionary hypotheses: sperm competition and sexual conflict.

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.

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.