An analysis on gene architecture in human and mouse genomes

A comparative genome analysis on exon-intron distribution profiles is performed for human and mouse genomes to deduce similarities and differences between them. Interestingly, both in human and mouse genomes, the total length in introns and intergenic DNA on each chromosome is significantly correlated to the chromosome size. The results presented provide a framework for understanding the nature and patterns of exon-intron length distributions, the constraints on them and their role in genome design and evolution.     

The fine-scale and complex architecture of human copy-number variation

Despite considerable excitement over the potential functional significance of copy-number variants (CNVs), we still lack knowledge of the fine-scale architecture of the large majority of CNV regions in the human genome. In this study, we used a high-resolution array-based comparative genomic hybridization (aCGH) platform that targeted known CNV regions of the human genome at approximately 1 kb resolution to interrogate the genomic DNAs of 30 individuals from four HapMap populations. Our results revealed that 1020 of 1153 CNV loci (88%) were actually smaller in size than what is recorded in the Database of Genomic Variants based on previously published studies. A reduction in size of more than 50% was observed for 876 CNV regions (76%). We conclude that the total genomic content of currently known common human CNVs is likely smaller than previously thought. In addition, approximately 8% of the CNV regions observed in multiple individuals exhibited genomic architectural complexity in the form of smaller CNVs within larger ones and CNVs with interindividual variation in breakpoints. Future association studies that aim to capture the potential influences of CNVs on disease phenotypes will need to consider how to best ascertain this previously uncharacterized complexity.     

Fine scale structural variants distinguish the genomes of Drosophila melanogaster and D. pseudoobscura

Background  

A primary objective of comparative genomics is to identify genomic elements of functional significance that contribute to phenotypic diversity. Complex changes in genome structure (insertions, duplications, rearrangements, translocations) may be widespread, and have important effects on organismal diversity. Any survey of genomic variation is incomplete without an assessment of structural changes.     

Conservation of fine-scale DNA marker order in the genomes of rice and the Triticeae.

DNA markers distribute over large chromosomal regions exhibit conservation of order (collinearity) in different cereal species, but it is not known whether this is maintained on a finer scale, i.e. < or = 2 cM. To address this, sets of two or more genetically linked DNA markers were localised to yeast artificial chromosomes containing rice DNA inserts. Linkage analysis of these DNA markers in barley revealed complete correspondence with their genetic order in rice, the distance between linked sequences on rice chromosomes being < 1.6 cM or < or = 1 + 10(6) bp (1 Mb). Thus, DNA markers separated in this range are collinear in rice, barley and, by inference, other members of the Triticeae. These results are discussed with respect to the use of rice as a key system for the isolation of cereal genes.     

The high-resolution architecture and structural dynamics of Bacillus spores

The capability to image single microbial cell surfaces at nanometer scale under native conditions would profoundly impact mechanistic and structural studies of pathogenesis, immunobiology, environmental resistance, and biotransformation. Here, using in vitro atomic force microscopy, we have directly visualized high-resolution native structures of bacterial endospores, including the exosporium and spore coats of four Bacillus species in air and water environments. Our results demonstrate that the mechanisms of spore coat self-assembly are similar to those described for inorganic and macromolecular crystallization. The dimensions of individual Bacillus atrophaeus spores decrease reversibly by 12% in response to a change in the environment from fully hydrated to air-dried state, establishing that the dormant spore is a dynamic physical structure. The interspecies distributions of spore length and width were determined for four species of Bacillus spores in water and air environments. The dimensions of individual spores differ significantly depending upon species, growth regimes, and environmental conditions. These findings may be useful in the reconstruction of environmental and physiological conditions during spore formation and for modeling the inhalation and dispersal of spores. This study provides a direct insight into molecular architecture and structural variability of bacterial endospores as a function of spatial and developmental organizational scales.     

The six hyaluronidase-like genes in the human and mouse genomes.

The human genome contains six hyaluronidase-like genes. Three genes (HYAL1, HYAL2 and HYAL3) are clustered on chromosome 3p21.3, and another two genes (HYAL4 and PH-20/SPAM1) and one expressed pseudogene (HYALP1) are similarly clustered on chromosome 7q31.3. The extensive homology between the different hyaluronidase genes suggests ancient gene duplication, followed by en masse block duplication, events that occurred before the emergence of modern mammals. Very recently we have found that the mouse genome also has six hyaluronidase-like genes that are also grouped into two clusters of three, in regions syntenic with the human genome. Surprisingly, the mouse ortholog of HYALP1 does not contain any mutations, and unlike its human counterpart may actually encode an active enzyme. Hyal-1 is the only hyaluronidase in mammalian plasma and urine, and is also found at high levels in major organs such as liver, kidney, spleen, and heart. A model is proposed suggesting that Hyal-2 and Hyal-1 are the major mammalian hyaluronidases in somatic tissues, and that they act in concert to degrade high molecular weight hyaluronan to the tetrasaccharide. Twenty-kDa hyaluronan fragments are generated at the cell surface in unique endocytic vesicles resulting from digestion by the glycosylphosphatidyl-inositol-anchored Hyal-2, transported intracellularly by an unknown process, and then further digested by Hyal-1. The two beta-exoglycosidases, beta-glucuronidase and beta-N-acetyl glucosaminidase, remove sugars from reducing termini of hyaluronan oligomers, and supplement the hyaluronidases in the catabolism of hyaluronan.     

