The optimization principle in phylogenetic analysis tends to give incorrect topologies when the number of nucleotides or amino acids used is small

In the maximum parsimony (MP) and minimum evolution (ME) methods of phylogenetic inference, evolutionary trees are constructed by searching for the topology that shows the minimum number of mutational changes required (M) and the smallest sum of branch lengths (S), respectively, whereas in the maximum likelihood (ML) method the topology showing the highest maximum likelihood (A) of observing a given data set is chosen. However, the theoretical basis of the optimization principle remains unclear. We therefore examined the relationships of M, S, and A for the MP, ME, and ML trees with those for the true tree by using computer simulation. The results show that M and S are generally greater for the true tree than for the MP and ME trees when the number of nucleotides examined (n) is relatively small, whereas A is generally lower for the true tree than for the ML tree. This finding indicates that the optimization principle tends to give incorrect topologies when n is small. To deal with this disturbing property of the optimization principle, we suggest that more attention should be given to testing the statistical reliability of an estimated tree rather than to finding the optimal tree with excessive efforts. When a reliability test is conducted, simplified MP, ME, and ML algorithms such as the neighbor-joining method generally give conclusions about phylogenetic inference very similar to those obtained by the more extensive tree search algorithms.

Telomeres in the mouse have large inter-chromosomal variations in the number of T2AG3 repeats

The ultra-long telomeres that have been observed in mice are not in accordance with the concept that critical telomere shortening is related to aging and immortalization. Here, we have used quantitative fluorescence in situ hybridization to estimate (T2AG3)n lengths of individual telomeres in various mouse strains. Telomere lengths were very heterogeneous, but specific chromosomes of bone marrow cells and skin fibroblasts from individual mice had similar telomere lengths. We estimate that the shortest telomeres are around 10 kb in length, indicating that each mouse cell has a few telomeres with (T2AG3)n lengths within the range of human telomeres. These short telomeres may be critical in limiting the replicative potential of murine cells.

Varying the number of telomere-bound proteins does not alter telomere length in tel1Δ cells

Yeast telomere DNA consists of a continuous, ≈330-bp tract of the heterogeneous repeat TG1-3 with irregularly spaced, high affinity sites for the protein Rap1p. Yeast monitor, or count, the number of telomeric Rap1p C termini in a negative feedback mechanism to modulate the length of the terminal TG1-3 repeats, and synthetic telomeres that tether Rap1p molecules adjacent to the TG1-3 tract cause wild-type cells to maintain a shorter TG1-3 tract. To identify trans-acting proteins required to count Rap1p molecules, these same synthetic telomeres were placed in two short telomere mutants: yku70Δ (which lack the yeast Ku70 protein) and tel1Δ (which lack the yeast ortholog of ATM). Although both mutants maintain telomeres with ≈100 bp of TG1-3, only yku70Δ cells maintained shorter TG1-3 repeats in response to internal Rap1p molecules. This distinct response to internal Rap1p molecules was not caused by a variation in Rap1p site density in the TG1-3 repeats as sequencing of tel1Δ and yku70Δ telomeres showed that both strains have only five to six Rap1p sites per 100-bp telomere. In addition, the tel1Δ short telomere phenotype was epistatic to the unregulated telomere length caused by deletion of the Rap1p C-terminal domain. Thus, the length of the TG1-3 repeats in tel1Δ cells was independent of the number of the Rap1p C termini at the telomere. These data indicate that tel1Δ cells use an alternative mechanism to regulate telomere length that is distinct from monitoring the number of telomere binding proteins.

A very limited number of keywords (main patterns) describes all sequences of the human variable heavy (VH) and κ (Vκ) domains

Sequences of the variable heavy (VH) and κ (Vκ) domains of Ig structures were divided into 21 fragments that correspond to strands, loops, or parts of these structural units of the variable domains. Amino acid sequences of fragments (termed “words”) were collected from the 1,172 human heavy and 668 human κ chains available in the Kabat database. Statistical analysis of words of 17 fragments was performed (fragments that comprise the complementary determining regions′ fragments will not be discussed in this paper). The number of different words (those with different residues in at least one position) ranged, for various fragments, from 11 to 75 in the κ chains, and from 23 to 189 in the heavy chains. The main result of this study is that very few keywords, or main patterns of words, were necessary to describe over 90% of the sequences (no more than two keywords per fragment in the κ and no more than five per fragment in the heavy chains). No identical keywords were found for different fragments of the variable domains. Keywords of aligned fragments of the VH and Vκ domains were different in all but two instances. Thus, knowing the keywords, one can determine whether any given small part of a sequence belongs to a heavy or κ chain and predict its precise localization in the sequence. In addition, by using all of the keywords obtained through analysis of the Kabat database, it was possible to describe completely the sequences of the human VH and Vκ germ-line segments.

