TRGT-denovo: accurate detection of de novo tandem repeat mutations.

Motivation: Identifying de novo tandem repeat (TR) mutations on a genome-wide scale is essential for understanding genetic variability and its implications in rare diseases. While PacBio HiFi sequencing data enhances the accessibility of the genome's TR regions for genotyping, simple de novo calling strategies often generate an excess of likely false positives, which can obscure true positive findings, particularly as the number of surveyed genomic regions increases. Results: We developed TRGT-denovo, a computational method designed to accurately identify all types of de novo TR mutations-including expansions, contractions, and compositional changes- within family trios. TRGT-denovo directly interrogates read evidence, allowing for the detection of subtle variations often overlooked in variant call format (VCF) files. TRGT-denovo improves the precision and specificity of de novo mutation (DNM) identification, reducing the number of de novo candidates by an order of magnitude compared to genotype-based approaches. In our experiments involving eight rare disease trios previously studied TRGT-denovo correctly reclassified all false positive DNM candidates as true negatives. Using an expanded repeat catalog, it identified new candidates, of which 95% (19/20) were experimentally validated, demonstrating its effectiveness in minimizing likely false positives while maintaining high sensitivity for true discoveries. Availability and implementation: Built in Rust, TRGT-denovo is available as source code and a pre-compiled Linux binary along with a user guide at: https://github.com/PacificBiosciences/trgt-denovo .

Identification of allele-specific KIV-2 repeats and impact on Lp(a) measurements for cardiovascular disease risk.

The abundance of Lp(a) protein holds significant implications for the risk of cardiovascular disease (CVD), which is directly impacted by the copy number (CN) of KIV-2, a 5.5 kbp sub-region. KIV-2 is highly polymorphic in the population and accurate analysis is challenging. In this study, we present the DRAGEN KIV-2 CN caller, which utilizes short reads. Data across 166 WGS show that the caller has high accuracy, compared to optical mapping and can further phase ~50% of the samples. We compared KIV-2 CN numbers to 24 previously postulated KIV-2 relevant SNVs, revealing that many are ineffective predictors of KIV-2 copy number. Population studies, including USA-based cohorts, showed distinct KIV-2 CN, distributions for European-, African-, and Hispanic-American populations and further underscored the limitations of SNV predictors. We demonstrate that the CN estimates correlate significantly with the available Lp(a) protein levels and that phasing is highly important.

Lower-than-expected linkage disequilibrium between tightly linked markers in humans suggests a role for gene conversion.

Understanding the pattern of linkage disequilibrium (LD) in the human genome is important both for successful implementation of disease-gene mapping approaches and for inferences about human demographic histories. Previous studies have examined LD between loci within single genes or confined genomic regions, which may not be representative of the genome; between loci separated by large distances, where little LD is seen; or in population groups that differ from one study to the next. We measured LD in a large set of locus pairs distributed throughout the genome, with loci within each pair separated by short distances (average 124 bp). Given current models of the history of the human population, nearly all pairs of loci at such short distances would be expected to show complete LD as a consequence of lack of recombination in the short interval. Contrary to this expectation, a significant fraction of pairs showed incomplete LD. A standard model of recombination applied to these data leads to an estimate of effective human population size of 110,000. This estimate is an order of magnitude higher than most estimates based on nucleotide diversity. The most likely explanation of this discrepancy is that gene conversion increases the apparent rate of recombination between nearby loci.

Sequence variation and linkage disequilibrium in the human T-cell receptor beta (TCRB) locus.

The T-cell receptor (TCR) plays a central role in the immune system, and > 90% of human T cells present a receptor that consists of the alpha TCR subunit (TCRA) and the beta subunit (TCRB). Here we report an analysis of 63 variable genes (BV), spanning 553 kb of TCRB that yielded 279 single-nucleotide polymorphisms (SNPs). Samples were drawn from 10 individuals and represent four populations-African American, Chinese, Mexican, and Northern European. We found nine variants that produce nonfunctional BV segments, removing those genes from the TCRB genomic repertoire. There was significant heterogeneity among population samples in SNP frequency (including the BV-inactivating sites), indicating the need for multiple-population samples for adequate variant discovery. In addition, we observed considerable linkage disequilibrium (LD) (r(2) > 0.1) over distances of approximately 30 kb in TCRB, and, in general, the distribution of r(2) as a function of physical distance was in close agreement with neutral coalescent simulations. LD in TCRB showed considerable spatial variation across the locus, being concentrated in "blocks" of LD; however, coalescent simulations of the locus illustrated that the heterogeneity of LD we observed in TCRB did not differ markedly from that expected from neutral processes. Finally, examination of the extended genotypes for each subject demonstrated homozygous stretches of >100 kb in the locus of several individuals. These results provide the basis for optimization of locuswide SNP typing in TCRB for studies of genotype-phenotype association.

An update on the science and therapy of obesity and its relationship to osteoarthritis.

Obesity and osteoarthritis are two commonly encountered clinical problems that can lead to significant physical and emotional disability. This report examines the association between obesity and osteoarthritis, and discusses potential mechanisms by which obesity influences osteoarthritis. Special attention is devoted to reviewing the molecular and genetic mechanisms that underlie the development of clinical obesity. Improved understanding of obesity will hopefully lead to improved treatment and subsequent amelioration of this important risk factor for osteoarthritis.

