Leveraging Distant Relatedness to Quantify Human Mutation and Gene-Conversion Rates

The rate at which human genomes mutate is a central biological parameter that has many implications for our ability to understand demographic and evolutionary phenomena. We present a method for inferring mutation and gene-conversion rates by using the number of sequence differences observed in identical-by-descent (IBD) segments together with a reconstructed model of recent population-size history. This approach is robust to, and can quantify, the presence of substantial genotyping error, as validated in coalescent simulations. We applied the method to 498 trio-phased sequenced Dutch individuals and inferred a point mutation rate of 1.66 × 10⁻⁸ per base per generation and a rate of 1.26 × 10⁻⁹ for <20 bp indels. By quantifying how estimates varied as a function of allele frequency, we inferred the probability that a site is involved in non-crossover gene conversion as 5.99 × 10⁻⁶. We found that recombination does not have observable mutagenic effects after gene conversion is accounted for and that local gene-conversion rates reflect recombination rates. We detected a strong enrichment of recent deleterious variation among mismatching variants found within IBD regions and observed summary statistics of local sharing of IBD segments to closely match previously proposed metrics of background selection; however, we found no significant effects of selection on our mutation-rate estimates. We detected no evidence of strong variation of mutation rates in a number of genomic annotations obtained from several recent studies. Our analysis suggests that a mutation-rate estimate higher than that reported by recent pedigree-based studies should be adopted in the context of DNA-based demographic reconstruction.     

Atomic Resolution Cryo-EM Structure of β-Galactosidase

1. Download high-res image (381KB)
2. Download full-size image

     

Morphing the feature-based multi-blocks of normative/healthy vertebral geometries to scoliosis vertebral geometries: development of personalized finite element models

Personalized Finite Element (FE) models and hexahedral elements are preferred for biomechanical investigations. Feature-based multi-block methods are used to develop anatomically accurate personalized FE models with hexahedral mesh. It is tedious to manually construct multi-blocks for large number of geometries on an individual basis to develop personalized FE models. Mesh-morphing method mitigates the aforementioned tediousness in meshing personalized geometries every time, but leads to element warping and loss of geometrical data. Such issues increase in magnitude when normative spine FE model is morphed to scoliosis-affected spinal geometry. The only way to bypass the issue of hex-mesh distortion or loss of geometry as a result of morphing is to rely on manually constructing the multi-blocks for scoliosis-affected spine geometry of each individual, which is time intensive. A method to semi-automate the construction of multi-blocks on the geometry of scoliosis vertebrae from the existing multi-blocks of normative vertebrae is demonstrated in this paper. High-quality hexahedral elements were generated on the scoliosis vertebrae from the morphed multi-blocks of normative vertebrae. Time taken was 3 months to construct the multi-blocks for normative spine and less than a day for scoliosis. Efforts taken to construct multi-blocks on personalized scoliosis spinal geometries are significantly reduced by morphing existing multi-blocks.     

Two Novel Human Members of an Emerging Mammalian Gene Family Related to Mono-ADP-Ribosylating Bacterial Toxins

Mono-ADP-ribosylation is one of the posttranslational protein modifications regulating cellular metabolism, e.g., nitrogen fixation, in prokaryotes. Several bacterial toxins mono-ADP-ribosylate and inactivate specific proteins in their animal hosts. Recently, two mammalian GPI-anchored cell surface enzymes with similar activities were cloned (designated ART1 and ART2). We have now identified six related expressed sequence tags (ESTs) in the public database and cloned the two novel human genes from which these are derived (designatedART3andART4). The deduced amino acid sequences of the predicted gene products show 28% sequence identity to one another and 32–41% identity vs the muscle and T cell enzymes. They contain signal peptide sequences characteristic of GPI anchorage. Southern Zoo blot analyses suggest the presence of related genes in other mammalian species. By PCR screening of somatic cell hybrids and byin situhybridization, we have mapped the two genes to human chromosomes 4p14–p15.1 and 12q13.2–q13.3. Northern blot analyses show that these genes are specifically expressed in testis and spleen, respectively. Comparison of genomic and cDNA sequences reveals a conserved exon/intron structure, with an unusually large exon encoding the predicted mature membrane proteins. Secondary structure prediction analyses indicate conserved motifs and amino acid residues consistent with a common ancestry of this emerging mammalian enzyme family and bacterial mono(ADP-ribosyl)transferases. It is possible that the four human gene family members identified so far represent the “tip of an iceberg,” i.e., a larger family of enzymes that influences the function of target proteins via mono-ADP-ribosylation.     