The number of unidentified amacrine cells in the mammalian retina

The three largest known populations of amacrine cells in the rabbit retina were stained with fluorescent probes in whole mounts and counted at a series of retinal eccentricities. The retinas were counterstained using a fluorescent DNA-binding molecule and the total number of nuclei in the inner nuclear layer were counted in confocal sections. From the total number of inner nuclear layer cells and the known fraction of them occupied by amacrine cells, the fraction of amacrine cells made up by the stained populations could be calculated. Starburst cells made up 3%, indoleamine-accumulating cells made up 4%, and AII cells made up 11% of all amacrine cells. By referring four smaller populations of amacrine cells to the number of indoleamine-accumulating cells, they were estimated to make up 4% of all amacrine cells. Thus, 78% of all amacrine cells in the rabbit’s retina are known only from isolated examples, if at all. This proportion is similar in the retinas of the mouse, cat, and monkey. It is likely that a substantial fraction of the local circuit neurons present in other regions of the central nervous system are also invisible as populations to current techniques.

Heat effects on DNA repair after ionising radiation: hyperthermia commonly increases the number of non-repaired double-strand breaks and structural rearrangements

After ionising radiation double-strand breaks (dsb) are lethal if not repaired or misrepaired. Cell killing is greatly enhanced by hyperthermia and it is questioned here whether heat not only affects dsb repair capacity but also fidelity in a chromosomal context. dsb repair experiments were designed so as to mainly score non-homologous end joining, while homologous recombination was largely precluded. Human male G0 fibroblasts were either preheated (45°C, 20 min) or not before X-irradiation. dsb induction and repair were measured by conventional gel electrophoresis and an assay combining restriction digestion using a rare cutting enzyme (NotI) and Southern hybridisation, which detects large chromosomal rearrangements (>100 kb). dsb induction rate in an X-chromosomal NotI fragment was 4.8 × 10–3 dsb/Gy/Mb. Similar values were found for the genome overall and also when cells were preheated. After 50 Gy, fibroblasts were competent to largely restore the original restriction fragment size. Five per cent of dsb remained non-rejoined and 14% were misrejoined. Correct restitution of restriction fragments occurred preferably during the first hour but continued at a slow rate for 12–16 h. In addition, dsb appeared to misrejoin throughout the entire repair period. After hyperthermia the fractions of non-rejoined and misrejoined dsb were similarly increased to 13 and 51%, respectively. It is suggested that heat increases the probability of dsb being incorrectly rejoined but it is not likely to interfere with one dsb repair pathway in particular.

Multiplex genotype determination at a large number of gene loci.

To facilitate large-scale genotype analysis, an efficient PCR-based multiplex approach has been developed. For simultaneously amplifying the target sequences at a large number of genetic loci, locus-specific primers containing 5' universal tails are used. Attaching the universal tails to the target sequences in the initial PCR steps allows replacement of all specific primers with a pair of primers identical to the universal tails and converts the multiplex amplification into "uniplex." Simultaneous amplification of 26 genetic loci with this approach is described. The multiplex amplification can be coupled with genotype determination. By incorporating a single-base mismatch between a primer and the template into the target sequences, a polymorphic site can be converted into a desirable restriction fragment length polymorphism when it is necessary. In this way, the allelic PCR products for the polymorphic loci can be discriminated by gel electrophoresis after restriction enzyme digestion. In this study, 32 loci were typed in such a multiplex way.

Cisplatin increases meiotic crossing-over in mice

Genetic mapping of traits and mutations in mammals is dependent upon linkage analysis. The resolution achieved by this method is related to the number of offspring that can be scored and position of crossovers near a gene. Higher precision mapping is obtained by expanding the collection of progeny from an appropriate cross, which in turn increases the number of potentially informative recombinants. A more efficient approach would be to increase the frequency of recombination, rather than the number of progeny. The anticancer drug cisplatin, which causes DNA strand breakage and is highly recombinogenic in some model organisms, was tested for its ability to induce germ-line recombination in mice. Males were exposed to cisplatin and mated at various times thereafter to monitor the number of crossovers inherited by offspring. We observed a striking increase on all three chromosomes examined and established a regimen that nearly doubled crossover frequency. The timing of the response indicated that the crossovers were induced at the early pachytene stage of meiosis I. The ability to increase recombination should facilitate genetic mapping and positional cloning in mice.

