Sequence variation in the Fanconi anemia gene FAA

Fanconi anemia (FA) is a genetically heterogeneous autosomal recessive syndrome associated with chromosomal instability, hypersensitivity to DNA crosslinking agents, and predisposition to malignancy. The gene for FA complementation group A (FAA) recently has been cloned. The cDNA is predicted to encode a polypeptide of 1,455 amino acids, with no homologies to any known protein that might suggest a function for FAA. We have used single-strand conformational polymorphism analysis to screen genomic DNA from a panel of 97 racially and ethnically diverse FA patients from the International Fanconi Anemia Registry for mutations in the FAA gene. A total of 85 variant bands were detected. Forty-five of the variants are probably benign polymorphisms, of which nine are common and can be used for various applications, including mapping studies for other genes in this region of chromosome 16q. Amplification refractory mutation system assays were developed to simplify their detection. Forty variants are likely to be pathogenic mutations. Seventeen of these are microdeletions/microinsertions associated with short direct repeats or homonucleotide tracts, a type of mutation thought to be generated by a mechanism of slipped-strand mispairing during DNA replication. A screening of 350 FA probands from the International Fanconi Anemia Registry for two of these deletions (1115–1118del and 3788–3790del) revealed that they are carried on about 2% and 5% of the FA alleles, respectively. 3788–3790del appears in a variety of ethnic groups and is found on at least two different haplotypes. We suggest that FAA is hypermutable, and that slipped-strand mispairing, a mutational mechanism recognized as important for the generation of germ-line and somatic mutations in a variety of cancer-related genes, including p53, APC, RB1, WT1, and BRCA1, may be a major mechanism for FAA mutagenesis.

Striking sequence similarity in inter- and intra-specific comparisons of class I SLG alleles from Brassica oleracea and Brassica campestris: Implications for the evolution and recognition mechanism

Self-incompatibility in Brassica is controlled by a single multi-allelic locus (S locus), which contains at least two highly polymorphic genes expressed in the stigma: an S glycoprotein gene (SLG) and an S receptor kinase gene (SRK). The putative ligand-binding domain of SRK exhibits high homology to the secretory protein SLG, and it is believed that SLG and SRK form an active receptor kinase complex with a self-pollen ligand, which leads to the rejection of self-pollen. Here, we report 31 novel SLG sequences of Brassica oleracea and Brassica campestris. Sequence comparisons of a large number of SLG alleles and SLG-related genes revealed the following points. (i) The striking sequence similarity observed in an inter-specific comparison (95.6% identity between SLG14 of B. oleracea and SLG25 of B. campestris in deduced amino acid sequence) suggests that SLG diversification predates speciation. (ii) A perfect match of the sequences in hypervariable regions, which are thought to determine S specificity in an intra-specific comparison (SLG8 and SLG46 of B. campestris) and the observation that the hypervariable regions of SLG and SRK of the same S haplotype were not necessarily highly similar suggests that SLG and SRK bind different sites of the pollen ligand and that they together determine S specificity. (iii) Comparison of the hypervariable regions of SLG alleles suggests that intragenic recombination, together with point mutations, has contributed to the generation of the high level of sequence variation in SLG alleles. Models for the evolution of SLG/SRK are presented.

Human Immunodeficiency Virus Reverse Transcriptase and Protease Sequence Database: an expanded data model integrating natural language text and sequence analysis programs

