A saturated map of common genetic variants associated with human height.

Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes1. Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel2) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries.

The power of genetic diversity in genome-wide association studies of lipids.

Increased blood lipid levels are heritable risk factors of cardiovascular disease with varied prevalence worldwide owing to different dietary patterns and medication use1. Despite advances in prevention and treatment, in particular through reducing low-density lipoprotein cholesterol levels2, heart disease remains the leading cause of death worldwide3. Genome-wideassociation studies (GWAS) of blood lipid levels have led to important biological and clinical insights, as well as new drug targets, for cardiovascular disease. However, most previous GWAS4-23 have been conducted in European ancestry populations and may have missed genetic variants that contribute to lipid-level variation in other ancestry groups. These include differences in allele frequencies, effect sizes and linkage-disequilibrium patterns24. Here we conduct a multi-ancestry, genome-wide genetic discovery meta-analysis of lipid levels in approximately 1.65 million individuals, including 350,000 of non-European ancestries. We quantify the gain in studying non-European ancestries and provide evidence to support the expansion of recruitment of additional ancestries, even with relatively small sample sizes. We find that increasing diversity rather than studying additional individuals of European ancestry results in substantial improvements in fine-mapping functional variants and portability of polygenic prediction (evaluated in approximately 295,000 individuals from 7 ancestry groupings). Modest gains in the number of discovered loci and ancestry-specific variants were also achieved. As GWAS expand emphasis beyond the identification of genes and fundamental biology towards the use of genetic variants for preventive and precision medicine25, we anticipate that increased diversity of participants will lead to more accurate and equitable26 application of polygenic scores in clinical practice.

A multi-layer functional genomic analysis to understand noncoding genetic variation in lipids.

A major challenge of genome-wide association studies (GWASs) is to translate phenotypic associations into biological insights. Here, we integrate a large GWAS on blood lipids involving 1.6 million individuals from five ancestries with a wide array of functional genomic datasets to discover regulatory mechanisms underlying lipid associations. We first prioritize lipid-associated genes with expression quantitative trait locus (eQTL) colocalizations and then add chromatin interaction data to narrow the search for functional genes. Polygenic enrichment analysis across 697 annotations from a host of tissues and cell types confirms the central role of the liver in lipid levels and highlights the selective enrichment of adipose-specific chromatin marks in high-density lipoprotein cholesterol and triglycerides. Overlapping transcription factor (TF) binding sites with lipid-associated loci identifies TFs relevant in lipid biology. In addition, we present an integrative framework to prioritize causal variants at GWAS loci, producing a comprehensive list of candidate causal genes and variants with multiple layers of functional evidence. We highlight two of the prioritized genes, CREBRF and RRBP1, which show convergent evidence across functional datasets supporting their roles in lipid biology.

A fast linkage method for population GWAS cohorts with related individuals.

Linkage analysis, a class of methods for detecting co-segregation of genomic segments and traits in families, was used to map disease-causing genes for decades before genotyping arrays and dense SNP genotyping enabled genome-wide association studies in population samples. Population samples often contain related individuals, but the segregation of alleles within families is rarely used because traditional linkage methods are computationally inefficient for larger datasets. Here, we describe Population Linkage, a novel application of Haseman-Elston regression as a method of moments estimator of variance components and their standard errors. We achieve additional computational efficiency by using modern methods for detection of IBD segments and variance component estimation, efficient preprocessing of input data, and minimizing redundant numerical calculations. We also refined variance component models to account for the biases in population-scale methods for IBD segment detection. We ran Population Linkage on four blood lipid traits in over 70,000 individuals from the HUNT and SardiNIA studies, successfully detecting 25 known genetic signals. One notable linkage signal that appeared in both was for low-density lipoprotein (LDL) cholesterol levels in the region near the gene APOE (LOD = 29.3, variance explained = 4.1%). This is the region where the missense variants rs7412 and rs429358, which together make up the ε2, ε3, and ε4 alleles each account for 2.4% and 0.8% of variation in circulating LDL cholesterol. Our results show the potential for linkage analysis and other large-scale applications of method of moments variance components estimation.

Regulatory variants in TCF7L2 are associated with thoracic aortic aneurysm.