Typing of urinary JC virus DNA offers a novel means of tracing human migrations

Although polyomavirus JC (JCV) is the proven pathogen of progressive multifocal leukoencephalopathy, the fatal demyelinating disease, this virus is ubiquitous as a usually harmless symbiote among human beings. JCV propagates in the adult kidney and excretes its progeny in urine, from which JCV DNA can readily be recovered. The main mode of transmission of JCV is from parents to children through long cohabitation. In this study, we collected a substantial number of urine samples from native inhabitants of 34 countries in Europe, Africa, and Asia. A 610-bp segment of JCV DNA was amplified from each urine sample, and its DNA sequence was determined. A worldwide phylogenetic tree subsequently constructed revealed the presence of nine subtypes including minor ones. Five subtypes (EU, Af2, B1, SC, and CY) occupied rather large territories that overlapped with each other at their boundaries. The entire Europe, northern Africa, and western Asia were the domain of EU, whereas the domain of Af2 included nearly all of Africa and southwestern Asia all the way to the northeastern edge of India. Partially overlapping domains in Asia were occupied by subtypes B1, SC, and CY. Of particular interest was the recovery of JCV subtypes in a pocket or pockets that were separated by great geographic distances from the main domains of those subtypes. Certain of these pockets can readily be explained by recent migrations of human populations carrying these subtypes. Overall, it appears that JCV genotyping promises to reveal previously unknown human migration routes: ancient as well as recent.

Quantitative insight into proliferation and differentiation of oligodendrocyte type 2 astrocyte progenitor cells in vitro

As part of our attempts at understanding fundamental principles that underlie the generation of nondividing terminally differentiated progeny from dividing precursor cells, we have developed approaches to a quantitative analysis of proliferation and differentiation of oligodendrocyte type 2 astrocyte (O-2A) progenitor cells at the clonal level. Owing to extensive previous studies of clonal differentiation in this lineage, O-2A progenitor cells represent an excellent system for such an analysis. Previous studies have resulted in two competing hypotheses; one of them suggests that progenitor cell differentiation is symmetric, the other hypothesis introduces an asymmetric process of differentiation. We propose a general model that incorporates both such extreme hypotheses as special cases. Our analysis of experimental data has shown, however, that neither of these extreme cases completely explains the observed kinetics of O-2A progenitor cell proliferation and oligodendrocyte generation in vitro. Instead, our results indicate that O-2A progenitor cells become competent for differentiation after they complete a certain number of critical mitotic cycles that represent a period of symmetric development. This number varies from clone to clone and may be thought of as a random variable; its probability distribution was estimated from experimental data. Those O-2A cells that have undergone the critical divisions then may differentiate into an oligodendrocyte in each of the subsequent mitotic cycles with a certain probability, thereby exhibiting the asymmetric type of differentiation.

Nuclear Pore Complex Number and Distribution throughout the Saccharomyces cerevisiae Cell Cycle by Three-Dimensional Reconstruction from Electron Micrographs of Nuclear Envelopes

The number of nuclear pore complexes (NPCs) in individual nuclei of the yeast Saccharomyces cerevisiae was determined by computer-aided reconstruction of entire nuclei from electron micrographs of serially sectioned cells. Nuclei of 32 haploid cells at various points in the cell cycle were modeled and found to contain between 65 and 182 NPCs. Morphological markers, such as cell shape and nuclear shape, were used to determine the cell cycle stage of the cell being examined. NPC number was correlated with cell cycle stage to reveal that the number of NPCs increases steadily, beginning in G1-phase, suggesting that NPC assembly occurs continuously throughout the cell cycle. However, the accumulation of nuclear envelope observed during the cell cycle, indicated by nuclear surface area, is not continuous at the same rate, such that the density of NPCs per unit area of nuclear envelope peaks in apparent S-phase cells. Analysis of the nuclear envelope reconstructions also revealed no preferred NPC-to-NPC distance. However, NPCs were found in large clusters over regions of the nuclear envelope. Interestingly, clusters of NPCs were most pronounced in early mitotic nuclei and were found to be associated with the spindle pole bodies, but the functional significance of this association is unknown.

The Malthusian parameter of ascents: What prevents the exponential increase of one’s ancestors?