Associations of <Emphasis Type="Italic">CTLA4</Emphasis> +49 A/G Dimorphism and HLA-DRB1*/DQB1* Alleles With Type 1 Diabetes from South India

The aim of present study was to elucidate the association of CTLA4 +49 A/G and HLA-DRB1*/DQB1* gene polymorphism in south Indian T1DM patients. The patients and controls (n?=?196 each) were enrolled for CTLA4 and HLA-DRB1*/DQB1* genotyping by RFLP/PCR-SSP methods. The increased frequencies of CTLA4 ‘AG’ (OR?=?1.99; p?=?0.001), ‘GG’ (OR?=?3.94; p?=?0.001) genotypes, and ‘G’ allele (OR?=?2.42; p?=?9.26?×?10^?8) were observed in patients. Reduced frequencies of ‘AA’ (OR?=?0.35; p?=?7.19?×?10^?7) and ‘A’ (OR?=?0.41; p?=?9.26?×?10^?8) in patients revealed protective association. Among HLA-DRB1*/DQB1* alleles, DRB1*04 (OR?=?3.29; p?=?1.0?×?10^?5), DRB1*03 (OR?=?2.81; p?=?1.9?×?10^?6), DQB1*02:01 (OR?=?2.93; p?=?1.65?×?10^?5), DQB1*02:02 (OR?=?3.38; p?=?0.0003), and DQB1*03:02 (OR?=?7.72; p?=?0.0003) were in susceptible association. Decreased frequencies of alleles, DRB1*15 (OR?=?0.32; p?=?2.55?×?10^?7), DRB1*10 (OR?=?0.45; p?=?0.002), DQB1*06:01 (OR?=?0.43; p?=?0.0001), and DQB1*05:02 (OR?=?0.28; p?=?2.1?×?10^?4) in patients were suggested protective association. The combination of DRB1*03+AG (OR?=?5.21; p?=?1.4?×?10^?6), DRB1*04+AG (OR?=?2.14; p?=?0.053), DRB1*04+GG (OR?=?5.21; p?=?0.036), DQB1*02:01+AG (OR?=?4.44; p?=?3.6?×?10^?5), DQB1*02:02+AG (OR?=?20.9; p?=?9.5?×?10^?4), and DQB1*02:02+GG (OR?=?4.06; p?=?0.036) revealed susceptible association. However, the combination of DRB1*10+AA (OR?=?0.35; p?=?0.003), DRB1*15+AA (OR?=?0.22; p?=?5.3?×?10^?7), DQB1*05:01+AA (OR?=?0.45; p?=?0.007), DQB1*05:02+AA (OR?=?0.17; p?=?1.7?×?10^?4), DQB1*06:01+AA (OR?=?0.40; p?=?0.002), and DQB1*06:02+AG (OR?=?0.34; p?=?0.001) showed decreased frequency in patients, suggesting protective association. In conclusion, CTLA4/HLA-DR/DQ genotypic combinations revealed strong susceptible/protective association toward T1DM in south India. A female preponderance in disease associations was also documented.     

Efficacy and specificity of antisense laminin chain-specific expression vectors in blocking laminin induction by TGFβ1: Effect of laminin blockade on TGFβ1-mediated cellular responses

