Highly specific protein sequence motifs for genome analysis

We present a method for discovering conserved sequence motifs from families of aligned protein sequences. The method has been implemented as a computer program called emotif (http://motif.stanford.edu/emotif). Given an aligned set of protein sequences, emotif generates a set of motifs with a wide range of specificities and sensitivities. emotif also can generate motifs that describe possible subfamilies of a protein superfamily. A disjunction of such motifs often can represent the entire superfamily with high specificity and sensitivity. We have used emotif to generate sets of motifs from all 7,000 protein alignments in the blocks and prints databases. The resulting database, called identify (http://motif.stanford.edu/identify), contains more than 50,000 motifs. For each alignment, the database contains several motifs having a probability of matching a false positive that range from 10−10 to 10−5. Highly specific motifs are well suited for searching entire proteomes, while generating very few false predictions. identify assigns biological functions to 25–30% of all proteins encoded by the Saccharomyces cerevisiae genome and by several bacterial genomes. In particular, identify assigned functions to 172 of proteins of unknown function in the yeast genome.

A gene family required for human germ cell development evolved from an ancient meiotic gene conserved in metazoans

The Deleted in AZoospermia (DAZ) genes encode potential RNA-binding proteins that are expressed exclusively in prenatal and postnatal germ cells and are strong candidates for human fertility factors. Here we report the identification of an additional member of the DAZ gene family, which we have called BOULE. With the identification of this gene, it is clear that the human DAZ gene family contains at least three members: DAZ, a Y-chromosome gene cluster that arose 30–40 million years ago and whose deletion is linked to infertility in men; DAZL, the “father” of DAZ, a gene that maps to human chromosome 3 and has homologs required for both female and male germ cell development in other organisms; and BOULE, a gene that we propose is the “grandfather” of DAZ and maps to human chromosome 2. Human and mouse BOULE resemble the invertebrate meiotic regulator Boule, the proposed ortholog of DAZ, in sequence and expression pattern and hence likely perform a similar meiotic function. In contrast, the previously identified human DAZ and DAZL are expressed much earlier than BOULE in prenatal germ stem cells and spermatogonia; DAZL also is expressed in female germ cells. These data suggest that homologs of the DAZ gene family can be grouped into two subfamilies (BOULE and DAZL) and that members of the DAZ family evolved from an ancestral meiotic regulator, Boule, to assume distinct, yet overlapping, functions in germ cell development.

Expression Patterns Conferred by Tyrosine/Dihydroxyphenylalanine Decarboxylase Promoters from Opium Poppy Are Conserved in Transgenic Tobacco1

Opium poppy (Papaver somniferum) contains a large family of tyrosine/dihydroxyphenylalanine decarboxylase (tydc) genes involved in the biosynthesis of benzylisoquinoline alkaloids and cell wall-bound hydroxycinnamic acid amides. Eight members from two distinct gene subfamilies have been isolated, tydc1, tydc4, tydc6, tydc8, and tydc9 in one group and tydc2, tydc3, and tydc7 in the other. The tydc8 and tydc9 genes were located 3.2 kb apart on one genomic clone, suggesting that the family is clustered. Transcripts for most tydc genes were detected only in roots. Only tydc2 and tydc7 revealed expression in both roots and shoots, and TYDC3 mRNAs were the only specific transcripts detected in seedlings. TYDC1, TYDC8, and TYDC9 mRNAs, which occurred in roots, were not detected in elicitor-treated opium poppy cultures. Expression of tydc4, which contains a premature termination codon, was not detected under any conditions. Five tydc promoters were fused to the β-glucuronidase (GUS) reporter gene in a binary vector. All constructs produced transient GUS activity in microprojectile-bombarded opium poppy and tobacco (Nicotiana tabacum) cell cultures. The organ- and tissue-specific expression pattern of tydc promoter-GUS fusions in transgenic tobacco was generally parallel to that of corresponding tydc genes in opium poppy. GUS expression was most abundant in the internal phloem of shoot organs and in the stele of roots. Select tydc promoter-GUS fusions were also wound induced in transgenic tobacco, suggesting that the basic mechanisms of developmental and inducible tydc regulation are conserved across plant species.

Family of MADS-Box Genes Expressed Early in Male and Female Reproductive Structures of Monterey Pine

Three MADS-box genes isolated from Monterey pine (Pinus radiata), PrMADS1, PrMADS2, and PrMADS3, are orthologs to members of the AGL2 and AGL6 gene subfamilies in Arabidopsis. These genes were expressed during early stages of pine shoot development in differentiating seed- and pollen-cone buds. Their transcripts were found within a group of cells that formed ovuliferous scale and microsporophyll primordia. Expression of PrMADS3 was also detected in a group of cells giving rise to needle primordia within differentiated vegetative buds, and in needle primordia.

The Three-Dimensional Structure of Aspergillus niger Pectin Lyase B at 1.7-Å Resolution1

The three-dimensional structure of Aspergillus niger pectin lyase B (PLB) has been determined by crystallographic techniques at a resolution of 1.7 Å. The model, with all 359 amino acids and 339 water molecules, refines to a final crystallographic R factor of 16.5%. The polypeptide backbone folds into a large right-handed cylinder, termed a parallel β helix. Loops of various sizes and conformations protrude from the central helix and probably confer function. The largest loop of 53 residues folds into a small domain consisting of three antiparallel β strands, one turn of an α helix, and one turn of a 310 helix. By comparison with the structure of Erwinia chrysanthemi pectate lyase C (PelC), the primary sequence alignment between the pectate and pectin lyase subfamilies has been corrected and the active site region for the pectin lyases deduced. The substrate-binding site in PLB is considerably less hydrophilic than the comparable PelC region and consists of an extensive network of highly conserved Trp and His residues. The PLB structure provides an atomic explanation for the lack of a catalytic requirement for Ca2+ in the pectin lyase family, in contrast to that found in the pectate lyase enzymes. Surprisingly, however, the PLB site analogous to the Ca2+ site in PelC is filled with a positive charge provided by a conserved Arg in the pectin lyases. The significance of the finding with regard to the enzymatic mechanism is discussed.

Isolation of a superfamily of candidate disease-resistance genes in soybean based on a conserved nucleotide-binding site.

The tobacco N and Arabidopsis RPS2 genes, among several recently cloned disease-resistance genes, share highly conserved structure, a nucleotide-binding site (NBS). Using degenerate oligonucleotide primers for the NBS region of N and RPS2, we have amplified and cloned the NBS sequences from soybean. Each of these PCR-derived NBS clones detected low-or moderate-copy soybean DNA sequences and belongs to 1 of 11 different classes. Sequence analysis showed that all PCR clones encode three motifs (P-loop, kinase-2, and kinase-3a) of NBS nearly identical to those in N and RPS2. The intervening region between P-loop and kinase-3a of the 11 classes has high (26% average) amino acid sequence similarity to the N gene although not as high (19% average) to RPS2. These 11 classes represent a superfamily of NBS-containing soybean genes that are homologous to N and RPS2. Each class or subfamily was assessed for its positional association with known soybean disease-resistance genes through near-isogenic line assays, followed by linkage analysis in F2 populations using restriction fragment length polymorphisms. Five of the 11 subfamilies have thus far been mapped to the vicinity of known soybean genes for resistance to potyviruses (Rsv1 and Rpv), Phytophthora root rot (Rps1, Rps2, and Rps3), and powdery mildew (rmd). The conserved N- or RPS2-homologous NBS sequences and their positional associations with mapped soybean-resistance genes suggest that a number of the soybean disease-resistance genes may belong to this superfamily. The candidate subfamilies of NBS-containing genes identified by genetic mapping should greatly facilitate the molecular cloning of disease-resistance genes.

Organization and expression of canine olfactory receptor genes.

Four members of the canine olfactory receptor gene family were characterized. The predicted proteins shared 40-64% identity with previously identified olfactory receptors. The four subfamilies identified in Southern hybridization experiments had as few as 2 and as many as 20 members. All four genes were expressed exclusively in olfactory epithelium. Expression of multiple members of the larger subfamilies was detected, suggesting that most if not all of the cross-hybridizing bands in genomic Southern blots represented actively transcribed olfactory receptor genes. Analysis of large DNA fragments using Southern blots of pulsed-field gels indicated that subfamily members were clustered together, and that two of the subfamilies were closely linked in the dog genome. Analysis of the four olfactory receptor gene subfamilies in 26 breeds of dog provided evidence that the number of genes per subfamily was stable in spite of differential selection on the basis of olfactory acuity in scent hounds, sight hounds, and toy breeds.

The primary structures of the Archaeon Halobacterium salinarium blue light receptor sensory rhodopsin II and its transducer, a methyl-accepting protein.

Recently, a large family of transducer proteins in the Archaeon Halobacterium salinarium was identified. On the basis of the comparison of the predicted structural domains of these transducers, three distinct subfamilies of transducers were proposed. Here we report isolation, complete gene sequences, and analysis of the encoded primary structures of transducer gene htrII, a member of family B, and its blue light receptor gene (sopII) of sensory rhodopsin II (SRII). The start codon ATG of the 714-bp sopII gene is one nucleotide beyond the termination codon TGA of the 2298-bp htrII gene. The deduced protein sequence of HtrII predicts a eubacterial chemotaxis transducer type with two hydrophobic membrane-spanning segments connecting sizable domains in the periplasm and cytoplasm. HtrII has a common feature with HtrI, the sensory rhodopsin I transducer; like HtrI, HtrII possesses a hydrophilic loop structure just after the second transmembrane segment. The C-terminal 299 residues (765 amino acid residues total) of HtrII show strong homology to the signaling and methylation domain of eubacterial transducer Tsr. The hydropathy plot of the primary structure of SRII indicates seven membrane-spanning alpha-helical segments, a characteristic feature of retinylidene proteins ("rhodopsins") from a widespread family of photoactive pigments. SRII shows high identity with SRI (42%), bacteriorhodopsin (BR) (32%), and halorhodopsin (24%). The crucial positions for retinal binding sites in these proteins are nearly identical, with the exception of Met-118 (numbering according to the mature BR sequence), which is replaced by Val in SRII. In BR, residues Asp-85 and Asp-96 are crucial in proton pumping. In SRII, the position corresponding to Asp-85 in BR is conserved, but the corresponding position of Asp-96 is replaced by an aromatic Tyr. Coexpression of the htrII and sopII genes restores SRII phototaxis to a mutant (Pho81) that contains a deletion in the htrI/sopI and insertion in htrII/sopII regions. This paper describes the first example that both HtrI and HtrII exist in the same halobacterial cell, confirming that different sensory rhodopsins SRI and SRII in the same organism have their own distinct transducers.

Tiggers and DNA transposon fossils in the human genome.

We report several classes of human interspersed repeats that resemble fossils of DNA transposons, elements that move by excision and reintegration in the genome, whereas previously characterized mammalian repeats all appear to have accumulated by retrotransposition, which involves an RNA intermediate. The human genome contains at least 14 families and > 100,000 degenerate copies of short (180-1200 bp) elements that have 14- to 25-bp terminal inverted repeats and are flanked by either 8 bp or TA target site duplications. We describe two ancient 2.5-kb elements with coding capacity, Tigger1 and -2, that closely resemble pogo, a DNA transposon in Drosophila, and probably were responsible for the distribution of some of the short elements. The deduced pogo and Tigger proteins are related to products of five DNA transposons found in fungi and nematodes, and more distantly, to the Tc1 and mariner transposases. They also are very similar to the major mammalian centromere protein CENP-B, suggesting that this may have a transposase origin. We further identified relatively low-copy-number mariner elements in both human and sheep DNA. These belong to two subfamilies previously identified in insect genomes, suggesting lateral transfer between diverse species.

Fate of biological control introductions: monitoring an Australian fungal pathogen of grasshoppers in North America.

In North America there are two generally recognized pathotypes (pathotypes 1 and 2) of the fungus Entomophaga grylli which show host-preferential infection of grasshopper subfamilies. Pathotype 3, discovered in Australia, has a broader grasshopper host range and was considered to be a good biocontrol agent. Between 1989 and 1991 pathotype 3 was introduced at two field sites in North Dakota. Since resting spores are morphologically indistinguishable among pathotypes, we used pathotype-specific DNA probes to confirm pathotype identification in E. grylli-infected grasshoppers collected at the release sites in 1992, 1993, and 1994. In 1992, up to 23% of E. grylli-infected grasshoppers of the subfamilies Melanoplinae, Oedipodinae, and Gomphocerinae were infected by pathotype 3, with no infections > 1 km from the release sites. In 1993, pathotype 3 infections declined to 1.7%. In 1994 grasshopper populations were low and no pathotype 3 infections were found. The frequency of pathotype 3 infection has declined to levels where its long-term survival in North America is questionable. Analyses of biocontrol releases are critical to evaluating the environmental risks associated with these ecological manipulations, and molecular probes are powerful tools for monitoring biocontrol releases.

The molecular role of the common gamma c subunit in signal transduction reveals functional asymmetry within multimeric cytokine receptor complexes.

The specific signal transduction function of the gamma c subunit in the interleukin (IL) 2, IL-4, IL-7, IL-9, and IL-15 receptor complexes remains undefined. The present structure-function analyses demonstrated that the entire cytoplasmic tail of gamma c could be functionally replaced in the IL-2 receptor (IL-2R) signaling complex by a severely truncated erythropoietin receptor cytoplasmic domain lacking tyrosine residues. Heterodimerization of IL-2R beta with either gamma c or the truncated erythropoietin receptor chain led to an array of specific signals normally derived from the native IL-2R despite the substitution of Janus kinase JAK2 for JAK3 in the receptor complex. These findings thus suggest a model in which the gamma c subunit serves as a common and generic "trigger" chain by providing a nonspecific Janus kinase for signaling program initiation, while signal specificity is determined by the unique "driver" subunit in each of the gamma c- containing receptor complexes. Furthermore, these results may have important functional implications for the asymmetric design of many cytokine receptor complexes and the evolutionary design of receptor subfamilies that share common trigger or driver subunits.

Genetics and the origin of European languages.

A new set of European genetic data has been analyzed to dissect independent patterns of geographic variation. The most important cause of European genetic variation has been confirmed to correspond to the migration of Neolithic farmers from the area of origin of agriculture in the Middle East. The next most important component of genetic variation is apparently associated with a north-south gradient possibly due to adaptation to cold climates but also to the differentiation of the Uralic and the Indo-European language-speaking people; however, the relevant correlations are not significantly different from zero after elimination of the spatial autocorrelation. The third component is highly correlated with the infiltration of the Yamna ("Kurgan") people, nomadic pastoralists who domesticated the horse and who have been claimed to have spread Indo-European languages to Europe; this association, which is statistically significant even when taking spatial autocorrelations into account, does not completely exclude the hypothesis of Indo-European as the language of Neolithic farmers. It is possible that both expansions were responsible for the spread of different subfamilies of Indo-European languages, but our genetic data cannot resolve their relative importance.

Morfologia, taxonomia e análise cladística da tribo Kingsleyini (Crustacea: Decapoda: Pseudothelphusidae)

Kingsleyini corresponde a uma das cinco tribos de Pseudothelphusidae, grupo exclusivamente americano de caranguejos de água doce. Atualmente a tribo inclui 59 espécies agrupadas em 13 gêneros, com distribuição associada aos rios, riachos e igarapés das bacias do Amazonas e do Orinoco. Desde a criação de Kingsleyini o aumento de novos táxons atribuídos a esta tribo não tem sido acompanhado por estudos cladísticos. No presente trabalho é realizada a análise cladística de Kingsleyini, acompanhada de uma revisão morfológica e taxonômica do grupo. Com este propósito, foram estudados espécimes de 60 espécies representantes das cinco tribos e duas subfamílias inclusas em Pseudothelphusidae. O material estudado se encontra depositado nas coleções carcinológicas de seis instituições e inclui os tipos nominais de 29 espécies. Na revisão morfológica foram descritas e ilustradas estruturas somáticas e sexuais da morfologia externa do grupo de estudo. Os estudos morfológicos foram auxiliados por técnicas de Microscopia Electrônica de Varredura (MEV) e cortes histológicos. A partir destas observações foram propostas modificações na terminologia utilizada para denominar as estruturas do primeiro apêndice sexual masculino (primeiro gonópodo) em Kingsleyini. A parte taxonômica deste trabalho inclui chaves de identificação, mapas de distribuição, listas sinonímicas e a descrição e ilustração do primeiro apêndice sexual masculino para a grande maioria das espécies examinadas, assim como a diagnose dos gêneros considerados monofiléticos. A análise filogenética foi realizada a partir de 92 caracteres obtidos de 57 táxons terminais: 49 terminais do grupo interno (Kingsleyini) e oito do grupo-externo (representantes dos demais Pseudothelphusidae). Como resultado da análise cladística foram obtidas seis hipóteses filogenéticas igualmente parcimoniosas: todas elas apoiam o monofiletismo de Kingsleyini e a exclusão do gênero Spirocarcinus da tribo. O monofiletismo dos gêneros Fredius, Kingsleya, Eudaniela e Rodriguezus também encontra-se sustentado em todas as hipóteses filogenéticas obtidas, enquanto que os gêneros Microthelphusa, Neopseudothelphusa, Orthothelphusa e Brasiliothelphusa revelaram-se parafiléticos.

Beta diversity at multiple hierarchical levels: explaining the high diversity of scarab beetles in tropical montane forests

Aim: High gamma diversity in tropical montane forests may be ascribed to high geographical turnover of community composition, resulting from population isolation that leads to speciation. We studied the evolutionary processes responsible for diversity and turnover in assemblages of tropical scarab beetles (Scarabaeidae) by assessing DNA sequence variation at multiple hierarchical levels. Location: A 300-km transect across six montane forests (900–1100 m) in Costa Rica. Methods: Assemblages of Scarabaeidae (subfamilies Dynastinae, Rutelinae, Melolonthinae) including 118 morphospecies and > 500 individuals were sequenced for the cox1 gene to establish species limits with a mixed Yule–coalescent method. A species-level phylogenetic tree was constructed from cox1 and rrnL genes. Total diversity and turnover among assemblages were then assessed at three hierarchical levels: haplotypes, species and higher clades. Results: DNA-based analyses showed high turnover among communities at all hierarchical levels. Turnover was highest at the haplotype level (community similarity 0.02–0.12) and decreased with each step of the hierarchy (species: 0.21–0.46; clades: 0.41–0.43). Both compositional and phylogenetic similarities of communities were geographically structured, but turnover was not correlated with distance among forests. When three major clades were investigated separately, communities of Dynastinae showed consistently higher alpha diversity, larger species ranges and lower turnover than Rutelinae and Melolonthinae. Main conclusions: Scarab communities of montane forests show evidence of evolutionary persistence of communities in relative isolation, presumably tracking suitable habitats elevationally to accommodate climatic changes. Patterns of diversity on all hierarchical levels seem to be determined by restricted dispersal, and differences in Dynastinae could be explained by their greater dispersal ability. Community-wide DNA sequencing across multiple lineages and hierarchical levels reveals the evolutionary processes that led to high beta diversity in tropical montane forests through time.

The Plasmodium serine-type SERA proteases display distinct expression patterns and non-essential in vivo roles during life cycle progression of the malaria parasite

Parasite proteases play key roles in several fundamental steps of the Plasmodium life cycle, including haemoglobin degradation, host cell invasion and parasite egress. Plasmodium exit from infected host cells appears to be mediated by a class of papain-like cysteine proteases called 'serine repeat antigens' (SERAs). A SERA subfamily, represented by Plasmodium falciparum SERA5, contains an atypical active site serine residue instead of a catalytic cysteine. Members of this SERAser subfamily are abundantly expressed in asexual blood stages, rendering them attractive drug and vaccine targets. In this study, we show by antibody localization and in vivo fluorescent tagging with the red fluorescent protein mCherry that the two P. berghei serine-type family members, PbSERA1 and PbSERA2, display differential expression towards the final stages of merozoite formation. Via targeted gene replacement, we generated single and double gene knockouts of the P. berghei SERAser genes. These loss-of-function lines progressed normally through the parasite life cycle, suggesting a specialized, non-vital role for serine-type SERAs in vivo. Parasites lacking PbSERAser showed increased expression of the cysteine-type PbSERA3. Compensatory mechanisms between distinct SERA subfamilies may thus explain the absence of phenotypical defect in SERAser disruptants, and challenge the suitability to develop potent antimalarial drugs based on specific inhibitors of Plasmodium serine-type SERAs.