The reason that the indefinite exponential increase in the number of one’s ancestors does not take place is found in the law of sibling interference, which can be expressed by the following simple equation:\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}\begin{matrix}{\mathit{N}}_{{\mathit{n}}} \enskip & \\ {\mathit{{\blacksquare}}} \enskip & \\ {\mathit{ASZ}} \enskip & \end{matrix} {\mathrm{\hspace{.167em}{\times}\hspace{.167em}2\hspace{.167em}=\hspace{.167em}}}{\mathit{N_{n+1},}}\end{equation*}\end{document} where Nn is the number of ancestors in the nth generation, ASZ is the average sibling size of these ancestors, and Nn+1 is the number of ancestors in the next older generation (n + 1). Accordingly, the exponential increase in the number of one’s ancestors is an initial anomaly that occurs while ASZ remains at 1. Once ASZ begins to exceed 1, the rate of increase in the number of ancestors is progressively curtailed, falling further and further behind the exponential increase rate. Eventually, ASZ reaches 2, and at that point, the number of ancestors stops increasing for two generations. These two generations, named AN SA and AN SA + 1, are the most critical in the ancestry, for one’s ancestors at that point come to represent all the progeny-produced adults of the entire ancestral population. Thereafter, the fate of one’s ancestors becomes the fate of the entire population. If the population to which one belongs is a successful, slowly expanding one, the number of ancestors would slowly decline as you move toward the remote past. This is because ABZ would exceed 2. Only when ABZ is less than 2 would the number of ancestors increase beyond the AN SA and AN SA + 1 generations. Since the above is an indication of a failing population on the way to extinction, there had to be the previous AN SA involving a far greater number of individuals for such a population. Simulations indicated that for a member of a continuously successful population, the AN SA ancestors might have numbered as many as 5.2 million, the AN SA generation being the 28th generation in the past. However, because of the law of increasingly irrelevant remote ancestors, only a very small fraction of the AN SA ancestors would have left genetic traces in the genome of each descendant of today.

Regulation of number and size of digits by posterior Hox genes: A dose-dependent mechanism with potential evolutionary implications

The proper development of digits, in tetrapods, requires the activity of several genes of the HoxA and HoxD homeobox gene complexes. By using a variety of loss-of-function alleles involving the five Hox genes that have been described to affect digit patterning, we report here that the group 11, 12, and 13 genes control both the size and number of murine digits in a dose-dependent fashion, rather than through a Hox code involving differential qualitative functions. A similar dose–response is observed in the morphogenesis of the penian bone, the baculum, which further suggests that digits and external genitalia share this genetic control mechanism. A progressive reduction in the dose of Hox gene products led first to ectrodactyly, then to olygodactyly and adactyly. Interestingly, this transition between the pentadactyl to the adactyl formula went through a step of polydactyly. We propose that in the distal appendage of polydactylous short-digited ancestral tetrapods, such as Acanthostega, the HoxA complex was predominantly active. Subsequent recruitment of the HoxD complex contributed to both reductions in digit number and increase in digit length. Thus, transition through a polydactylous limb before reaching and stabilizing the pentadactyl pattern may have relied, at least in part, on asynchronous and independent changes in the regulation of HoxA and HoxD gene complexes.

Number needed to screen: development of a statistic for disease screening

Objectives: To develop the number needed to screen, a new statistic to overcome inappropriate national strategies for disease screening. Number needed to screen is defined as the number of people that need to be screened for a given duration to prevent one death or adverse event.

Comprehensive copy number and gene expression profiling of the 17q23 amplicon in human breast cancer

The biological significance of DNA amplification in cancer is thought to be due to the selection of increased expression of a single or few important genes. However, systematic surveys of the copy number and expression of all genes within an amplified region of the genome have not been performed. Here we have used a combination of molecular, genomic, and microarray technologies to identify target genes for 17q23, a common region of amplification in breast cancers with poor prognosis. Construction of a 4-Mb genomic contig made it possible to define two common regions of amplification in breast cancer cell lines. Analysis of 184 primary breast tumors by fluorescence in situ hybridization on tissue microarrays validated these results with the highest amplification frequency (12.5%) observed for the distal region. Based on GeneMap'99 information, 17 known genes and 26 expressed sequence tags were localized to the contig. Analysis of genomic sequence identified 77 additional transcripts. A comprehensive analysis of expression levels of these transcripts in six breast cancer cell lines was carried out by using complementary DNA microarrays. The expression patterns varied from one cell line to another, and several overexpressed genes were identified. Of these, RPS6KB1, MUL, APPBP2, and TRAP240 as well as one uncharacterized expressed sequence tag were located in the two common amplified regions. In summary, comprehensive analysis of the 17q23 amplicon revealed a limited number of highly expressed genes that may contribute to the more aggressive clinical course observed in breast cancer patients with 17q23-amplified tumors.