The HIV Reverse Transcriptase and Protease Sequence Database is an on-line relational database that catalogs evolutionary and drug-related sequence variation in the human immunodeficiency virus (HIV) reverse transcriptase (RT) and protease enzymes, the molecular targets of anti-HIV therapy (http://hivdb.stanford.edu). The database contains a compilation of nearly all published HIV RT and protease sequences, including submissions from International Collaboration databases and sequences published in journal articles. Sequences are linked to data about the source of the sequence sample and the antiretroviral drug treatment history of the individual from whom the isolate was obtained. During the past year 3500 sequences have been added and the data model has been expanded to include drug susceptibility data on sequenced isolates. Database content has also been integrated with didactic text and the output of two sequence analysis programs.

Characterization of the CHD family of proteins

The murine gene CHD1 (MmCHD1) was previously isolated in a search for proteins that bound a DNA promoter element. The presence of chromo (chromatin organization modifier) domains and an SNF2-related helicase/ATPase domain led to speculation that this gene regulated chromatin structure or gene transcription. This study describes the cloning and characterization of three novel human genes related to MmCHD1. Examination of sequence databases produced several more related genes, most of which were not known to be similar to MmCHD1, yielding a total of 12 highly conserved CHD genes from organisms as diverse as yeast and mammals. The major region of sequence variation is in the C-terminal part of the protein, a region with DNA-binding activity in MmCHD1. Targeted deletion of ScCHD1, the sole Saccharomyces cerevesiae CHD gene, was performed with deletion strains being less sensitive than wild type to the cytotoxic effect of 6-azauracil. This finding suggested that enhanced transcriptional arrest at RNA polymerase II pause sites due to 6-azauracil-induced nucleotide pool depletion was reduced in the deletion strain and that ScCHD1 inhibited transcription. This observation, along with the known roles of other proteins with chromo or SNF2-related helicase/ATPase domains, suggests that alteration of gene expression by CHD genes might occur by modifications of chromatin structure, with altered access of the transcriptional apparatus to its chromosomal DNA template.

Lysis of HIV-1-infected cells and inhibition of viral replication by universal receptor T cells

Increasing evidence suggests that HIV-1-specific cytotoxic T lymphocytes (CTLs) are a key host immune response to HIV-1 infection. Generation of CTL responses for prevention or therapy of HIV-1 infection has several intrinsic technical barriers such as antigen expression and presentation, the varying HLA restrictions between different individuals, and the potential for viral escape by sequence variation or surface molecule alteration on infected cells. A strategy to circumvent these limitations is the construction of a chimeric T cell receptor containing human CD4 or HIV-1-specific Ig sequences linked to the signaling domain of the T cell receptor ζ chain (universal T cell receptor). CD8+ CTLs transduced with this universal receptor can then bind and lyse infected cells that express surface HIV-1 gp120. We evaluated the ability of universal-receptor-bearing CD8+ cells from a seronegative donor to lyse acutely infected cells and inhibit HIV-1 replication in vitro. The kinetics of lysis and efficiency of inhibition were comparable to that of naturally occurring HIV-1-specific CTL clones isolated from infected individuals. Further study will be required to determine the utility of these cells as a therapeutic strategy in vivo.

Population genetics of Ice Age brown bears

The Pleistocene was a dynamic period for Holarctic mammal species, complicated by episodes of glaciation, local extinctions, and intercontinental migration. The genetic consequences of these events are difficult to resolve from the study of present-day populations. To provide a direct view of population genetics in the late Pleistocene, we measured mitochondrial DNA sequence variation in seven permafrost-preserved brown bear (Ursus arctos) specimens, dated from 14,000 to 42,000 years ago. Approximately 36,000 years ago, the Beringian brown bear population had a higher genetic diversity than any extant North American population, but by 15,000 years ago genetic diversity appears similar to the modern day. The older, genetically diverse, Beringian population contained sequences from three clades now restricted to local regions within North America, indicating that current phylogeographic patterns may provide misleading data for evolutionary studies and conservation management. The late Pleistocene phylogeographic data also indicate possible colonization routes to areas south of the Cordilleran ice sheet.

Mitochondrial diversity and the origins of African and European cattle.

The nature of domestic cattle origins in Africa are unclear as archaeological data are relatively sparse. The earliest domesticates were humpless, or Bos taurus, in morphology and may have shared a common origin with the ancestors of European cattle in the Near East. Alternatively, local strains of the wild ox, the aurochs, may have been adopted by peoples in either continent either before or after cultural influence from the Levant. This study examines mitochondrial DNA displacement loop sequence variation in 90 extant bovines drawn from Africa, Europe, and India. Phylogeny estimation and analysis of molecular variance verify that sequences cluster significantly into continental groups. The Indian Bos indicus samples are most markedly distinct from the others, which is indicative of a B. taurus nature for both European and African ancestors. When a calibration of sequence divergence is performed using comparisons with bison sequences and an estimate of 1 Myr since the Bison/Bos Leptobos common ancestor, estimates of 117-275,000 B.P. and 22-26,000 B.P. are obtained for the separation between Indians and others and between African and European ancestors, respectively. As cattle domestication is thought to have occurred approximately 10,000 B.P., these estimates suggest the domestication of genetically discrete aurochsen strains as the origins of each continental population. Additionally, patterns of variation that are indicative of population expansions (probably associated with the domestication process) are discernible in Africa and Europe. Notably, the genetic signatures of these expansions are clearly younger than the corresponding signature of African/European divergence.

DNA sequence insertion and evolutionary variation in gene regulation.

Current evidence on the long-term evolutionary effect of insertion of sequence elements into gene regions is reviewed, restricted to cases where a sequence derived from a past insertion participates in the regulation of expression of a useful gene. Ten such examples in eukaryotes demonstrate that segments of repetitive DNA or mobile elements have been inserted in the past in gene regions, have been preserved, sometimes modified by selection, and now affect control of transcription of the adjacent gene. Included are only examples in which transcription control was modified by the insert. Several cases in which merely transcription initiation occurred in the insert were set aside. Two of the examples involved the long terminal repeats of mammalian endogenous retroviruses. Another two examples were control of transcription by repeated sequence inserts in sea urchin genomes. There are now six published examples in which Alu sequences were inserted long ago into human gene regions, were modified, and now are central in control/enhancement of transcription. The number of published examples of Alu sequences affecting gene control has grown threefold in the last year and is likely to continue growing. Taken together, all of these examples show that the insertion of sequence elements in the genome has been a significant source of regulatory variation in evolution.

Evolution of plastid gene rps2 in a lineage of hemiparasitic and holoparasitic plants: Many losses of photosynthesis and complex patterns of rate variation

The plastid genomes of some nonphotosynthetic parasitic plants have experienced an extreme reduction in gene content and an increase in evolutionary rate of remaining genes. Nothing is known of the dynamics of these events or whether either is a direct outcome of the loss of photosynthesis. The parasitic Scrophulariaceae and Orobanchaceae, representing a continuum of heterotrophic ability ranging from photosynthetic hemiparasites to nonphotosynthetic holoparasites, are used to investigate these issues. We present a phylogenetic hypothesis for parasitic Scrophulariaceae and Orobanchaceae based on sequences of the plastid gene rps2, encoding the S2 subunit of the plastid ribosome. Parasitic Scrophulariaceae and Orobanchaceae form a monophyletic group in which parasitism can be inferred to have evolved once. Holoparasitism has evolved independently at least five times, with certain holoparasitic lineages representing single species, genera, and collections of nonphotosynthetic genera. Evolutionary loss of the photosynthetic gene rbcL is limited to a subset of holoparasitic lineages, with several holoparasites retaining a full length rbcL sequence. In contrast, the translational gene rps2 is retained in all plants investigated but has experienced rate accelerations in several hemi- as well as holoparasitic lineages, suggesting that there may be substantial molecular evolutionary changes to the plastid genome of parasites before the loss of photosynthesis. Independent patterns of synonymous and nonsynonymous rate acceleration in rps2 point to distinct mechanisms underlying rate variation in different lineages. Parasitic Scrophulariaceae (including the traditional Orobanchaceae) provide a rich platform for the investigation of molecular evolutionary process, gene function, and the evolution of parasitism.

Allelic variation of a dehydrin gene cosegregates with chilling tolerance during seedling emergence

Dehydrins (DHNs, LEA D-11) are plant proteins present during environmental stresses associated with dehydration or low temperatures and during seed maturation. Functions of DHNs have not yet been defined. Earlier, we hypothesized that a ≈35-kDa DHN and membrane properties that reduce electrolyte leakage from seeds confer chilling tolerance during seedling emergence of cowpea (Vigna unguiculata L. Walp.) in an additive and independent manner. Evidence for this hypothesis was not rigorous because it was based on correlations of presence/absence of the DHN and slow electrolyte leakage with chilling tolerance in closely related cowpea lines that have some other genetic differences. Here, we provide more compelling genetic evidence for involvement of the DHN in chilling tolerance of cowpea. We developed near-isogenic lines by backcrossing. We isolated and determined the sequence of a cDNA corresponding to the ≈35-kDa DHN and used gene-specific oligonucleotides derived from it to test the genetic linkage between the DHN presence/absence trait and the DHN structural gene. We tested for association between the DHN presence/absence trait and both low-temperature seed emergence and electrolyte leakage. We show that allelic differences in the Dhn structural gene map to the same position as the DHN protein presence/absence trait and that the presence of the ≈35-kDa DHN is indeed associated with chilling tolerance during seedling emergence, independent of electrolyte leakage effects. Two types of allelic variation in the Dhn gene were identified in the protein-coding region, deletion of one Φ-segment from the DHN-negative lines and two single amino acid substitutions.

Comparative studies on mammalian Hoxc8 early enhancer sequence reveal a baleen whale-specific deletion of a cis-acting element

Variations in regulatory regions of developmental control genes have been implicated in the divergence of axial morphologies. To find potentially significant changes in cis-regulatory regions, we compared nucleotide sequences and activities of mammalian Hoxc8 early enhancers. The nucleotide sequence of the early enhancer region is extremely conserved among mammalian clades, with five previously described cis-acting elements, A–E, being invariant. However, a 4-bp deletion within element C of the Hoxc8 early enhancer sequence is observed in baleen whales. When assayed in transgenic mouse embryos, a baleen whale enhancer (unlike other mammalian enhancers) directs expression of the reporter gene to more posterior regions of the neural tube but fails to direct expression to posterior mesoderm. We suggest that regulation of Hoxc8 in baleen whales differs from other mammalian species and may be associated with variation in axial morphology.

STACK: Sequence Tag Alignment and Consensus Knowledgebase

STACK is a tool for detection and visualisation of expressed transcript variation in the context of developmental and pathological states. The datasystem organises and reconstructs human transcripts from available public data in the context of expression state. The expression state of a transcript can include developmental state, pathological association, site of expression and isoform of expressed transcript. STACK consensus transcripts are reconstructed from clusters that capture and reflect the growing evidence of transcript diversity. The comprehensive capture of transcript variants is achieved by the use of a novel clustering approach that is tolerant of sub-sequence diversity and does not rely on pairwise alignment. This is in contrast with other gene indexing projects. STACK is generated at least four times a year and represents the exhaustive processing of all publicly available human EST data extracted from GenBank. This processed information can be explored through 15 tissue-specific categories, a disease-related category and a whole-body index and is accessible via WWW at http://www.sanbi.ac.za/Dbases.html. STACK represents a broadly applicable resource, as it is the only reconstructed transcript database for which the tools for its generation are also broadly available (http://www.sanbi.ac.za/CODES).

The nucleotide changes governing cuticular hydrocarbon variation and their evolution in Drosophila melanogaster

The cuticular hydrocarbon (CH) pheromones in Drosophila melanogaster exhibit strong geographic variation. African and Caribbean populations have a high ratio of 5,9 heptacosadiene/7,11 heptacosadiene (the “High” CH type), whereas populations from all other areas have a low ratio (“Low” CH type). Based on previous genetic mapping, DNA markers were developed that localized the genetic basis of this CH polymorphism to within a 13-kb region. We then carried out a hierarchical search for diagnostic nucleotide sites starting with four lines, and increasing to 24 and 43 lines from a worldwide collection. Within the 13-kb region, only one variable site shows a complete concordance with the CH phenotype. This is a 16-bp deletion in the 5′ region of a desaturase gene (desat2) that was recently suggested to be responsible for the CH polymorphism on the basis of its expression [Dallerac, R., Labeur, C., Jallon, J.-M., Knipple, D. C., Roelofs, W. L. & Wicker-Thomas, C. (2000) Proc. Natl. Acad. Sci. 97, 9449–9454]. The cosmopolitan Low type is derived from the ancestral High type, and DNA sequence variations suggest that the former spread worldwide with the aid of positive selection. Whether this CH variation could be a component of the sexual isolation between Zimbabwe and other cosmopolitan populations remains an interesting and unresolved question.

Naturally occurring variation in copia expression is due to both element (cis) and host (trans) regulatory variation.

Significant differences in levels of copia [Drosophila long terminal repeat (LTR) retrotransposon] expression exist among six species representing the Drosophila melanogaster species complex (D. melanogaster, Drosophila mauritiana, Drosophila simulans, Drosophila sechellia, Drosophila yakuba, and Drosophila erecta) and a more distantly related species (Drosophila willistoni). These differences in expression are correlated with major size variation mapping to putative regulatory regions of the copia 5' LTR and adjacent untranslated leader region (ULR). Sequence analysis indicates that these size variants were derived from a series of regional duplication events. The ability of the copia LTR-ULR size variants to drive expression of a bacterial chloramphenicol acetyltransferase reporter gene was tested in each of the seven species. The results indicate that both element-encoded (cis) and host-genome-encoded (trans) genetic differences are responsible for the variability in copia expression within and between Drosophila species. This finding indicates that models purporting to explain the dynamics and distribution of retrotransposons in natural populations must consider the potential impact of both element-encoded and host-genome-encoded regulatory variation to be valid. We propose that interelement selection among retrotransposons may provide a molecular drive mechanism for the evolution of eukaryotic enhancers which can be subsequently distributed throughout the genome by retrotransposition.