Evolutionary pattern of protein architecture in mammal and fruit fly genomes

Mutations, which can alter amino acid constitution, contribute greatly to protein evolution. However, little is reported of their pattern during protein structural evolution. We investigated the distribution of non-synonymous single nucleotide polymorphisms (nsSNPs) and insertions/deletions (indels) along mammal and fruit fly proteins. We found the nsSNPs (and d(N)) and indels increased in protein boundary regions, and this pattern is inversely correlated with the distribution of protein domain density. Additionally, synonymous substitutions (and d(S)) are reduced in 5' and 3' regions, indicating more variable protein boundaries, compared with central interior. All evidence suggests that the inner part of coding sequences (CDSs) is comparatively conserved, whereas the 5' and 3' regions, with higher evolution rates, are more variable. We assumed that due to greater frequencies of nsSNPs and indels in adaptive regions of CDSs it could be easier to ultimately alter, gain, or lose amino acids, thus becoming the front line of protein evolution.     

The primary structure of genetic variants of mouse hemoglobin

The primary structures of the globins from CE/J, DBA/2J, and a stock of Potter's mice were determined to identify the amino acid substitutions associated with the unique isoelectric focusing patterns of these hemoglobins. In addition, the primary structures of the globins from MOL III and PERU mice were studied in search of amino acid substitutions that may not be detected by isoelectric focusing. CE/J hemoglobin contains a unique kind of globin called chain 5. It differs from the single kind of globin (chain 1) in C57BL/6 by having alanine rather than glycine at position 78. DAB/2J hemoglobin has two kinds of globins: one half is like chain 5 and the other half is like chain 1. The hemoglobin from Potter's stock of Mus musculus molossinus also contains chains 1 and 5, but they are expressed at different levels, i.e., 80% chain 1 and 20% chain 5. MOL III hemoglobin has a single kind of globin identical to that in C57BL/6, and PERU hemoglobin contains approximately 40% chain 1 and 60% chain 4. Chains 1 and 4 have different amino acids at positions 25, 62, and 68. These studies confirm that mouse hemoglobins separable by isoelectric focusing, but not by other means of electrophoresis, have substitutions of neutrally charged amino acids in their chains.This research was sponsored by the Office of Health and Environmental Research, U.S. Department of Energy, under Contract W-740r-eng-26 with the Union Carbide Corporation.     

Phosphoglucomutase electrophoretic variants in the mouse

The phosphoglucomutase (PGM) electrophoretic phenotype of the mouse (Mus musculus) consists of several distinct components which can be grouped into two major zones designated PGM-1 and PGM-2. Evidence presented here indicates that each zone is controlled by a single genetic locus denoted Pgm-1 and Pgm-2, respectively. Two variant forms segregated at the Pgm-1 locus. They were codominantly expressed and inherited as alleles at an autosomal locus. The alleles were termed Pgm-1 ^a (fast) and Pgm-1 ^b (slow). These alleles were separately fixed in a number of inbred strains of mice. Preliminary evidence based on wild mouse phenotypes indicates that variant forms also exist for PGM-2 which are inherited as alleles at an autosomal locus. Genetic linkage relationships have not been determined for these loci. PGM-1 variants and PGM-2 were expressed in mouse fibroblasts in vitro.Supported by U.S. Public Health Service grants GM-09966 and GM-07249 from General Medical Sciences and 5 F2 HD-35,531 from Child Health and Human Development; and Atomic Energy Commission contract AT(30-1)-3671.Postdoctoral Fellow of the U.S. Public Health Service.     

Histone variants in mouse centromeric chromatin

Summary Highly purified centromeric heterochromatin was isolated from mouse liver nuclei and the pattern of core histone variants was analyzed. In comparison with total chromatin, the centromeric heterochromatin of young animals was characterized by (1) enrichment in the replication-dependent variants H2A₁, H2B₂ and H3₂, (2) reduced amount of the minor variant H2A_z and (3) absence of ubiquitinated molecules of H2A. This specific variant pattern changed upon ageing as a result of accumulation of replacement variants so that in adult animals both chromatin preparations exhibited similar pattern for H2A and H2B, while the difference in the profile of H3 variants was preserved.     

