A first trial of retrospective collaboration for positional cloning in complex inheritance: Assay of the cytokine region on chromosome 5 by the Consortium on Asthma Genetics (COAG)

The central problem of complex inheritance is to map oligogenes for disease susceptibility, integrating linkage and association over samples that differ in several ways. Combination of evidence over multiple samples with 1,037 families supports loci contributing to asthma susceptibility in the cytokine region on 5q [maximum logarithm of odds (lod) = 2.61 near IL-4], but no evidence for atopy. The principal problems with retrospective collaboration on linkage appear to have been solved, providing far more information than a single study. A multipoint lod table evaluated at commonly agreed reference loci is required for both collaboration and metaanalysis, but variations in ascertainment, pedigree structure, phenotype definition, and marker selection are tolerated. These methods are invariant with statistical methods that increase the power of lods and are applicable to all diseases, motivating collaboration rather than competition. In contrast to linkage, positional cloning by allelic association has yet to be extended to multiple samples, a prerequisite for efficient combination with linkage and the greatest current challenge to genetic epidemiology.

The proliferative and antiapoptotic effects of substance P are facilitated by formation of a β-arrestin-dependent scaffolding complex

A requirement for scaffolding complexes containing internalized G protein-coupled receptors and β-arrestins in the activation and subcellular localization of extracellular signal-regulated kinases 1 and 2 (ERK1/2) has recently been proposed. However, the composition of these complexes and the importance of this requirement for function of ERK1/2 appear to differ between receptors. Here we report that substance P (SP) activation of neurokinin-1 receptor (NK1R) stimulates the formation of a scaffolding complex comprising internalized receptor, β-arrestin, src, and ERK1/2 (detected by gel filtration, immunoprecipitation, and immunofluorescence). Inhibition of complex formation, by expression of dominant-negative β-arrestin or a truncated NK1R that fails to interact with β-arrestin, inhibits both SP-stimulated endocytosis of the NK1R and activation of ERK1/2, which is required for the proliferative and antiapoptotic effects of SP. Thus, formation of a β-arrestin-containing complex facilitates the proliferative and antiapoptotic effects of SP, and these effects of SP could be diminished in cells expressing truncated NK1R corresponding to a naturally occurring variant.

Mode matches and their locations in the hydrophobic free energy sequences of peptide ligands and their receptor eigenfunctions

Patterns in sequences of amino acid hydrophobic free energies predict secondary structures in proteins. In protein folding, matches in hydrophobic free energy statistical wavelengths appear to contribute to selective aggregation of secondary structures in “hydrophobic zippers.” In a similar setting, the use of Fourier analysis to characterize the dominant statistical wavelengths of peptide ligands’ and receptor proteins’ hydrophobic modes to predict such matches has been limited by the aliasing and end effects of short peptide lengths, as well as the broad-band, mode multiplicity of many of their frequency (power) spectra. In addition, the sequence locations of the matching modes are lost in this transformation. We make new use of three techniques to address these difficulties: (i) eigenfunction construction from the linear decomposition of the lagged covariance matrices of the ligands and receptors as hydrophobic free energy sequences; (ii) maximum entropy, complex poles power spectra, which select the dominant modes of the hydrophobic free energy sequences or their eigenfunctions; and (iii) discrete, best bases, trigonometric wavelet transformations, which confirm the dominant spectral frequencies of the eigenfunctions and locate them as (absolute valued) moduli in the peptide or receptor sequence. The leading eigenfunction of the covariance matrix of a transmembrane receptor sequence locates the same transmembrane segments seen in n-block-averaged hydropathy plots while leaving the remaining hydrophobic modes unsmoothed and available for further analyses as secondary eigenfunctions. In these receptor eigenfunctions, we find a set of statistical wavelength matches between peptide ligands and their G-protein and tyrosine kinase coupled receptors, ranging across examples from 13.10 amino acids in acid fibroblast growth factor to 2.18 residues in corticotropin releasing factor. We find that the wavelet-located receptor modes in the extracellular loops are compatible with studies of receptor chimeric exchanges and point mutations. A nonbinding corticotropin-releasing factor receptor mutant is shown to have lost the signatory mode common to the normal receptor and its ligand. Hydrophobic free energy eigenfunctions and their transformations offer new quantitative physical homologies in database searches for peptide-receptor matches.

Advances in random matrix theory, zeta functions, and sphere packing

