Characterization of the porcine AMPK alpha 2 catalytic subunitgene (PRKAA2): genomic structure, polymorphism detection and association study.

AMP-activated protein kinase (AMPK), known as a key regulator of cellular energy homeostasis, plays an important role in regulation of glucose and lipid metabolism, and protein synthesis in mammals. The characterization of porcine PRKAA2 encoding the alpha 2 catalytic subunit of AMPK is reported in this study. PRKAA2 was assigned to porcine chromosome 6q by analysis of radiation hybrids (IMpRH panel), and its genomic structure was determined by BAC sequencing. PRKAA2 spans more than 62 kb and consists of nine exons and eight introns. A total of 25 polymorphisms were identified by re-sequencing approximately 7 kb, including all the exons, exon-intron boundaries and 5' and 3' gene flanking regions using twelve founder animals of a Mangalitsa x Piétrain intercross. Neither of two single nucleotide polymorphisms (SNPs) found in the coding region caused an amino acid substitution. Two SNPs (NM_214266.1: c.236+142A>G and NM_214266.1: c.630C>T) in PRKAA2 were genotyped in the Mangalitsa x Piétrain F(2) cross (n = 589) and two commercial populations [Piétrain (n = 1173) and German Landrace (n = 536)] and evaluated for association with traits of interest (muscle development and fat deposition). Single SNP and haplotype analyses revealed weak associations between the PRKAA2 genotypes and loin muscle area in the investigated populations.

Genome sequence, comparative analysis, and population genetics of the domestic horse.

We report a high-quality draft sequence of the genome of the horse (Equus caballus). The genome is relatively repetitive but has little segmental duplication. Chromosomes appear to have undergone few historical rearrangements: 53% of equine chromosomes show conserved synteny to a single human chromosome. Equine chromosome 11 is shown to have an evolutionary new centromere devoid of centromeric satellite DNA, suggesting that centromeric function may arise before satellite repeat accumulation. Linkage disequilibrium, showing the influences of early domestication of large herds of female horses, is intermediate in length between dog and human, and there is long-range haplotype sharing among breeds.

Interspecies comparison of gene structure and computational analysis of gene regulation of 17beta-hydroxysteroid dehydrogenase type 1.