Characterizing complex structural variation in germline and somatic genomes

Genome structural variation (SV) is a major source of genetic diversity in mammals and a hallmark of cancer. Although SV is typically defined by its canonical forms (duplication, deletion, insertion, inversion and translocation), recent breakpoint mapping studies have revealed a surprising number of 'complex' variants that evade simple classification. Complex SVs are defined by clustered breakpoints that arose through a single mutation but cannot be explained by one simple end-joining or recombination event. Some complex variants exhibit profoundly complicated rearrangements between distinct loci from multiple chromosomes, whereas others involve more subtle alterations at a single locus. These diverse and unpredictable features present a challenge for SV mapping experiments. Here, we review current knowledge of complex SV in mammals, and outline techniques for identifying and characterizing complex variants using next-generation DNA sequencing.     

Immunological and structural relatedness of isozymes and genetic variants of 3-phosphoglycerate kinase from the mouse

Isozymes (PGK-1 and PGK-2) and genetic variants (PGK-2A, PGK-2B, and PGK-2C) of 3-phosphoglycerate kinase were purified by affinity chromatography using an 8-(6-aminohexyl)-amino-ATP-Sepharose column as the key step. Antisera raised against purified PGK-1 and PGK-2A were tested for specificity and cross-reactivity by application of double immunodiffusion and enzyme immunoinactivation methods. By double immunodiffusion, no precipitin lines were observed between anti-PGK-2A and PGK-1, but a weak cross-reactivity between anti-PGK-1 and PGK-2A was detected. In addition to specific inhibition of PGK-1 and PGK-2A by their respective antisera, anti-PGK-1 was shown to inhibit PGK-2 activity at high antiserum concentrations, whereas no inhibition of PGK-1 activity by anti-PGK-2A was observed. The amino acid compositions of PGK-1 and PGK-2 revealed a certain degree of homology. However, tryptic peptide maps showed no obvious similarity in the peptide spots between these two 3-phosphoglycerate kinase isozymes. Three electrophoretic variants of PGK-2 were compared biochemically and immunologically. PGK-2C from C57L/J mice, a low activity variant, was shown to be the result of a structural gene mutation that affects the active site of the enzyme.     

The structural organization of mouse metaphase chromosomes

The binding of highly purified anti-nucleoside antibodies to mouse (Mus musculus) metaphase chromosomes was studied by an immunofluorescence technique. The chromosomal DNA was denatured by one of two selective denaturation procedures because these antibodies reacted with single stranded but not native DNA. After ultraviolet irradiation (UV), which produced single stranded regions primarily in AT rich DNA, the binding of antiadenosine (anti-A) produced a pattern of fluorescent bands similar to that produced by quinacrine (Q-bands). Additional foci of bright fluorescence were observed at the centrometric (C-band) regions, which are known to contain AT rich satellite DNA. After photooxidation, which produced single stranded regions in GC rich DNA, the binding of anti-A produced a fluorescent banding pattern similar to the R-banding pattern seen after thermal denaturation and staining with coriphosphine O. After photooxidation, R-band patterns were also obtained with anti-cytidine (anti-C) and anti-5-methylcytidine (anti-M). After either UV irradiation or photooxidation, anti-M, but not anti-C, showed intense binding to the C-band regions of mouse chromosomes. — These findings led to the following conclusions: (1) Antibody banding patterns reflect the presence of a class of AT rich, GC poor DNA in chromosome regions which show bright quinacrine fluorescence and in the regions that contain the AT rich satellite DNA. (2) The alternate, quinacrine dull regions contain a relatively GC rich class of DNA which appears to be more highly methylated than the AT rich DNA in the Q-bright bands, but not the AT rich satellite DNA in the Q-dull C-bands. (3) 5-Methylcytosine residues occur in a sequence of mouse satellite DNA that contains both adjacent pyrimidines and guanine residues. The basic repeating unit of mouse satellite DNA is known to contain the sequence 5-GAAAAATGA-3 (Biro et al., 1975). Therefore, assuming the antibodies used could detect single bases in denatured DNA, the methylated sequence in mouse satellite DNA      

Spatial separation of parental genomes in preimplantation mouse embryos

We have used two different experimental approaches to demonstrate topological separation of parental genomes in preimplantation mouse embryos: mouse eggs fertilized with 5-bromodeoxyuridine (BrdU)-labeled sperm followed by detection of BrdU in early diploid embryos, and differential heterochromatin staining in mouse interspecific hybrid embryos. Separation of chromatin according to parental origin was preserved up to the four-cell embryo stage and then gradually disappeared. In F1 hybrid animals, genome separation was also observed in a proportion of somatic cells. Separate nuclear compartments during preimplantation development, when extreme chromatin remodelling occurs, and possibly in some differentiated cell types, may be associated with epigenetic reprogramming.     