Over four hundred years ago, Sir Walter Raleigh asked his mathematical assistant to find formulas for the number of cannonballs in regularly stacked piles. These investigations aroused the curiosity of the astronomer Johannes Kepler and led to a problem that has gone centuries without a solution: why is the familiar cannonball stack the most efficient arrangement possible? Here we discuss the solution that Hales found in 1998. Almost every part of the 282-page proof relies on long computer verifications. Random matrix theory was developed by physicists to describe the spectra of complex nuclei. In particular, the statistical fluctuations of the eigenvalues (“the energy levels”) follow certain universal laws based on symmetry types. We describe these and then discuss the remarkable appearance of these laws for zeros of the Riemann zeta function (which is the generating function for prime numbers and is the last special function from the last century that is not understood today.) Explaining this phenomenon is a central problem. These topics are distinct, so we present them separately with their own introductory remarks.

Essential role for complex N-glycans in forming an organized layer of bronchial epithelium.

Mice lacking the complex subset of N-glycans due to inactivation of the Mgat1 gene die at mid-gestation, making it difficult to identify specific biological functions for this class of cell surface carbohydrates. To circumvent this embryonic lethality and to uncover tissue-specific functions for complex N-glycans, WW6 embryonic stem cells with inactivated Mgat1 alleles were tracked in chimeric embryos. The Mgat1 gene encodes N-acetylglucosaminyltransferase I (Glc-NAc-TI; EC 2.4.1.101), the transferase that initiates the synthesis of complex N-glycans. WW6 cells carry an inert beta-globin transgene that allows their identification in chimeras by DNA-DNA in situ hybridization. Independent Mgat1-/- and Mgat1+/- mutant WW6 isolates contributed like parent WW6 cells to the tissues of embryonic day (E) 10.5 to E16.5 chimeras. However, a cell type-specific difference was observed in lung. Homozygous null Mgat1-/- WW6 cells did not contribute to the epithelial layer in more than 99% bronchi. This deficiency was corrected by transfection of a Mgat1 transgene. Interestingly, heterozygous Mgat1+/- WW6 cells were also deficient in populating the layer of bronchial epithelium. Furthermore, examination of lung bud in E9.5 Mgat1-/- mutant embryos showed complete absence of an organized epithelial cell layer in the bronchus. Thus, complex N-glycans are required to form a morphologically recognizable bronchial epithelium, revealing an in vivo, cell type-specific function for this class of N-glycans.

Complex dynamics of multilocus systems subjected to cyclical selection.

Earlier we have shown that oscillations with a long period ("supercycles") may arise in two-locus systems experiencing cyclical selection with a short period. However, this mode of complex limiting behavior appeared to be possible for narrow ranges of parameters. Here we demonstrate that a multilocus system subjected to stabilizing selection with cyclically moving optimum can generate ubiquitous complex limiting behavior including supercycles, T-cycles, and chaotic-like phenomena. This mode of multilocus dynamics far exceeds the potential attainable under ordinary selection models resulting in simple behavior. It may represent a novel evolutionary mechanism increasing genetic diversity over long-term time periods.

Genes encoding the same three subunits of respiratory complex II are present in the mitochondrial DNA of two phylogenetically distant eukaryotes.

Although mitochondrial DNA is known to encode a limited number (<20) of the polypeptide components of respiratory complexes I, III, IV, and V, genes for components of complex II [succinate dehydrogenase (ubiquinone); succinate:ubiquinone oxidoreductase, EC 1.3.5.1] are conspicuously lacking in mitochondrial genomes so far characterized. Here we show that the same three subunits of complex II are encoded in the mitochondrial DNA of two phylogenetically distant eukaryotes, Porphyra purpurea (a photosynthetic red alga) and Reclinomonas americana (a heterotrophic zooflagellate). These complex II genes, sdh2, sdh3, and sdh4, are homologs, respectively, of Escherichia coli sdhB, sdhC, and sdhD. In E. coli, sdhB encodes the iron-sulfur subunit of succinate dehydrogenase (SDH), whereas sdhC and sdhD specify, respectively, apocytochrome b558 and a hydrophobic 13-kDa polypeptide, which together anchor SDH to the inner mitochondrial membrane. Amino acid sequence similarities indicate that sdh2, sdh3, and sdh4 were originally encoded in the protomitochondrial genome and have subsequently been transferred to the nuclear genome in most eukaryotes. The data presented here are consistent with the view that mitochondria constitute a monophyletic lineage.