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.     

A genetic linkage map of Theobroma cacao L.

A linkage map of the cocoa genome comprising 193 loci has been constructed. These loci consist of 5 isozymes, 101 cDNA/RFLPs, 4 loci from genes of known function, 55 genomic DNA/RFLPs and 28 RAPDs. A population of 100 individuals derived from a cross between two heterozygous genotypes was used. Segregation analyses were performed with the JoinMap program. Ten linkage groups, which putatively correspond to the ten gametic chromosomes of cocoa, were identified. The map covers a total length of 759 cM with a 3.9 cM average distance between 2 markers. A small fraction (9%) of the markers deviated significantly from the expected Mendelian ratios.     

An improved synthesis of 1-methylcyclopropanol using the Kulinkovich reaction

An improved process for the preparation of 1-methylcyclopropanol using the Kulinkovich reaction is described. The use of titanium tetramethoxide as catalyst resulted in minimal side product formation. The reaction, isolation and purification procedures were optimized so they can be easily implemented in multi-purpose equipment.     

Orthosteric Binding of ρ-Da1a,a Natural Peptide of Snake Venom Interacting Selectively with the α1A-Adrenoceptor

ρ-Da1a is a three-finger fold toxin from green mamba venom that is highly selective for the α_1A-adrenoceptor. This toxin has atypical pharmacological properties, including incomplete inhibition of ³H-prazosin or ¹²⁵I-HEAT binding and insurmountable antagonist action. We aimed to clarify its mode of action at the α_1A-adrenoceptor. The affinity (pKi 9.26) and selectivity of ρ-Da1a for the α_1A-adrenoceptor were confirmed by comparing binding to human adrenoceptors expressed in eukaryotic cells. Equilibrium and kinetic binding experiments were used to demonstrate that ρ-Da1a, prazosin and HEAT compete at the α_1A-adrenoceptor. ρ-Da1a did not affect the dissociation kinetics of ³H-prazosin or ¹²⁵I-HEAT, and the IC₅₀ of ρ-Da1a, determined by competition experiments, increased linearly with the concentration of radioligands used, while the residual binding by ρ-Da1a remained stable. The effect of ρ-Da1a on agonist-stimulated Ca²⁺ release was insurmountable in the presence of phenethylamine- or imidazoline-type agonists. Ten mutations in the orthosteric binding pocket of the α_1A-adrenoceptor were evaluated for alterations in ρ-Da1a affinity. The D106^3.32A and the S188^5.42A/S192^5.46A receptor mutations reduced toxin affinity moderately (6 and 7.6 times, respectively), while the F86^2.64A, F288^6.51A and F312^7.39A mutations diminished it dramatically by 18- to 93-fold. In addition, residue F86^2.64 was identified as a key interaction point for ¹²⁵I-HEAT, as the variant F86^2.64A induced a 23-fold reduction in HEAT affinity. Unlike the M1 muscarinic acetylcholine receptor toxin MT7, ρ-Da1a interacts with the human α_1A-adrenoceptor orthosteric pocket and shares receptor interaction points with antagonist (F86^2.64, F288^6.51 and F312^7.39) and agonist (F288^6.51 and F312^7.39) ligands. Its selectivity for the α_1A-adrenoceptor may result, at least partly, from its interaction with the residue F86^2.64, which appears to be important also for HEAT binding.     

Phenome-Wide Association Studies on a Quantitative Trait: Application to TPMT Enzyme Activity and Thiopurine Therapy in Pharmacogenomics

Phenome-Wide Association Studies (PheWAS) investigate whether genetic polymorphisms associated with a phenotype are also associated with other diagnoses. In this study, we have developed new methods to perform a PheWAS based on ICD-10 codes and biological test results, and to use a quantitative trait as the selection criterion. We tested our approach on thiopurine S-methyltransferase (TPMT) activity in patients treated by thiopurine drugs. We developed 2 aggregation methods for the ICD-10 codes: an ICD-10 hierarchy and a mapping to existing ICD-9-CM based PheWAS codes. Eleven biological test results were also analyzed using discretization algorithms. We applied these methods in patients having a TPMT activity assessment from the clinical data warehouse of a French academic hospital between January 2000 and July 2013. Data after initiation of thiopurine treatment were analyzed and patient groups were compared according to their TPMT activity level. A total of 442 patient records were analyzed representing 10,252 ICD-10 codes and 72,711 biological test results. The results from the ICD-9-CM based PheWAS codes and ICD-10 hierarchy codes were concordant. Cross-validation with the biological test results allowed us to validate the ICD phenotypes. Iron-deficiency anemia and diabetes mellitus were associated with a very high TPMT activity (p = 0.0004 and p = 0.0015, respectively). We describe here an original method to perform PheWAS on a quantitative trait using both ICD-10 diagnosis codes and biological test results to identify associated phenotypes. In the field of pharmacogenomics, PheWAS allow for the identification of new subgroups of patients who require personalized clinical and therapeutic management.     

Non-Invasive Determination of Left Ventricular Workload in Patients with Aortic Stenosis Using Magnetic Resonance Imaging and Doppler Echocardiography

Early detection and accurate estimation of aortic stenosis (AS) severity are the most important predictors of successful long-term outcomes in patients. Current clinical parameters used for evaluation of the AS severity have several limitations including flow dependency. Estimation of AS severity is specifically challenging in patients with low-flow and low transvalvular pressure gradient conditions. A proper diagnosis in these patients needs a comprehensive evaluation of the left ventricle (LV) hemodynamic loads. This study has two objectives: (1) developing a lumped-parameter model to describe the ventricular-valvular-arterial interaction and to estimate the LV stroke work (SW); (2) introducing and validating a new index, the normalized stroke work (N-SW), to assess the global hemodynamic load imposed on the LV. N-SW represents the global hemodynamic load that the LV faces for each unit volume of blood ejected. The model uses a limited number of parameters which all can be measured non-invasively using current clinical imaging modalities. The model was first validated by comparing its calculated flow waveforms with the ones measured using Cardiovascular Magnetic Resonance (CMR) in 49 patients and 8 controls. A very good correlation and concordance were found throughout the cycle (median root mean square: 12.21 mL/s) and between the peak values (r = 0.98; SEE = 0.001, p<0.001). The model was then used to determine SW using the parameters measured with transthoracic Doppler-echocardiography (TTE) and CMR. N-SW showed very good correlations with a previously-validated index of global hemodynamic load, the valvular arterial impedance (), using data from both imaging modalities (TTE: r = 0.82, SEE = 0.01, p<0.001; CMR: r = 0.74, SEE = 0.01, p<0.001). Furthermore, unlike , N-SW was almost independent from variations in the flow rate. This study suggests that considering N-SW may provide incremental diagnostic and prognostic information, beyond what standard indices of stenosis severity and provide, particularly in patients with low LV outflow.     

New branched amino acids for high affinity dendrimeric DC-SIGN ligands

A branched amino acid was synthesized from methyl glucopyranoside; this amino acid presents three amino groups protected by Fmoc and one acid group and can be used in classic peptide synthesis. In parallel, similar azido terminated blocks were synthesized.Successive coupling reaction and deprotection afforded dendrimers with up to 27 azido functional groups. As an example of application, d-mannose and l-fucose residues were linked through CuAAC coupling and resulting glycodendrimers were evaluated in their interaction with DC-SIGN using SPR competition assay.