An analysis of strategies for discovery of single-nucleotide polymorphisms.

Strategies for the discovery of single-nucleotide polymorphisms (SNPs) can be characterized by the number of individuals in the discovery sample, and by the minimal required number of observations of each allele. We examine the effect of different strategies on two key properties of the resulting SNP collection: (1) the probability that a SNP with a given population allele frequency is detected; and (2) the allele-frequency distribution of the discovered SNPs. We show that strategies that accept all polymorphic sites lead to collections with a high fraction of SNPs with rare minor alleles, particularly in expanded populations. Such SNPs have a low probability of replication in a second sample. We discuss how to tailor a discovery strategy to the desired properties of a SNP collection.

Acute reversible hypoxemia in systemic lupus erythematosus.

OBJECTIVE: To determine the frequency of unexplained reversible hypoxemia in patients with systemic lupus erythematosus and to assess the relation between hypoxemia and elevated plasma levels of complement split products. DESIGN: Cohort study. SETTING: Inpatient and outpatient facilities of the New York University Medical Center/Bellevue Hospital and the Hospital for Joint Diseases. PATIENTS: Case patients were 22 patients hospitalized with disease exacerbation and no evidence of parenchymal lung disease on chest roentgenogram. Four patients with stable disease were followed in the outpatient clinic, and five healthy normal volunteers served as controls. MEASUREMENTS: Plasma levels of complement split products (C3a, factor Bb fragment), alveolar-arterial (A-a) Po2 gradients, and pulmonary function were measured. MAIN RESULTS: Nine episodes of hypoxemia or hypocapnia (mean A-a gradient, 30.4 +/- 4.8 mm Hg) or both (despite normal chest roentgenogram results) were noted in six hospitalized patients (group 1). Gas exchange improved within 72 hours of steroid therapy (mean A-a gradient, 11.6 +/- 4.3 mm Hg; P less than 0.01). These patients had an elevated initial mean C3a level (938.4 +/- 246.8 ng/mL) that decreased within 72 hours (407.8 +/- 80.9 ng/mL; P less than 0.01), concomitant with improved oxygenation. Ventilation-perfusion scans, obtained for four of six group 1 patients, excluded pulmonary emboli. Four hospitalized patients (group 2) had a normal A-a gradient (mean, 7.5 +/- 2.7 mm Hg). The mean C3a level of this group (358.3 +/- 39.2 ng/mL) was lower than that of group 1 (P less than 0.05). Four patients with stable disease (group 3) had a mean A-a gradient and a mean C3a level of 3.3 +/- 2.7 mm Hg and 237.8 +/- 105.7 ng/mL, respectively, similar to values found in five normal volunteers, in whom the mean A-a gradient was 3.7 +/- 1.7 mm Hg and the mean C3a level was 124.8 +/- 9.2 ng/mL. CONCLUSION: A syndrome of reversible hypoxemia, unassociated with parenchymal lung disease, is unexpectedly common in acutely ill, hospitalized patients with systemic lupus erythematosus. The pathogenesis of this syndrome is unclear, although the data are compatible with the hypothesis that hypoxemia may be related to pulmonary leukoaggregation.

Methotrexate inhibits neutrophil function by stimulating adenosine release from connective tissue cells.

Although commonly used to control a variety of inflammatory diseases, the mechanism of action of a low dose of methotrexate remains a mystery. Methotrexate accumulates intracellularly where it may interfere with purine metabolism. Therefore, we determined whether a 48-hr pretreatment with methotrexate affected adenosine release from [14C]adenine-labeled human fibroblasts and umbilical vein endothelial cells. Methotrexate significantly increased adenosine release by fibroblasts from 4 +/- 1% to 31 +/- 6% of total purine released (EC50, 1 nM) and by endothelial cells from 24 +/- 4% to 42 +/- 7%. Methotrexate-enhanced adenosine release from fibroblasts was further increased to 51 +/- 4% (EC50, 6 nM) and from endothelial cells was increased to 58 +/- 5% of total purine released by exposure to stimulated (fMet-Leu-Phe at 0.1 microM) neutrophils. The effect of methotrexate on adenosine release was not due to cytotoxicity since cells treated with maximal concentrations of methotrexate took up [14C]adenine and released 14C-labeled purine (a measure of cell injury) in a manner identical to control cells. Methotrexate treatment of fibroblasts dramatically inhibited adherence to fibroblasts by both unstimulated neutrophils (IC50, 9 nM) and stimulated neutrophils (IC50, 13 nM). Methotrexate treatment inhibited neutrophil adherence by enhancing adenosine release from fibroblasts since digestion of extracellular adenosine by added adenosine deaminase completely abrogated the effect of methotrexate on neutrophil adherence without, itself, affecting adherence. One hypothesis that explains the effect of methotrexate on adenosine release is that, by inhibition of 5-aminoimidazole-4-carboxamide ribonucleotide (AICAR) transformylase, methotrexate induces the accumulation of AICAR, the nucleoside precursor of which (5-aminoimidazole-4-carboxamide ribonucleoside referred to hereafter as acadesine) has previously been shown to cause adenosine release from ischemic cardiac tissue. We found that acadesine also promotes adenosine release from and inhibits neutrophil adherence to connective tissue cells. The observation that the antiinflammatory actions of methotrexate are due to the capacity of methotrexate to induce adenosine release may form the basis for the development of an additional class of antiinflammatory drugs.