Thoracic aortic aneurysm (TAA) is characterized by dilation of the aortic root or ascending/descending aorta. TAA is a heritable disease that can be potentially life threatening. While 10%-20% of TAA cases are caused by rare, pathogenic variants in single genes, the origin of the majority of TAA cases remains unknown. A previous study implicated common variants in FBN1 with TAA disease risk. Here, we report a genome-wide scan of 1,351 TAA-affected individuals and 18,295 control individuals from the Cardiovascular Health Improvement Project and Michigan Genomics Initiative at the University of Michigan. We identified a genome-wide significant association with TAA for variants within the third intron of TCF7L2 following replication with meta-analysis of four additional independent cohorts. Common variants in this locus are the strongest known genetic risk factor for type 2 diabetes. Although evidence indicates the presence of different causal variants for TAA and type 2 diabetes at this locus, we observed an opposite direction of effect. The genetic association for TAA colocalizes with an aortic eQTL of TCF7L2, suggesting a functional relationship. These analyses predict an association of higher expression of TCF7L2 with TAA disease risk. In vitro, we show that upregulation of TCF7L2 is associated with BCL2 repression promoting vascular smooth muscle cell apoptosis, a key driver of TAA disease.

Identification of cell type specific ACE2 modifiers by CRISPR screening.

SARS-CoV-2 infection is initiated by binding of the viral spike protein to its receptor, ACE2, on the surface of host cells. ACE2 expression is heterogeneous both in vivo and in immortalized cell lines, but the molecular pathways that govern ACE2 expression remain unclear. We now report high-throughput CRISPR screens for functional modifiers of ACE2 surface abundance. In liver-derived HuH7 cells, we identified 35 genes whose disruption was associated with a change in the surface abundance of ACE2. Enriched among these ACE2 regulators were established transcription factors, epigenetic regulators, and functional networks. We further characterized individual HuH7 cell lines with disruption of SMAD4, EP300, PIAS1, or BAMBI and found these genes to regulate ACE2 at the mRNA level and to influence cellular susceptibility to SARS-CoV-2 infection. Orthogonal screening of lung-derived Calu-3 cells revealed a distinct set of ACE2 modifiers comprised of ACE2, KDM6A, MOGS, GPAA1, and UGP2. Collectively, our findings clarify the host factors involved in SARS-CoV-2 entry, highlight the cell type specificity of ACE2 regulatory networks, and suggest potential targets for therapeutic development.

Genome-scale CRISPR screening for modifiers of cellular LDL uptake.

Hypercholesterolemia is a causal and modifiable risk factor for atherosclerotic cardiovascular disease. A critical pathway regulating cholesterol homeostasis involves the receptor-mediated endocytosis of low-density lipoproteins into hepatocytes, mediated by the LDL receptor. We applied genome-scale CRISPR screening to query the genetic determinants of cellular LDL uptake in HuH7 cells cultured under either lipoprotein-rich or lipoprotein-starved conditions. Candidate LDL uptake regulators were validated through the synthesis and secondary screening of a customized library of gRNA at greater depth of coverage. This secondary screen yielded significantly improved performance relative to the primary genome-wide screen, with better discrimination of internal positive controls, no identification of negative controls, and improved concordance between screen hits at both the gene and gRNA level. We then applied our customized gRNA library to orthogonal screens that tested for the specificity of each candidate regulator for LDL versus transferrin endocytosis, the presence or absence of genetic epistasis with LDLR deletion, the impact of each perturbation on LDLR expression and trafficking, and the generalizability of LDL uptake modifiers across multiple cell types. These findings identified several previously unrecognized genes with putative roles in LDL uptake and suggest mechanisms for their functional interaction with LDLR.

Genome-wide association study of cardiac troponin I in the general population.

