Linkage analysis in the presence of errors I: complex-valued recombination fractions and complex phenotypes

Linkage is a phenomenon that correlates the genotypes of loci, rather than the phenotypes of one locus to the genotypes of another. It is therefore necessary to convert the observed trait phenotypes into trait-locus genotypes, which can then be analyzed for coinheritance with marker-locus genotypes. However, if the mode of inheritance of the trait is not known accurately, this conversion can often result in errors in the inferred trait-locus genotypes, which, in turn, can lead to the misclassification of the recombination status of meioses. As a result, the recombination fraction can be overestimated in two-point analysis, and false exclusions of the true trait locus can occur in multipoint analysis. We propose a method that increases the robustness of multipoint analysis to errors in the mode of inheritance assumptions of the trait, by explicitly allowing for misclassification of trait-locus genotypes. To this end, the definition of the recombination fraction is extended to the complex plane, as Theta=straight theta+straightepsiloni; theta is the recombination fraction between actual ("real") genotypes of marker and trait loci, and straightepsilon is the probability of apparent but false ("imaginary") recombinations between the actual and inferred trait-locus genotypes. "Complex" multipoint LOD scores are proven to be stochastically equivalent to conventional two-point LOD scores. The greater robustness to modeling errors normally associated with two-point analysis can thus be extended to multiple two-point analysis and multipoint analysis. The use of complex-valued recombination fractions also allows the stochastic equivalence of "model-based" and "model-free" methods to be extended to multipoint analysis.     

Two-trait-locus linkage analysis: a powerful strategy for mapping complex genetic traits.

Recent advances in molecular biology have provided geneticists with ever-increasing numbers of highly polymorphic genetic markers that have made possible linkage mapping of loci responsible for many human diseases. However, nearly all diseases mapped to date follow clear Mendelian, single-locus segregation patterns. In contrast, many common familial diseases such as diabetes, psoriasis, several forms of cancer, and schizophrenia are familial and appear to have a genetic component but do not exhibit simple Mendelian transmission. More complex models are required to explain the genetics of these important diseases. In this paper, we explore two-trait-locus, two-marker-locus linkage analysis in which two trait loci are mapped simultaneously to separate genetic markers. We compare the utility of this approach to standard one-trait-locus, one-marker-locus linkage analysis with and without allowance for heterogeneity. We also compare the utility of the two-trait-locus, two-marker-locus analysis to two-trait-locus, one-marker-locus linkage analysis. For common diseases, pedigrees are often bilineal, with disease genes entering via two or more unrelated pedigree members. Since such pedigrees often are avoided in linkage studies, we also investigate the relative information content of unilineal and bilineal pedigrees. For the dominant-or-recessive and threshold models that we consider, we find that two-trait-locus, two-marker-locus linkage analysis can provide substantially more linkage information, as measured by expected maximum lod score, than standard one-trait-locus, one-marker-locus methods, even allowing for heterogeneity, while, for a dominant-or-dominant generating model, one-locus models that allow for heterogeneity extract essentially as much information as the two-trait-locus methods. For these three models, we also find that bilineal pedigrees provide sufficient linkage information to warrant their inclusion in such studies. We also discuss strategies for assessing the significance of the two linkages assumed in two-trait-locus, two-marker-locus models.     

Domain orientation in the inactive response regulator Mycobacterium tuberculosis MtrA provides a barrier to activation

The structure of MtrA, an essential gene product for the human pathogen Mycobacterium tuberculosis, has been solved to a resolution of 2.1 A. MtrA is a member of the OmpR/PhoB family of response regulators and represents the fourth family member for which a structure of the protein in its inactive state has been determined. As is true for all OmpR/PhoB family members, MtrA possesses an N-terminal regulatory domain and a C-terminal winged helix-turn-helix DNA-binding domain, with phosphorylation of the regulatory domain modulating the activity of the protein. In the inactive form of MtrA, these two domains form an extensive interface that is composed of the alpha4-beta5-alpha5 face of the regulatory domain and the C-terminal end of the positioning helix, the trans-activation loop, and the recognition helix of the DNA-binding domain. This domain orientation suggests a mechanism of mutual inhibition by the two domains. Activation of MtrA would require a disruption of this interface to allow the alpha4-beta5-alpha5 face of the regulatory domain to form the intermolecule interactions that are associated with the active state and to allow the recognition helix to interact with DNA. Furthermore, the interface appears to stabilize the inactive conformation of MtrA, potentially reducing the rate of phosphorylation of the N-terminal domain. This combination of effects may form a switch, regulating the activity of MtrA. The domain orientation exhibited by MtrA also provides a rationale for the variation in linker length that is observed within the OmpR/PhoB family of response regulators.     

Rapid protein-folding assay using green fluorescent protein.