Transforming growth factorβ1 (TGFβ1) elicits a multitude of cellular responses from the epithelial-derived human colon cancer Moser cells. TGFβ1 induces the expression of laminin and fibronectin, and previous studies show that the induction of fibronectin is functionally associated with the regulation of carcinoembryonic antigen (CEA) expression by TGFβ1 (Huang and Chakrabarty, 1994, J Biol Chem 269:28764–28768). In this study we constructed antisense laminin chain-specific expression vectors and determined their efficacy in blocking the expression and the induction of the large multichain laminin molecule by TGFβ1. We also determined the functional role of laminin in several TGFβ1-mediated responses: growth inhibition, downmodulation of anchorage-independent growth, and cellular invasion. Expression of either antisense laminin chain A, B1, or B2 RNA resulted in a downmodulation of endogenous laminin mRNA expression and blocked the induction of laminin protein by TGFβ1 without affecting the induction of other adhesion molecules such as fibronectin or CEA. It is concluded that antisense RNA directed to only one of the laminin chains was sufficient to disrupt the induction of the complex laminin molecule in quite a specific manner. Expression of antisense laminin RNA downregulated cellular adhesion to extracellular matrix (ECM) laminin and blocked the ability of TGFβ1 to upmodulate adhesion to ECM laminin. Expression of antisense laminin RNA, however, did not alter the downregulating effect of TGFβ1 on cellular proliferation, anchorage-independent growth, or cellular invasion, suggesting that the induction of laminin did not play a significant functional role in these TGFβ1-mediated cellular responses. It is likely that other adhesion pathways may be involved in mediating the action of TGFβ1 in this cell line. J. Cell. Physiol. 178:296–303, 1999. © 1999 Wiley-Liss, Inc.     

Stoichiometry and absolute quantification of proteins with mass spectrometry using fluorescent and isotope-labeled concatenated peptide standards

We have explored a general approach for the determination of absolute amounts and the relative stoichiometry of proteins in a mixture using fluorescence and mass spectrometry. We engineered a gene to express green fluorescent protein (GFP) with a synthetic fusion protein (GAB-GFP) in Escherichia coli to function as a spectroscopic standard for the quantification of an analogous stable isotope-labeled, non-fluorescent fusion protein (GAB*) and for the quantification and stoichiometric analysis of purified transducin, a heterotrimeric G-protein complex. Both GAB-GFP and GAB* contain concatenated sequences of specific proteotypic peptides that are derived from the alpha, beta, and gamma protein subunits of transducin and that are each flanked by spacer regions that maintain the native proteolytic properties for these peptide fragments. Spectroscopic quantification of GAB-GFP provided a molar scale for mass spectrometric ratios from tryptic peptides of GAB* and defined molar responses for mass spectrometric signal intensities from a purified transducin complex. The stoichiometry of transducin subunits alpha, beta, and gamma was measured to be 1:1.1:1.15 over a 5-fold range of labeled internal standard with a relative standard deviation of 9%. Fusing a unique genetically coded spectroscopic signal element with concatenated proteotypic peptides provides a powerful method to accurately quantify and determine the relative stoichiometry of multiple proteins present in complexes or mixtures that cannot be readily assessed using classical gravimetric, enzymatic, or antibody-based technologies.     

Extended polypeptide linkers establish the spatial architecture of a pyruvate dehydrogenase multienzyme complex

Icosahedral pyruvate dehydrogenase (PDH) enzyme complexes are molecular machines consisting of a central E2 core decorated by a shell of peripheral enzymes (E1 and E3) found localized at a distance of approximately 75-90 A from the core. Using a combination of biochemical, biophysical, and cryo-electron microscopic techniques, we show here that the gap between the E2 core and the shell of peripheral enzymes is maintained by the flexible but extended conformation adopted by 60 linker polypeptides that radiate outwards from the inner E2 core, irrespective of the E1 or E3 occupancy. The constancy of the gap is thus not due to protein-protein interactions in the outer protein shell. The extended nature of the E2 inner-linker regions thereby creates the restricted annular space in which the lipoyl domains of E2 that carry catalytic intermediates shuttle between E1, E2, and E3 active sites, while their conformational flexibility facilitates productive encounters.