The genetic architecture of selection response. Inferences from fine-scale mapping of bristle number quantitative trait loci in Drosophila melanogaster.

Quantitative trait loci (QTL) affecting responses and correlated responses to selection for abdominal and sternopleural bristle number have been mapped with high resolution to the X and third chromosomes. Advanced intercross recombinant isogenic chromosomes were constructed from high and low selection lines in an unselected inbred background, and QTL were detected using composite interval mapping and high density transposable element marker maps. We mapped a total of 26 bristle number QTL with large effects, which were in or immediately adjacent to intervals previously inferred to contain bristle number QTL on these chromosomes. The QTL contributing to response to selection for high bristle number were not the same as those contributing to response to selection for low bristle number, suggesting that distributions of allelic effects per locus may be asymmetrical. Correlated responses were more often attributable to loose linkage than pleiotropy or close linkage. Bristle number QTL mapping to the same locations have been inferred in studies with different parental strains. Of the 26 QTL, 20 mapped to locations consistent with candidate genes affecting peripheral nervous system development and/or bristle number. This facilitates determining the molecular basis of quantitative variation and allele frequencies by associating molecular variation at the candidate genes with phenotypic variation in bristle number in samples of alleles from nature.     

An automated comparative analysis of 17 complete microbial genomes

MOTIVATION: As sequenced genomes become larger and sequencing becomes faster, there is a need to develop accurate automated genome comparison techniques and databases to facilitate derivation of genome functionality; identification of enzymes, putative operons and metabolic pathways; and to derive phylogenetic classification of microbes. RESULTS: This paper extends an automated pair-wise genome comparison technique (Bansal et al., Math. Model. Sci. Comput., 9, 1-23, 1998, Bansal and Bork, in First International Workshop of Declarative Languages, Springer, pp. 275-289, 1999) used to identify orthologs and gene groups to derive orthologous genes in a group of genomes and to identify genes with conserved functionality. Seventeen microbial genomes archived at ftp://ncbi.nlm.nih.gov/genbank/genomes have been compared using the automated technique. Data related to orthologs, gene groups, gene duplication, gene fusion, orthologs with conserved functionality, and genes specifically orthologous to Escherichia coli and pathogens has been presented and analyzed. AVAILABILITY: A prototype database is available at ftp://www.mcs.kent.edu/arvind/intellibio / orthos.html. The software is free for academic research under an academic license. The detailed database for every microbial genome in NCBI is commercially available through intellibio software and consultancy corporation (Web site: http://www.mcs.kent.edu/?rvind/intellibio . html). CONTACT: arvind@mcs.kent.edu.     

Correlation between sequence conservation and structural thermodynamics of microRNA precursors from human,mouse, and chicken genomes

Background

Previous studies have shown that microRNA precursors (pre-miRNAs) have considerably more stable secondary structures than other native RNAs (tRNA, rRNA, and mRNA) and artificial RNA sequences. However, pre-miRNAs with ultra stable secondary structures have not been investigated. It is not known if there is a tendency in pre-miRNA sequences towards or against ultra stable structures? Furthermore, the relationship between the structural thermodynamic stability of pre-miRNA and their evolution remains unclear.

Results

We investigated the correlation between pre-miRNA sequence conservation and structural stability as measured by adjusted minimum folding free energies in pre-miRNAs isolated from human, mouse, and chicken. The analysis revealed that conserved and non-conserved pre-miRNA sequences had structures with similar average stabilities. However, the relatively ultra stable and unstable pre-miRNAs were more likely to be non-conserved than pre-miRNAs with moderate stability. Non-conserved pre-miRNAs had more G+C than A+U nucleotides, while conserved pre-miRNAs contained more A+U nucleotides. Notably, the U content of conserved pre-miRNAs was especially higher than that of non-conserved pre-miRNAs. Further investigations showed that conserved and non-conserved pre-miRNAs exhibited different structural element features, even though they had comparable levels of stability.

Conclusions

We proposed that there is a correlation between structural thermodynamic stability and sequence conservation for pre-miRNAs from human, mouse, and chicken genomes. Our analyses suggested that pre-miRNAs with relatively ultra stable or unstable structures were less favoured by natural selection than those with moderately stable structures. Comparison of nucleotide compositions between non-conserved and conserved pre-miRNAs indicated the importance of U nucleotides in the pre-miRNA evolutionary process. Several characteristic structural elements were also detected in conserved pre-miRNAs.