Formation of the chromophore of green fluorescent protein (GFP) depends on the correct folding of the protein. We constructed a "folding reporter" vector, in which a test protein is expressed as an N-terminal fusion with GFP. Using a test panel of 20 proteins, we demonstrated that the fluorescence of Escherichia coli cells expressing such GFP fusions is related to the productive folding of the upstream protein domains expressed alone. We used this fluorescent indicator of protein folding to evolve proteins that are normally prone to aggregation during expression in E. coli into closely related proteins that fold robustly and are fully soluble and functional. This approach to improving protein folding does not require functional assays for the protein of interest and provides a simple route to improving protein folding and expression by directed evolution.     

Benefits of recombinant adeno-associated virus (rAAV)-mediated insulinlike growth factor I (IGF-I) overexpression for the long-term reconstruction of human osteoarthritic cartilage by modulation of the IGF-I axis

Administration of therapeutic genes to human osteoarthritic (OA) cartilage is a potential approach to generate effective, durable treatments against this slow, progressive disorder. Here, we tested the ability of recombinant adeno-associated virus (rAAV)-mediated overexpression of human insulinlike growth factor (hIGF)-I to reproduce an original surface in human OA cartilage in light of the pleiotropic activities of the factor. We examined the proliferative, survival and anabolic effects of the rAAV-hIGF-I treatment in primary human normal and OA chondrocytes in vitro and in explant cultures in situ compared with control (reporter) vector delivery. Efficient, prolonged IGF-I secretion via rAAV stimulated the biological activities of OA chondrocytes in all the systems evaluated over extended periods of time, especially in situ, where it allowed for the long-term reconstruction of OA cartilage (at least for 90 d). Remarkably, production of high, stable amounts of IGF-I in OA cartilage using rAAV advantageously modulated the expression of central effectors of the IGF-I axis by downregulating IGF-I inhibitors (IGF binding protein [IGFBP]-3 and IGFBP4) while up-regulating key potentiators (IGFBP5, the IGF-I receptor and downstream mitogen-activated protein kinase/extracellular signal-regulated kinase 1/2 [MAPK/ERK-1/2] and phosphatidylinisitol-3/Akt [PI3K/Akt] signal transduction pathways), probably explaining the enhanced responsiveness of OA cartilage to IGF-I treatment. These findings show the benefits of directly providing an IGF-I sequence to articular cartilage via rAAV for the future treatment of human osteoarthritis.     

Recent Advances in GFP Folding Reporter and Split-GFP Solubility Reporter Technologies. Application to Improving the Folding and Solubility of Recalcitrant Proteins from Mycobacterium tuberculosis

We have improved our green fluorescent protein (GFP) folding reporter technology [Waldo et al., (1999) Nat. Biotechnol. 17, 691–695] to evolve recalcitrant proteins from Mycobacterium tuberculosis. The target protein is inserted into the scaffolding of the GFP, eliminating false-positive artifacts caused by expression of truncated protein variants from internal cryptic ribosome binding sites in the target RNA. In parallel, we have developed a new quantitative fluorescent protein tagging and detection system based on micro-domains of GFP. This split-GFP system, which works both in vivo and in vitro, is amenable to high-throughput assays of protein expression and solubility [Cabantous et al., (2005) Nat. Biotechnol. 23, 102–107]. Together, the GFP folding reporter and split-GFP technologies offer a comprehensive system for manipulating and improving protein folding and solubility.     

Subfamily-Specific Adaptations in the Structures of Two Penicillin-Binding Proteins from Mycobacterium tuberculosis

Beta-lactam antibiotics target penicillin-binding proteins including several enzyme classes essential for bacterial cell-wall homeostasis. To better understand the functional and inhibitor-binding specificities of penicillin-binding proteins from the pathogen, Mycobacterium tuberculosis, we carried out structural and phylogenetic analysis of two predicted D,D-carboxypeptidases, Rv2911 and Rv3330. Optimization of Rv2911 for crystallization using directed evolution and the GFP folding reporter method yielded a soluble quadruple mutant. Structures of optimized Rv2911 bound to phenylmethylsulfonyl fluoride and Rv3330 bound to meropenem show that, in contrast to the nonspecific inhibitor, meropenem forms an extended interaction with the enzyme along a conserved surface. Phylogenetic analysis shows that Rv2911 and Rv3330 belong to different clades that emerged in Actinobacteria and are not represented in model organisms such as Escherichia coli and Bacillus subtilis. Clade-specific adaptations allow these enzymes to fulfill distinct physiological roles despite strict conservation of core catalytic residues. The characteristic differences include potential protein-protein interaction surfaces and specificity-determining residues surrounding the catalytic site. Overall, these structural insights lay the groundwork to develop improved beta-lactam therapeutics for tuberculosis.