Circulating cardiac troponin proteins are associated with structural heart disease and predict incident cardiovascular disease in the general population. However, the genetic contribution to cardiac troponin I (cTnI) concentrations and its causal effect on cardiovascular phenotypes are unclear. We combine data from two large population-based studies, the Trøndelag Health Study and the Generation Scotland Scottish Family Health Study, and perform a genome-wide association study of high-sensitivity cTnI concentrations with 48 115 individuals. We further use two-sample Mendelian randomization to investigate the causal effects of circulating cTnI on acute myocardial infarction (AMI) and heart failure (HF). We identified 12 genetic loci (8 novel) associated with cTnI concentrations. Associated protein-altering variants highlighted putative functional genes: CAND2, HABP2, ANO5, APOH, FHOD3, TNFAIP2, KLKB1 and LMAN1. Phenome-wide association tests in 1688 phecodes and 83 continuous traits in UK Biobank showed associations between a genetic risk score for cTnI and cardiac arrhythmias, metabolic and anthropometric measures. Using two-sample Mendelian randomization, we confirmed the non-causal role of cTnI in AMI (5948 cases, 355 246 controls). We found indications for a causal role of cTnI in HF (47 309 cases and 930 014 controls), but this was not supported by secondary analyses using left ventricular mass as outcome (18 257 individuals). Our findings clarify the biology underlying the heritable contribution to circulating cTnI and support cTnI as a non-causal biomarker for AMI in the general population. Using genetically informed methods for causal inference helps inform the role and value of measuring cTnI in the general population.

A Systems Approach for Assessing Low Back Pain Care Quality in Veterans Health Administration Chiropractic Visits: A Cross-Sectional Analysis.

OBJECTIVE: The purpose of this study was to explore a systemwide process for assessing components of low back pain (LBP) care quality in Veterans Health Administration (VHA) chiropractic visits using electronic health record (EHR) data. METHODS: We performed a cross-sectional quality improvement project. We randomly sampled 1000 on-station VHA chiropractic initial visits occurring from October 1, 2017, to September 30, 2018, for patients with no such visits within the prior 12 months. Characteristics of LBP visits were extracted from VHA national EHR data via structured data queries and manual chart review. We developed quality indicators for history and/or examination and treatment procedures using previously published literature and calculated frequencies of visits meeting these indicators. Visits meeting our history and/or examination and treatment indicators were classified as "high-quality" visits. We performed a regression analysis to assess associations between demographic/clinical characteristics and visits meeting our quality criteria. RESULTS: There were 592 LBP visits identified. Medical history, physical examination, and neurologic examination were documented in 76%, 77%, and 63% of all LBP visits, respectively. Recommended treatments, such as any manipulation, disease-specific education/advice, and therapeutic exercise, occurred in 75%, 69%, and 40% of chronic visits (n = 383), respectively. In acute/subacute visits (n = 37), any manipulation (92%), manual soft tissue therapy (57%), and disease-specific advice/education (54%) occurred most frequently. Female patients and those with a neck pain comorbid diagnosis were significantly less likely to have a "high-quality" visit, while other regression associations were non-significant. CONCLUSION: This study explored a systemwide process for assessing components of care quality in VHA chiropractic visits for LBP. These results produced a potential framework for uniform assessment of care quality in VHA chiropractic visits for LBP and highlight potential areas for improvements in LBP care quality assessments.

Genetic loci and prioritization of genes for kidney function decline derived from a meta-analysis of 62 longitudinal genome-wide association studies.

Estimated glomerular filtration rate (eGFR) reflects kidney function. Progressive eGFR-decline can lead to kidney failure, necessitating dialysis or transplantation. Hundreds of loci from genome-wide association studies (GWAS) for eGFR help explain population cross section variability. Since the contribution of these or other loci to eGFR-decline remains largely unknown, we derived GWAS for annual eGFR-decline and meta-analyzed 62 longitudinal studies with eGFR assessed twice over time in all 343,339 individuals and in high-risk groups. We also explored different covariate adjustment. Twelve genome-wide significant independent variants for eGFR-decline unadjusted or adjusted for eGFR-baseline (11 novel, one known for this phenotype), including nine variants robustly associated across models were identified. All loci for eGFR-decline were known for cross-sectional eGFR and thus distinguished a subgroup of eGFR loci. Seven of the nine variants showed variant-by-age interaction on eGFR cross section (further about 350,000 individuals), which linked genetic associations for eGFR-decline with age-dependency of genetic cross-section associations. Clinically important were two to four-fold greater genetic effects on eGFR-decline in high-risk subgroups. Five variants associated also with chronic kidney disease progression mapped to genes with functional in-silico evidence (UMOD, SPATA7, GALNTL5, TPPP). An unfavorable versus favorable nine-variant genetic profile showed increased risk odds ratios of 1.35 for kidney failure (95% confidence intervals 1.03-1.77) and 1.27 for acute kidney injury (95% confidence intervals 1.08-1.50) in over 2000 cases each, with matched controls). Thus, we provide a large data resource, genetic loci, and prioritized genes for kidney function decline, which help inform drug development pipelines revealing important insights into the age-dependency of kidney function genetics.