17Beta-hydroxysteroid dehydrogenase type 1 (HSD17B1) is a key enzyme of 17beta-estradiol biosynthesis, and in rodents is additionally involved in testosterone biosynthesis. The human HSD17B1 gene, located on chromosome 17q12-21, is duplicated in tandem, with the 3'-copy being the functional gene. Here we show by sequencing the gene from a diverse set of related species that this duplication is of very recent evolutionary origin, having occurred in the common ancestor of Hominoidae (apes and humans) while being absent in the closely related Old World monkeys (Macaca) and the outgroup species Tupaia belangeri and Mus musculus. By computational analysis of the conserved regulatory elements in the 5'-untranslated (5'-UTR) and putative promoter region of the HSD17B1 gene and, where present, pseudogene, across our broad sample of species we can show significant differences that might point to the origin of the divergent substrate specificity of human and rodent HSD17B1 and highlight potential functionally relevant differences in regulatory patterns in different evolutionary lineages.

DNA sequence and comparative analysis of chimpanzee chromosome 22.

Human-chimpanzee comparative genome research is essential for narrowing down genetic changes involved in the acquisition of unique human features, such as highly developed cognitive functions, bipedalism or the use of complex language. Here, we report the high-quality DNA sequence of 33.3 megabases of chimpanzee chromosome 22. By comparing the whole sequence with the human counterpart, chromosome 21, we found that 1.44% of the chromosome consists of single-base substitutions in addition to nearly 68,000 insertions or deletions. These differences are sufficient to generate changes in most of the proteins. Indeed, 83% of the 231 coding sequences, including functionally important genes, show differences at the amino acid sequence level. Furthermore, we demonstrate different expansion of particular subfamilies of retrotransposons between the lineages, suggesting different impacts of retrotranspositions on human and chimpanzee evolution. The genomic changes after speciation and their biological consequences seem more complex than originally hypothesized.

Toward a catalog of human genes and proteins: sequencing and analysis of 500 novel complete protein coding human cDNAs.

With the complete human genomic sequence being unraveled, the focus will shift to gene identification and to the functional analysis of gene products. The generation of a set of cDNAs, both sequences and physical clones, which contains the complete and noninterrupted protein coding regions of all human genes will provide the indispensable tools for the systematic and comprehensive analysis of protein function to eventually understand the molecular basis of man. Here we report the sequencing and analysis of 500 novel human cDNAs containing the complete protein coding frame. Assignment to functional categories was possible for 52% (259) of the encoded proteins, the remaining fraction having no similarities with known proteins. By aligning the cDNA sequences with the sequences of the finished chromosomes 21 and 22 we identified a number of genes that either had been completely missed in the analysis of the genomic sequences or had been wrongly predicted. Three of these genes appear to be present in several copies. We conclude that full-length cDNA sequencing continues to be crucial also for the accurate identification of genes. The set of 500 novel cDNAs, and another 1000 full-coding cDNAs of known transcripts we have identified, adds up to cDNA representations covering 2%--5 % of all human genes. We thus substantially contribute to the generation of a gene catalog, consisting of both full-coding cDNA sequences and clones, which should be made freely available and will become an invaluable tool for detailed functional studies.

Novel features in a combined polyketide synthase/non-ribosomal peptide synthetase: the myxalamid biosynthetic gene cluster of the myxobacterium Stigmatella aurantiaca Sga15.

BACKGROUND: Myxobacteria have been well established as a potent source for natural products with biological activity. They produce a considerable variety of compounds which represent typical polyketide structures with incorporated amino acids (e.g. the epothilons, the myxothiazols and the myxalamids). Several of these secondary metabolites are effective inhibitors of the electron transport via the respiratory chain and have been widely used. Molecular cloning and characterization of the genes governing the biosynthesis of these structures is of considerable interest, because such information adds to the limited knowledge as to how polyketide synthases (PKSs) and non-ribosomal peptide synthetases (NRPSs) interact and how they might be manipulated in order to form novel antibiotics. RESULTS: A DNA region of approximately 50000 base pairs from Stigmatella aurantiaca Sga15 was sequenced and shown by gene disruption to be involved in myxalamid biosynthesis. Sequence analysis reveals that the myxalamids are formed by a combined PKS/NRPS system. The terminal NRPS MxaA extends the assembled polyketide chain of the myxalamids with alanine. MxaA contains an N-terminal domain with homology to NAD binding proteins, which is responsible during the biogenesis for a novel type of reductive chain release giving rise to the 2-amino-propanol moiety of the myxalamids. The last module of the PKS reveals an unprecedented genetic organization; it is encoded on two genes (mxaB1 and mxaB2), subdividing the domains of one module from each other. A sequence comparison of myxobacterial acyl-transferase domains with known systems from streptomycetes and bacilli reveals that consensus sequences proposed to be specific for methylmalonyl-CoA and malonyl-CoA are not always reliable. CONCLUSIONS: The complete biosynthetic gene cluster of the myxalamid-type electron transport inhibitor from S. aurantiaca Sga15 has been cloned and analyzed. It represents one of the few examples of combined PKS/NRPS systems, the analysis and manipulation of which has the potential to generate novel hybrid structures via combinatorial biosynthesis (e.g. via module-swapping techniques). Additionally, a new type of reductive release from PKS/NRPS systems is described.

Sequence and analysis of chromosome 3 of the plant Arabidopsis thaliana.

Arabidopsis thaliana is an important model system for plant biologists. In 1996 an international collaboration (the Arabidopsis Genome Initiative) was formed to sequence the whole genome of Arabidopsis and in 1999 the sequence of the first two chromosomes was reported. The sequence of the last three chromosomes and an analysis of the whole genome are reported in this issue. Here we present the sequence of chromosome 3, organized into four sequence segments (contigs). The two largest (13.5 and 9.2 Mb) correspond to the top (long) and the bottom (short) arms of chromosome 3, and the two small contigs are located in the genetically defined centromere. This chromosome encodes 5,220 of the roughly 25,500 predicted protein-coding genes in the genome. About 20% of the predicted proteins have significant homology to proteins in eukaryotic genomes for which the complete sequence is available, pointing to important conserved cellular functions among eukaryotes.

The myxochelin iron transport regulon of the myxobacterium Stigmatella aurantiaca Sg a15.

The biosynthetic gene cluster of the myxochelin-type iron chelator was cloned from Stigmatella aurantiaca Sg a15 and characterized. This catecholate siderophore was only known from two other myxobacteria. The biosynthetic genes of 2,3-dihydroxybenzoic acid are located in the cluster (mxcC-mxcF). Two molecules of 2, 3-dihydroxybenzoic acid are activated and condensed with lysine in a unique way by a protein homologous to nonribosomal peptide synthetases (MxcG). Inactivation of mxcG, which encodes an adenylation domain for lysine, results in a myxochelin negative mutant unable to grow under iron-limiting conditions. Growth could be restored by adding Fe3+, myxochelin A or B to the medium. Inactivation of mxcD leads to the same phenotype. A new type of reductive release from nonribosomal peptide synthetases of the 2, 3-dihydroxybenzoic acid bis-amide of lysine from MxcG, catalyzed by a protein domain with homology to NAD(P) binding sites, is discussed. The product of a gene, encoding a protein similar to glutamate-1-semialdehyde 2,1-aminomutases (mxcL), is assumed to transaminate the aldehyde that is proposed as an intermediate. Further genes encoding proteins homologous to typical iron utilization and iron uptake polypeptides are reported.

The DNA sequence of human chromosome 21.

Chromosome 21 is the smallest human autosome. An extra copy of chromosome 21 causes Down syndrome, the most frequent genetic cause of significant mental retardation, which affects up to 1 in 700 live births. Several anonymous loci for monogenic disorders and predispositions for common complex disorders have also been mapped to this chromosome, and loss of heterozygosity has been observed in regions associated with solid tumours. Here we report the sequence and gene catalogue of the long arm of chromosome 21. We have sequenced 33,546,361 base pairs (bp) of DNA with very high accuracy, the largest contig being 25,491,867 bp. Only three small clone gaps and seven sequencing gaps remain, comprising about 100 kilobases. Thus, we achieved 99.7% coverage of 21q. We also sequenced 281,116 bp from the short arm. The structural features identified include duplications that are probably involved in chromosomal abnormalities and repeat structures in the telomeric and pericentromeric regions. Analysis of the chromosome revealed 127 known genes, 98 predicted genes and 59 pseudogenes.

Identification of an L-dopa decarboxylase gene from Sorangium cellulosum So ce90.

During a screening program intended to identify genes encoding enzymes typical for secondary metabolism in Sorangium cellulosum So ce90, an aromatic amino acid decarboxylase gene (ddc) was detected. Expression of ddc in Escherichia coli and subsequent enzyme assays with cell-free extracts confirmed the proposed function derived from amino acid sequence comparisons. In contrast to other aromatic amino acid decarboxylases of eukaryotic origin, the S. cellulosum Ddc converted only L-dihydroxy phenylalanine. This is the first report of a gene encoding an L-dihydroxy phenylalanine decarboxylase in bacteria.

New lessons for combinatorial biosynthesis from myxobacteria. The myxothiazol biosynthetic gene cluster of Stigmatella aurantiaca DW4/3-1.

The biosynthetic mta gene cluster responsible for myxothiazol formation from the fruiting body forming myxobacterium Stigmatella aurantiaca DW4/3-1 was sequenced and analyzed. Myxothiazol, an inhibitor of the electron transport via the bc(1)-complex of the respiratory chain, is biosynthesized by a unique combination of several polyketide synthases (PKS) and nonribosomal peptide synthetases (NRPS), which are activated by the 4'-phosphopantetheinyl transferase MtaA. Genomic replacement of a fragment of mtaB and insertion of a kanamycin resistance gene into mtaA both impaired myxothiazol synthesis. Genes mtaC and mtaD encode the enzymes for bis-thiazol(ine) formation and chain extension on one pure NRPS (MtaC) and on a unique combination of PKS and NRPS (MtaD). The genes mtaE and mtaF encode PKSs including peptide fragments with homology to methyltransferases. These methyltransferase modules are assumed to be necessary for the formation of the proposed methoxy- and beta-methoxy-acrylate intermediates of myxothiazol biosynthesis. The last gene of the cluster, mtaG, again resembles a NRPS and provides insight into the mechanism of the formation of the terminal amide of myxothiazol. The carbon backbone of an amino acid added to the myxothiazol-acid is assumed to be removed via an unprecedented module with homology to monooxygenases within MtaG.

The complete nucleotide sequence of the egg drop syndrome virus: an intermediate between mastadenoviruses and aviadenoviruses.

The complete nucleotide sequence of an avian adenovirus, the egg drop syndrome (EDS) virus, was determined. The total genome length is 33,213 nucleotides, resulting in a molecular weight of 21.9 x 10(6). The GC content is only 42.5%. Between map units 3.5 and 76.9, the distribution of open reading frames with homology to known genes is similar to that reported for other mammalian and avian adenoviruses. However, no homologies to adenovirus genes such as E1A, pIX, pV, and E3 could be found. Outside this region, several open reading frames were identified without any obvious homology to known adenovirus proteins. In the region organized similarly as other adenoviral genomes, most homologies were found to an ovine adenovirus (OAV strain 287). The highest level of amino acid identity was found for the hexon proteins of EDS and OAV. The virus-associated RNA (VA RNA) was identified thanks to the homology with the VA RNA of fowl adenovirus serotype 1 (FAV1). Similarities with FAV1 were also found in the fiber protein. Our results demonstrate that the avian EDS virus represents an intermediate between mammalian and avian adenoviruses. The nucleotide sequence and genomic organization of the EDS virus reflect the heterogeneity of the aviadenovirus genus and the Adenoviridae family.

The nucleotide sequence of Saccharomyces cerevisiae chromosome IV.

The complete DNA sequence of the yeast Saccharomyces cerevisiae chromosome IV has been determined. Apart from chromosome XII, which contains the 1-2 Mb rDNA cluster, chromosome IV is the longest S. cerevisiae chromosome. It was split into three parts, which were sequenced by a consortium from the European Community, the Sanger Centre, and groups from St Louis and Stanford in the United States. The sequence of 1,531,974 base pairs contains 796 predicted or known genes, 318 (39.9%) of which have been previously identified. Of the 478 new genes, 225 (28.3%) are homologous to previously identified genes and 253 (32%) have unknown functions or correspond to spurious open reading frames (ORFs). On average there is one gene approximately every two kilobases. Superimposed on alternating regional variations in G+C composition, there is a large central domain with a lower G+C content that contains all the yeast transposon (Ty) elements and most of the tRNA genes. Chromosome IV shares with chromosomes II, V, XII, XIII and XV some long clustered duplications which partly explain its origin.

Grafting of discontinuous sites: a protein modeling strategy.

A strategy for modeling continuous as well as discontinuous sites in protein structures has been developed. Central to this modeling strategy is the search algorithm of FITSITE, a program to search a given target structure for suitable combinations of backbone positions mirroring as closely as possible the geometric relationships of a source structural motif of interest. All target sites detected by FITSITE are further refined to mimic the source geometry. The sidechain rotamer library concept fails to precisely describe side chains involved in coordinative bonding (e.g, metal binding sites). Therefore an algorithm using detailed data-base bonding parameter information was applied for the side-chain construction. The FITSITE program and the subsequent processing of the program output are presented in a test case. The Rop protein, a four-helix bundle structure, served as the target protein. It was searched for candidate sites to model a variety of metal binding sites, with structures extracted from Brookhaven Protein Database entries. The preliminary protein models were investigated for structural overlaps with neighboring residues by interactive computer graphics; if required, additional changes were performed. A set of parameters for energy minimization with AMBER (including metal ions) was developed, and the completed Rop variants were energy minimized. Finally, 12 potentially metal binding Rop variants were selected for production via genetic engineering.

The 'shortmer' approach to nucleic acid sequence analysis. I: Computer simulation of sequencing projects to find economical primer sets.

In principle it is most economical to sequence large DNA fragments consecutively ('primer walking'), provided there is an immediate supply of sequencing primers. To solve the problem of primer supply we previously suggested generating a bank of short oligonucleotide primers ('shortmers'). In every sequencing reaction shortmers would have to be selected from this bank that are suitable to hybridize adjacently on the sequencing template. After their ligation the shortmers would form a long, and hence more specific, primer in the subsequent sequencing reaction. In the present study a computer simulation of large sequencing projects revealed a reduced set of approximately 12,000 selected octanucleotides (out of all 65,536) retaining maximum priming flexibility and minimum redundant information on the simulated sequence analyses. Establishing routine protocols for nucleic acid sequencing following the shortmer approach will abolish the tightest bottleneck of the consecutive sequencing route (primer supply) and hence may render this general scheme more attractive than the shotgun sequencing scheme. A twofold (or more) speed-up of genome sequencing projects by the shortmer approach may be assumed.

Phospholipid-hydroperoxide glutathione peroxidase. Genomic DNA, cDNA, and deduced amino acid sequence.

The complete amino acid sequence of the selenoprotein phospholipid-hydroperoxide glutathione peroxidase (PHGPX) from pig heart has been deduced from the corresponding genomic DNA, the cDNA covering the coding region, and by sequencing the N terminus of the protein. The maximum length of the peptide chain derived from the cDNA amounts to 170 amino acid residues. By protein sequencing the N-terminal residues methionine and cysteine of the deduced sequence were found to be cleaved. The molecular mass of 19,671 Da obtained by laser desorption mass spectroscopy, however, significantly exceeds the mean molecular mass of 19,257.09 calculated for the sequence 3-170 of PHGPX, thus indicating posttranscriptional modification. In contrast to glutathione peroxidase (GPX) the coding area of the PHGPX gene is composed of seven exons. Only the amino acid sequences encoded by the third and fifth exon are highly homologous to GPX sequences. The amino acid residues selenocysteine, tryptophan, and glutamine forming the catalytic site in bovine GPX are conserved in homologous positions of PHGPX, whereas the arginine residues presumed to bind GSH in GPX are not. Gaps in the PHGPX sequence correspond to subunit interaction sites of the tetrameric GPX. The data suggest an identical catalytic mechanism of the selenoperoxidases, a less stringent substrate specificity of PHGPX, and explain the monomeric nature of PHGPX. As in other selenoproteins, the selenocysteine residue of PHGPX is encoded by UGA. The 3'-untranslated region (UTR) of the PHGPX shows a limited consensus with that of GPX and 5'-deiodinase, where it was shown to be responsible for the decoding of UGA as selenocysteine. The 3'-UTR of PHGPX can form a stem/loop as in other mammalian selenoprotein genes. The 5'-UTR and the first intron of the PHGPX gene contain a variety of putative regulatory elements indicating hormonal control.

A phylogeny of New World Deltocephalus-like leafhopper genera based on mitochondrial 16S ribosomal DNA sequences.

The Deltocephalus-like leafhoppers are a group of grass-feeding insects that share common forms in the male genitalia, chiefly a fused aedeagal-connective structure. The group contains 27 New World genera. The phylogenetic relationships among these genera and other genera in the tribe are unclear because fusion of the aedeagus and connective, although clearly apomorphic, has occurred in both related and unrelated groups. An independent set of characters is required to test the hypothesis that the Deltocephalus-like genera are monophyletic and to derive relationships among these genera. We examined 562 sites in the 3'-end of the mitochondrial 16S ribosomal DNA from 21 species representing 19 New World genera. Eight species in closely related groups or tribes were also analyzed. The region was sequenced directly from PCR-amplified DNA of field-collected and museum-preserved dried specimens. The transversion rate was twice as high as the transition rate; however, 74% of these were adenine-thymine transversions. Analysis of secondary structure indicated that substitution rates were equal in stem and loop forming regions of the processed rRNA. Using either parsimony or distance analysis, derived phylogenies suggested that, with the exception of the genus Cabrulus, the Deltocephalus-like leafhoppers are genera within a monophyletic group.

Analysis of cytokine mRNA levels in interleukin-4-transgenic mice by quantitative polymerase chain reaction.

Interleukin (IL)-4-transgenic mice were used as a model system to study the consequences of low levels of IL-4 expression for the expression of other cytokines examined by quantitative polymerase chain reaction (PCR). For this purpose, a plasmid was constructed which contains, in tandem array, 5' and 3' primer sequences specific for the cytokine genes IL-1 to IL-6, tumor necrosis factor (TNF), lymphotoxin (LT), interferon (IFN)-gamma and beta-actin. During co-amplification, target and control DNA compete for the primers and the amount of PCR product is proportional to the amount of input DNA. Competitive PCR was performed first to adjust the cDNA to be compared to identical concentrations of beta-actin cDNA and subsequently to determine cytokine mRNA levels from spleen cells of normal and IL-4-transgenic animals. The sensitivity of this approach was demonstrated by the capability to detect a twofold difference in IL-4 mRNA levels between IL-4-transgenic heterozygous and homozygous animals. Upon lipopolysaccharide activation, the IL-4 transgene which is expressed essentially in B lymphocytes was induced approximately 50-fold. Several cytokine mRNA such as those coding for IL-5, IL-6, IFN-gamma and also the IL-4 receptor were found to be up-regulated in IL-4-transgenic mice, whereas IL-1, IL-2, IL-3, TNF and LT mRNA levels did not seemed to be influenced by IL-4. A possible functional significance of the elevated IFN-gamma mRNA was demonstrated by showing that (a) CD23 expression was not increased, and (b) Mac-1+ cells were markedly increased in the spleen of transgenic mice.

Effect of a single 3'-methylene phosphonate linkage on the conformation of an A-DNA octamer double helix.

The three-dimensional structure of the self-complementary DNA octamer d(GCCCGpGGC) has been determined in the crystalline state using X-ray diffraction data to a nominal resolutoin of 2.12 measured from a very small crystal at DESY, Hamburg. The structure was refined with stereochemical restraints to an R value of 17.1%. d(GCCCGpGGC), containing one single 3'-methylene phosphonate linkage (denoted p), forms an A-DNA double helix with strict dyad symmetry, that is distinct from canonical A-DNA by a wide open major groove and a small average base-pair inclination against the helix axis. The conformation of the unmodified control d(GCCCGGGC) is known from an X-ray analysis of isomorphous crystals (Heinemann et al. (1987) Nucleic Acids Res. 15, 9531-9550). Comparison of the two structures reveals only minor conformational differences, most notably in the pucker of the reduced deoxyribose. It is suggested that oligonucleotides with charged 3'-methylene phosphonate groups may form stable duplexes with complementary DNA or RNA strands rendering them candidates for use as gene-regulatory antisense probes.