Copy Number Variant Analysis and Genome-wide Association Study Identify Loci with Large Effect for Vesicoureteral Reflux.

BACKGROUND: Vesicoureteral reflux (VUR) is a common, familial genitourinary disorder, and a major cause of pediatric urinary tract infection (UTI) and kidney failure. The genetic basis of VUR is not well understood. METHODS: A diagnostic analysis sought rare, pathogenic copy number variant (CNV) disorders among 1737 patients with VUR. A GWAS was performed in 1395 patients and 5366 controls, of European ancestry. RESULTS: Altogether, 3% of VUR patients harbored an undiagnosed rare CNV disorder, such as the 1q21.1, 16p11.2, 22q11.21, and triple X syndromes ((OR, 3.12; 95% CI, 2.10 to 4.54; P=6.35×10-8) The GWAS identified three study-wide significant and five suggestive loci with large effects (ORs, 1.41-6.9), containing canonical developmental genes expressed in the developing urinary tract (WDPCP, OTX1, BMP5, VANGL1, and WNT5A). In particular, 3.3% of VUR patients were homozygous for an intronic variant in WDPCP (rs13013890; OR, 3.65; 95% CI, 2.39 to 5.56; P=1.86×10-9). This locus was associated with multiple genitourinary phenotypes in the UK Biobank and eMERGE studies. Analysis of Wnt5a mutant mice confirmed the role of Wnt5a signaling in bladder and ureteric morphogenesis. CONCLUSIONS: These data demonstrate the genetic heterogeneity of VUR. Altogether, 6% of patients with VUR harbored a rare CNV or a common variant genotype conferring an OR >3. Identification of these genetic risk factors has multiple implications for clinical care and for analysis of outcomes in VUR.

Genome-wide Study of Atrial Fibrillation Identifies Seven Risk Loci and Highlights Biological Pathways and Regulatory Elements Involved in Cardiac Development.

Atrial fibrillation (AF) is a common cardiac arrhythmia and a major risk factor for stroke, heart failure, and premature death. The pathogenesis of AF remains poorly understood, which contributes to the current lack of highly effective treatments. To understand the genetic variation and biology underlying AF, we undertook a genome-wide association study (GWAS) of 6,337 AF individuals and 61,607 AF-free individuals from Norway, including replication in an additional 30,679 AF individuals and 278,895 AF-free individuals. Through genotyping and dense imputation mapping from whole-genome sequencing, we tested almost nine million genetic variants across the genome and identified seven risk loci, including two novel loci. One novel locus (lead single-nucleotide variant [SNV] rs12614435; p = 6.76 × 10-18) comprised intronic and several highly correlated missense variants situated in the I-, A-, and M-bands of titin, which is the largest protein in humans and responsible for the passive elasticity of heart and skeletal muscle. The other novel locus (lead SNV rs56202902; p = 1.54 × 10-11) covered a large, gene-dense chromosome 1 region that has previously been linked to cardiac conduction. Pathway and functional enrichment analyses suggested that many AF-associated genetic variants act through a mechanism of impaired muscle cell differentiation and tissue formation during fetal heart development.

Phylogeny, sequence-typing and virulence profile of uropathogenic Escherichia coli (UPEC) strains from Pakistan.

BACKGROUND: Escherichia coli lineage ST131 predominates across various spectra of extra-intestinal infections, including urinary tract infection (UTI). The distinctive resistance profile, diverse armamentarium of virulence factors and rapid global dissemination of ST131 E. coli makes it an intriguing pathogen. However, not much is known about the prevalence and genetic attributes of ST131 lineage in Pakistan. METHODS: We estimated prevalence and genetic attributes of E. coli ST131 isolates causing UTI among 155 randomly selected samples. Samples were analyzed for phylogenetic grouping, O-typing and fumC/fimH typing. Isolates were further tested for the ESBL and virulence factors using PCR. RESULTS: Overall, 59% of the UPEC isolates belonged to the phylogenetic group B2, followed by D = 28%, B1 = 8% and A = 5%. Among 18 different Sequence-types, ST131 was the dominant lineage (n = 71; 46%) out of which 72% of the isolates were assigned to the phylogenetic group B2, while 61% adhered to the serogroup O25b. FumC/fimH typing confirmed 49% of the ST131 as H30 sub-types. In this study, significant numbers of the identified ST131 isolates were MDR and 42% showed ESBL phenotypes, out of which 37% carried bla-CTX-M-15. Moreover, different virulence factors were detected in following percentages: fimH,155(100%), iutA 86 (55%), feoB 76 (49%), papC 75 (48%), papGII 70 (45%), kpsMTII 40 (26%), papEF 37 (24%), fyuA 37 (24%), usp 22 (14%), papA 20 (13%), sfa/foc20 (13%), hlyA 18 (12%), afa 15 (10%), cdtB 11 (7%), papGI 6 (4%), papGIII 6 (4%), kpsMTIII 4 (3%) and bmaE2 (1%). CONCLUSION: Conclusively, this study provides important insight into the genetic and virulence attributes of pandemic MDR ST131 strains involved in UTIs. It also highlights higher prevalence of ST131-O25b-H30 UPEC isolates in patients, which was previously unreported from this part of globe.

Predicting Displaceable Water Sites Using Mixed-Solvent Molecular Dynamics.

Water molecules are an important factor in protein-ligand binding. Upon binding of a ligand with a protein's surface, waters can either be displaced by the ligand or may be conserved and possibly bridge interactions between the protein and ligand. Depending on the specific interactions made by the ligand, displacing waters can yield a gain in binding affinity. The extent to which binding affinity may increase is difficult to predict, as the favorable displacement of a water molecule is dependent on the site-specific interactions made by the water and the potential ligand. Several methods have been developed to predict the location of water sites on a protein's surface, but the majority of methods are not able to take into account both protein dynamics and the interactions made by specific functional groups. Mixed-solvent molecular dynamics (MixMD) is a cosolvent simulation technique that explicitly accounts for the interaction of both water and small molecule probes with a protein's surface, allowing for their direct competition. This method has previously been shown to identify both active and allosteric sites on a protein's surface. Using a test set of eight systems, we have developed a method using MixMD to identify conserved and displaceable water sites. Conserved sites can be determined by an occupancy-based metric to identify sites which are consistently occupied by water even in the presence of probe molecules. Conversely, displaceable water sites can be found by considering the sites which preferentially bind probe molecules. Furthermore, the inclusion of six probe types allows the MixMD method to predict which functional groups are capable of displacing which water sites. The MixMD method consistently identifies sites which are likely to be nondisplaceable and predicts the favorable displacement of water sites that are known to be displaced upon ligand binding.

MixMD Probeview: Robust Binding Site Prediction from Cosolvent Simulations.

Mixed-solvent molecular dynamics (MixMD) is a cosolvent simulation technique for identifying binding hotspots and specific favorable interactions on a protein's surface. MixMD studies have the ability to identify these biologically relevant sites by examining the occupancy of the cosolvent over the course of the simulation. However, previous MixMD analysis required a great deal of manual inspection to identify relevant sites. To address this limitation, we have developed MixMD Probeview as a plugin for the freely available, open-source version of the molecular visualization program PyMOL. MixMD Probeview incorporates two analysis procedures: (1) identifying and ranking whole binding sites and (2) identifying and ranking local maxima for each probe type. These functionalities were validated using four common benchmark proteins, including two with both active and allosteric sites. In addition, three different cosolvent procedures were compared to examine the impact of including more than one cosolvent in the simulations. For all systems tested, MixMD Probeview successfully identified known active and allosteric sites based on the total occupancy of neutral probe molecules. As an easy-to-use PyMOL plugin, we expect that MixMD Probeview will facilitate identification and analysis of binding sites from cosolvent simulations performed on a wide range of systems.

Identifying binding hot spots on protein surfaces by mixed-solvent molecular dynamics: HIV-1 protease as a test case.

Mixed-solvent molecular dynamics (MixMD) simulations use full protein flexibility and competition between water and small organic probes to achieve accurate hot-spot mapping on protein surfaces. In this study, we improved MixMD using human immunodeficiency virus type-1 protease (HIVp) as the test case. We used three probe-water solutions (acetonitrile-water, isopropanol-water, and pyrimidine-water), first at 50% w/w concentration and later at 5% v/v. Paradoxically, better mapping was achieved by using fewer probes; 5% simulations gave a superior signal-to-noise ratio and far fewer spurious hot spots than 50% MixMD. Furthermore, very intense and well-defined probe occupancies were observed in the catalytic site and potential allosteric sites that have been confirmed experimentally. The Eye site, an allosteric site underneath the flap of HIVp, has been confirmed by the presence of a 5-nitroindole fragment in a crystal structure. MixMD also mapped two additional hot spots: the Exo site (between the Gly16-Gly17 and Cys67-Gly68 loops) and the Face site (between Glu21-Ala22 and Val84-Ile85 loops). The Exo site was observed to overlap with crystallographic additives such as acetate and dimethyl sulfoxide that are present in different crystal forms of the protein. Analysis of crystal structures of HIVp in different symmetry groups has shown that some surface sites are common interfaces for crystal contacts, which means that they are surfaces that are relatively easy to desolvate and complement with organic molecules. MixMD should identify these sites; in fact, their occupancy values help establish a solid cut-off where "druggable" sites are required to have higher occupancies than the crystal-packing faces.

Multi-trait analysis characterizes the genetics of thyroid function and identifies causal associations with clinical implications.

To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

Evaluating rates of chiropractic use and utilization by patient sex within the United States Veterans Health Administration: a serial cross-sectional analysis.

BACKGROUND: Within the United States Veterans Health Administration (VHA), the number of patients using healthcare services has increased over the past several decades. Females make up a small proportion of overall patients within the VHA; however, this proportion is growing rapidly. Previous studies have described rates of VHA chiropractic use; however, no prior study assessed differences in use or utilization rates between male and female veterans. The purpose of this study was to assess rates of use and utilization of chiropractic care by sex among VHA patients receiving care at VHA facilities with on-station chiropractic clinics. METHODS: A serial cross-sectional analysis of VHA national electronic health record data was conducted in Fall 2021 for fiscal year (FY) 2005-2021. The cohort population was defined as VHA facilities with on-station chiropractic clinics, and facilities were admitted to the cohort after the first FY with a minimum of 500 on-station chiropractic visits. Variables extracted included counts of unique users of any VHA on-station facility outpatient services, unique users of VHA on-station facility chiropractic services, number of chiropractic visits, and sex. To calculate use, we determined the proportion of patients of each sex who received chiropractic services to the total patients of the same sex receiving any outpatient care within each facility. To calculate utilization, we determined the number of chiropractic care visits per patient per fiscal year. A linear mixed effects model was applied to examine the difference in chiropractic care utilization by sex. RESULTS: The percentage of female VHA on-station chiropractic patients increased from 11.7 to 17.7% from FY2005-FY2021. Among VHA facilities with on-station chiropractic care, the percentage of female VHA healthcare users who used chiropractic care (mean = 2.3%) was greater than the percentage of male VHA healthcare users who used chiropractic care (mean = 1.1%). Rates of chiropractic utilization by sex among VHA facilities with on-station chiropractic clinics were slightly higher for females (median = 4.3 visits per year, mean = 4.9) compared to males (median = 4.1 visits per year, mean = 4.6). CONCLUSION: We report higher use and utilization of VHA chiropractic care by females compared with males, yet for both sexes rates were lower than in the private US healthcare system. This highlights the need for further assessment of the determinants and outcomes of VHA chiropractic care.

Genome-wide association study of thoracic aortic aneurysm and dissection in the Million Veteran Program.

The current understanding of the genetic determinants of thoracic aortic aneurysms and dissections (TAAD) has largely been informed through studies of rare, Mendelian forms of disease. Here, we conducted a genome-wide association study (GWAS) of TAAD, testing ~25 million DNA sequence variants in 8,626 participants with and 453,043 participants without TAAD in the Million Veteran Program, with replication in an independent sample of 4,459 individuals with and 512,463 without TAAD from six cohorts. We identified 21 TAAD risk loci, 17 of which have not been previously reported. We leverage multiple downstream analytic methods to identify causal TAAD risk genes and cell types and provide human genetic evidence that TAAD is a non-atherosclerotic aortic disorder distinct from other forms of vascular disease. Our results demonstrate that the genetic architecture of TAAD mirrors that of other complex traits and that it is not solely inherited through protein-altering variants of large effect size.