Genetic diversity fuels gene discovery for tobacco and alcohol use.

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury1-4. These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries5. Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction.

An intact S-layer is advantageous to Clostridioides difficile within the host.

Clostridioides difficile is responsible for substantial morbidity and mortality in antibiotically-treated, hospitalised, elderly patients, in which toxin production correlates with diarrhoeal disease. While the function of these toxins has been studied in detail, the contribution of other factors, including the paracrystalline surface layer (S-layer), to disease is less well understood. Here, we highlight the essentiality of the S-layer in vivo by reporting the recovery of S-layer variants, following infection with the S-layer-null strain, FM2.5. These variants carry either correction of the original point mutation, or sequence modifications which restored the reading frame, and translation of slpA. Selection of these variant clones was rapid in vivo, and independent of toxin production, with up to 90% of the recovered C. difficile population encoding modified slpA sequence within 24 h post infection. Two variants, subsequently named FM2.5varA and FM2.5varB, were selected for study in greater detail. Structural determination of SlpA from FM2.5varB indicated an alteration in the orientation of protein domains, resulting in a reorganisation of the lattice assembly, and changes in interacting interfaces, which might alter function. Interestingly, variant FM2.5varB displayed an attenuated, FM2.5-like phenotype in vivo compared to FM2.5varA, which caused disease severity more comparable to that of R20291. Comparative RNA sequencing (RNA-Seq) analysis of in vitro grown isolates revealed large changes in gene expression between R20291 and FM2.5. Downregulation of tcdA/tcdB and several genes associated with sporulation and cell wall integrity may account for the reported attenuated phenotype of FM2.5 in vivo. RNA-seq data correlated well with disease severity with the more virulent variant, FM2.5varA, showing s similar profile of gene expression to R20291 in vitro, while the attenuated FM2.5varB showed downregulation of many of the same virulence associated traits as FM2.5. Cumulatively, these data add to a growing body of evidence that the S-layer contributes to C. difficile pathogenesis and disease severity.

Stratifying macrophages based on their infectious burden identifies novel host targets for intervention during Crohn's disease associated adherent-invasive Escherichia coli infection.

Bacterial infection is a dynamic process resulting in a heterogenous population of infected and uninfected cells. These cells respond differently based on their bacterial load and duration of infection. In the case of infection of macrophages with Crohn's disease (CD) associated adherent-invasive Escherichia coli (AIEC), understanding the drivers of pathogen success may allow targeting of cells where AIEC replicate to high levels. Here we show that stratifying immune cells based on their bacterial load identifies novel pathways and therapeutic targets not previously associated with AIEC when using a traditional homogeneous infected population approach. Using flow cytometry-based cell sorting we stratified cells into those with low or high intracellular pathogen loads, or those which were bystanders to infection. Immune cells transcriptomics revealed a diverse response to the varying levels of infection while pathway analysis identified novel intervention targets that were directly related to increasing intracellular AIEC numbers. Chemical inhibition of identified targets reduced AIEC intracellular replication or inhibited secretion of tumour necrosis factor alpha (TNFα), a key cytokine associated with AIEC infection. Our results have identified new avenues of intervention in AIEC infection that may also be applicable to CD through the repurposing of already available inhibitors. Additionally, they highlight the applicability of immune cell stratification post-infection as an effective approach for the study of microbial pathogens.

Machine learning to identify endometrial biomarkers predictive of pregnancy success following artificial insemination in dairy cows.

The objective was to identify a set of genes whose transcript abundance is predictive of a cow's ability to become pregnant following artificial insemination (AI). Endometrial epithelial cells from the uterine body were collected for RNA sequencing using the cytobrush method from 193 first-service Holstein cows at estrus prior to AI (day 0). A group of 253 first-service cows not used for cytobrush collection were controls. There was no effect of cytobrush collection on pregnancy outcomes at day 30 or 70 or on pregnancy loss between day 30 and 70. There were 2 upregulated and 214 downregulated genes (FDR < 0.05, absolute fold change >2-fold) for cows pregnant at day 30 versus those that were not pregnant. Functional terms overrepresented in the downregulated genes included those related to immune and inflammatory responses. Machine learning for fertility biomarkers with the R package BORUTA resulted in identification of 57 biomarkers that predicted pregnancy outcome at day 30 with an average accuracy of 77%. Thus, machine learning can identify predictive biomarkers of pregnancy in endometrium with high accuracy. Moreover, sampling of endometrial epithelium using the cytobrush can help understand functional characteristics of the endometrium at AI without compromising cow fertility. Functional characteristics of the genes comprising the set of biomarkers is indicative that a major determinant of cow fertility, at least for first insemination after calving, is immune status of the uterus, which, in turn, is likely to reflect the previous history of uterine disease.

Differences in transcriptional changes in psoriasis and psoriatic arthritis skin with immunoglobulin gene enrichment in psoriatic arthritis.

OBJECTIVES: Approximately 20% of people with psoriasis develop PsA. Although genetic, clinical and environmental risk factors have been identified, it is not known why some people with psoriasis develop PsA. The skin disease is traditionally considered the same in both. This study compares transcriptional changes in psoriasis and PsA skin for the first time. METHODS: Skin biopsies were collected from healthy controls (HC), and uninvolved and lesional skin from patients with PsA. Bulk tissue sequencing was performed and analysed using the pipeline Searchlight 2.0. Transcriptional changes in PsA skin were compared with existing sequencing data from participants with psoriasis without PsA (GSE121212). Psoriasis and PsA datasets could not be directly compared as different analysis methods were used. Data from participants with PsA in the GSE121212 dataset were used for validation. RESULTS: Skin samples from 9 participants with PsA and 9 HC were sequenced, analysed and compared with available transcriptomic data for 16 participants with psoriasis compared with 16 HC. Uninvolved skin in psoriasis shared transcriptional changes with lesional skin in psoriasis, but uninvolved skin in PsA did not. Most transcriptional changes in psoriasis and PsA lesional skin were shared, but immunoglobulin genes were upregulated in PsA lesional skin specifically. The transcription factor POU2F1, which regulates immunoglobulin gene expression, was enriched in PsA lesional skin. This was confirmed in the validation cohort. CONCLUSIONS: Immunoglobulin genes are upregulated in PsA but not in psoriasis skin lesions. This may have implications for the spread from the cutaneous compartment to other tissues.

Differential regulation of IL-17A and IL-17F via STAT5 contributes to psoriatic disease.

BACKGROUND: IL-17A plays a pivotal pathogenic role in several immune-mediated inflammatory diseases. Despite sharing 50% sequence homology with IL-17A, the role of IL-17F remains less clear. Clinical findings suggest that dual inhibition of IL-17A and IL-17F in psoriatic disease is more efficacious than IL-17A inhibition alone, positing a pathogenic role for IL-17F. OBJECTIVE: We characterized the regulation of IL-17A and IL-17F in psoriatic disease. METHODS: Using both in vitro systems and lesional skin tissue from patients, we interrogated the chromosomal, transcriptional, and protein expression landscape of IL-17A+ and IL-17F+ TH17 cells. Alongside established assays such as single-cell RNA sequencing, we developed a novel cytokine-capture technique that was combined with chromatin immunoprecipitation sequencing and RNA sequencing. RESULTS: We confirm a preferential elevation of IL-17F over IL-17A in psoriatic disease and show that expression of each isoform predominantly occurs in distinct cell populations. The expression of both IL-17A and IL-17F exhibited a high degree of plasticity, with the balance between the 2 isoforms influenced by proinflammatory signaling and by anti-inflammatory drugs such as methylprednisolone. This plasticity was reflected in a broad H3K4me3 region at the IL17A-F locus, while opposing effects of STAT5/IL-2 signaling were observed for each of the 2 genes. Functionally, higher IL17F expression was linked to greater cell proliferation. CONCLUSION: There are key differences in the regulation of IL-17A and IL-17F in psoriatic disease, leading to distinct inflammatory cell populations. As such, we propose that both IL-17A and IL-17F neutralization may be required to maximally inhibit IL-17-driven pathology.

Novel computational biology modeling system can accurately forecast response to neoadjuvant therapy in early breast cancer.

BACKGROUND: Generalizable population-based studies are unable to account for individual tumor heterogeneity that contributes to variability in a patient's response to physician-chosen therapy. Although molecular characterization of tumors has advanced precision medicine, in early-stage and locally advanced breast cancer patients, predicting a patient's response to neoadjuvant therapy (NAT) remains a gap in current clinical practice. Here, we perform a study in an independent cohort of early-stage and locally advanced breast cancer patients to forecast tumor response to NAT and assess the stability of a previously validated biophysical simulation platform. METHODS: A single-blinded study was performed using a retrospective database from a single institution (9/2014-12/2020). Patients included: ≥ 18 years with breast cancer who completed NAT, with pre-treatment dynamic contrast enhanced magnetic resonance imaging. Demographics, chemotherapy, baseline (pre-treatment) MRI and pathologic data were input into the TumorScope Predict (TS) biophysical simulation platform to generate predictions. Primary outcomes included predictions of pathological complete response (pCR) versus residual disease (RD) and final volume for each tumor. For validation, post-NAT predicted pCR and tumor volumes were compared to actual pathological assessment and MRI-assessed volumes. Predicted pCR was pre-defined as residual tumor volume ≤ 0.01 cm3 (≥ 99.9% reduction). RESULTS: The cohort consisted of eighty patients; 36 Caucasian and 40 African American. Most tumors were high-grade (54.4% grade 3) invasive ductal carcinomas (90.0%). Receptor subtypes included hormone receptor positive (HR+)/human epidermal growth factor receptor 2 positive (HER2+, 30%), HR+/HER2- (35%), HR-/HER2+ (12.5%) and triple negative breast cancer (TNBC, 22.5%). Simulated tumor volume was significantly correlated with post-treatment radiographic MRI calculated volumes (r = 0.53, p = 1.3 × 10-7, mean absolute error of 6.57%). TS prediction of pCR compared favorably to pathological assessment (pCR: TS n = 28; Path n = 27; RD: TS n = 52; Path n = 53), for an overall accuracy of 91.2% (95% CI: 82.8% - 96.4%; Clopper-Pearson interval). Five-year risk of recurrence demonstrated similar prognostic performance between TS predictions (Hazard ratio (HR): - 1.99; 95% CI [- 3.96, - 0.02]; p = 0.043) and clinically assessed pCR (HR: - 1.76; 95% CI [- 3.75, 0.23]; p = 0.054). CONCLUSION: We demonstrated TS ability to simulate and model tumor in vivo conditions in silico and forecast volume response to NAT across breast tumor subtypes.

The relevance of rich club regions for functional outcome post-stroke is enhanced in women.

This study aimed to investigate the influence of stroke lesions in predefined highly interconnected (rich-club) brain regions on functional outcome post-stroke, determine their spatial specificity and explore the effects of biological sex on their relevance. We analyzed MRI data recorded at index stroke and ~3-months modified Rankin Scale (mRS) data from patients with acute ischemic stroke enrolled in the multisite MRI-GENIE study. Spatially normalized structural stroke lesions were parcellated into 108 atlas-defined bilateral (sub)cortical brain regions. Unfavorable outcome (mRS > 2) was modeled in a Bayesian logistic regression framework. Effects of individual brain regions were captured as two compound effects for (i) six bilateral rich club and (ii) all further non-rich club regions. In spatial specificity analyses, we randomized the split into "rich club" and "non-rich club" regions and compared the effect of the actual rich club regions to the distribution of effects from 1000 combinations of six random regions. In sex-specific analyses, we introduced an additional hierarchical level in our model structure to compare male and female-specific rich club effects. A total of 822 patients (age: 64.7[15.0], 39% women) were analyzed. Rich club regions had substantial relevance in explaining unfavorable functional outcome (mean of posterior distribution: 0.08, area under the curve: 0.8). In particular, the rich club-combination had a higher relevance than 98.4% of random constellations. Rich club regions were substantially more important in explaining long-term outcome in women than in men. All in all, lesions in rich club regions were associated with increased odds of unfavorable outcome. These effects were spatially specific and more pronounced in women.

Comprehensive analyses of 723 transcriptomes enhance genetic and biological interpretations for complex traits in cattle.

By uniformly analyzing 723 RNA-seq data from 91 tissues and cell types, we built a comprehensive gene atlas and studied tissue specificity of genes in cattle. We demonstrated that tissue-specific genes significantly reflected the tissue-relevant biology, showing distinct promoter methylation and evolution patterns (e.g., brain-specific genes evolve slowest, whereas testis-specific genes evolve fastest). Through integrative analyses of those tissue-specific genes with large-scale genome-wide association studies, we detected relevant tissues/cell types and candidate genes for 45 economically important traits in cattle, including blood/immune system (e.g., CCDC88C) for male fertility, brain (e.g., TRIM46 and RAB6A) for milk production, and multiple growth-related tissues (e.g., FGF6 and CCND2) for body conformation. We validated these findings by using epigenomic data across major somatic tissues and sperm. Collectively, our findings provided novel insights into the genetic and biological mechanisms underlying complex traits in cattle, and our transcriptome atlas can serve as a primary source for biological interpretation, functional validation, studies of adaptive evolution, and genomic improvement in livestock.

Fibroblast growth factor 9 (FGF9)-mediated neurodegeneration: Implications for progressive multiple sclerosis?

AIMS: Fibroblast growth factor (FGF) signalling is dysregulated in multiple sclerosis (MS) and other neurological and psychiatric conditions, but there is little or no consensus as to how individual FGF family members contribute to disease pathogenesis. Lesion development in MS is associated with increased expression of FGF1, FGF2 and FGF9, all of which modulate remyelination in a variety of experimental settings. However, FGF9 is also selectively upregulated in major depressive disorder (MDD), prompting us to speculate it may also have a direct effect on neuronal function and survival. METHODS: Transcriptional profiling of myelinating cultures treated with FGF1, FGF2 or FGF9 was performed, and the effects of FGF9 on cortical neurons investigated using a combination of transcriptional, electrophysiological and immunofluorescence microscopic techniques. The in vivo effects of FGF9 were explored by stereotactic injection of adeno-associated viral (AAV) vectors encoding either FGF9 or EGFP into the rat motor cortex. RESULTS: Transcriptional profiling of myelinating cultures after FGF9 treatment revealed a distinct neuronal response with a pronounced downregulation of gene networks associated with axonal transport and synaptic function. In cortical neuronal cultures, FGF9 also rapidly downregulated expression of genes associated with synaptic function. This was associated with a complete block in the development of photo-inducible spiking activity, as demonstrated using multi-electrode recordings of channel rhodopsin-transfected rat cortical neurons in vitro and, ultimately, neuronal cell death. Overexpression of FGF9 in vivo resulted in rapid loss of neurons and subsequent development of chronic grey matter lesions with neuroaxonal reduction and ensuing myelin loss. CONCLUSIONS: These observations identify overexpression of FGF9 as a mechanism by which neuroaxonal pathology could develop independently of immune-mediated demyelination in MS. We suggest targeting neuronal FGF9-dependent pathways may provide a novel strategy to slow if not halt neuroaxonal atrophy and loss in MS, MDD and potentially other neurodegenerative diseases.

Functional annotation of regulatory elements in cattle genome reveals the roles of extracellular interaction and dynamic change of chromatin states in rumen development during weaning.

We profiled landscapes of bovine regulatory elements and explored dynamic changes of chromatin states in rumen development during weaning. The regulatory elements (15 chromatin states) and their coordinated activities in cattle were defined through genome-wide profiling of four histone modifications, CTCF-binding, DNA accessibility, DNA methylation, and transcriptome in rumen epithelial tissues. Each chromatin state presented specific enrichment for sequence ontology, methylation, trait-associated variants, transcription, gene expression-associated variants, selection signatures, and evolutionarily conserved elements. During weaning, weak enhancers and flanking active transcriptional start sites (TSS) were the most dynamic chromatin states and occurred in tandem with significant variations in gene expression and DNA methylation, significantly associated with stature, production, and reproduction economic traits. By comparing with in vitro cultured epithelial cells and in vivo rumen tissues, we showed the commonness and uniqueness of these results, especially the roles of cell interactions and mitochondrial activities in tissue development.

Black-White Differences in Ischemic Stroke Risk Factor Burden in Young Adults.

BACKGROUND AND PURPOSE: Although the US Black population has a higher incidence of stroke compared with the US White population, few studies have addressed Black-White differences in the contribution of vascular risk factors to the population burden of ischemic stroke in young adults. METHODS: A population-based case-control study of early-onset ischemic stroke, ages 15 to 49 years, was conducted in the Baltimore-Washington DC region between 1992 and 2007. Risk factor data was obtained by in-person interview in both cases and controls. The prevalence, odds ratio, and population-attributable risk percent (PAR%) of smoking, diabetes, and hypertension was determined among Black patients and White patients, stratified by sex. RESULTS: The study included 1044 cases and 1099 controls. Of the cases, 47% were Black patients, 54% were men, and the mean (±SD) age was 41.0 (±6.8) years. For smoking, the population-attributable risk percent were White men 19.7%, White women 32.5%, Black men 10.1%, and Black women 23.8%. For diabetes, the population-attributable risk percent were White men 10.5%, White women 7.4%, Black men 17.2%, and Black women 13.4%. For hypertension, the population-attributable risk percent were White men 17.2%, White women 19.3%, Black men 45.8%, and Black women 26.4%. CONCLUSIONS: Modifiable vascular risk factors account for a large proportion of ischemic stroke in young adults. Cigarette smoking was the strongest contributor to stroke among White patients while hypertension was the strongest contributor to stroke among Black patients. These results support early primary prevention efforts focused on smoking cessation and hypertension detection and treatment.

Genome-wide recombination map construction from single sperm sequencing in cattle.

BACKGROUND: Meiotic recombination is one of the important phenomena contributing to gamete genome diversity. However, except for human and a few model organisms, it is not well studied in livestock, including cattle. RESULTS: To investigate their distributions in the cattle sperm genome, we sequenced 143 single sperms from two Holstein bulls. We mapped meiotic recombination events at high resolution based on phased heterozygous single nucleotide polymorphism (SNP). In the absence of evolutionary selection pressure in fertilization and survival, recombination events in sperm are enriched near distal chromosomal ends, revealing that such a pattern is intrinsic to the molecular mechanism of meiosis. Furthermore, we further validated these findings in single sperms with results derived from sequencing its family trio of diploid genomes and our previous studies of recombination in cattle. CONCLUSIONS: To our knowledge, this is the first large-scale single sperm whole-genome sequencing effort in livestock, which provided useful information for future studies of recombination, genome instability, and male infertility.

Factors influencing the initiation of adjuvant endocrine therapy in patients with estrogen receptor-positive ductal carcinoma in situ: a single institution experience.

PURPOSE: This study evaluates whether race, socioeconomic status, insurance type, oncological provider type, and prior cancer treatment are associated with the suggestion and acceptance of hormonal therapy in patients with estrogen receptor (ER)-positive Ductal carcinoma in situ (DCIS). This study also assesses whether disparities exist pertaining to prescription of such medications. METHODS: This single-center retrospective study included 111 patients diagnosed with DCIS between 2020 and 2021. Data collected included race, type of insurance, prescribing providers, and socioeconomic status. We used zip codes to identify the poverty levels in these areas as published in the United States Census Bureau and stratified the patients into quartiles accordingly. Chi-Square statistics were used to calculate significance levels. RESULTS: There was no significant correlation between the intake of hormonal therapy and race (p = 0.60), insurance (p = 0.50), socioeconomic (p = 0.58), or providers (p = 0.99). 79.3% of women were offered endocrine therapy. Of those who were offered endocrine therapy, 70.8% accepted. Of patients not on hormonal therapy, 45.8% were not recommended the medications by their provider, and 54.2% declined treatment when offered. CONCLUSION: In this study, patients' demographics and providers were not associated with adjuvant hormonal therapy initiation in DCIS. Our results show that abstaining from endocrine therapy in DCIS patients is both due to lack of provider recommendation and patient rejection of these medications. The wide variation in hormonal therapy treatment among ER-positive DCIS patients suggests a need for improved provider-patient communication regarding the risks and benefits of endocrine therapy in order to ensure a shared decision-making process.

Highly accurate response prediction in high-risk early breast cancer patients using a biophysical simulation platform.

PURPOSE: Pathologic complete response (pCR) to neoadjuvant chemotherapy (NAC) in early breast cancer (EBC) is largely dependent on breast cancer subtype, but no clinical-grade model exists to predict response and guide selection of treatment. A biophysical simulation of response to NAC has the potential to address this unmet need. METHODS: We conducted a retrospective evaluation of a biophysical simulation model as a predictor of pCR. Patients who received standard NAC at the University of Chicago for EBC between January 1st, 2010 and March 31st, 2020 were included. Response was predicted using baseline breast MRI, clinicopathologic features, and treatment regimen by investigators who were blinded to patient outcomes. RESULTS: A total of 144 tumors from 141 patients were included; 59 were triple-negative, 49 HER2-positive, and 36 hormone-receptor positive/HER2 negative. Lymph node disease was present in half of patients, and most were treated with an anthracycline-based regimen (58.3%). Sensitivity and specificity of the biophysical simulation for pCR were 88.0% (95% confidence interval [CI] 75.7 - 95.5) and 89.4% (95% CI 81.3 - 94.8), respectively, with robust results regardless of subtype. In patients with predicted pCR, 5-year event-free survival was 98%, versus 79% with predicted residual disease (log-rank p = 0.01, HR 4.57, 95% CI 1.36 - 15.34). At a median follow-up of 5.4 years, no patients with predicted pCR experienced disease recurrence. CONCLUSION: A biophysical simulation model accurately predicts pCR and long-term outcomes from baseline MRI and clinical data, and is a promising tool to guide escalation/de-escalation of NAC.

Single nucleotide variations in ZBTB46 are associated with post-thrombolytic parenchymal haematoma.

Haemorrhagic transformation is a complication of recombinant tissue-plasminogen activator treatment. The most severe form, parenchymal haematoma, can result in neurological deterioration, disability, and death. Our objective was to identify single nucleotide variations associated with a risk of parenchymal haematoma following thrombolytic therapy in patients with acute ischaemic stroke. A fixed-effect genome-wide meta-analysis was performed combining two-stage genome-wide association studies (n = 1904). The discovery stage (three cohorts) comprised 1324 ischaemic stroke individuals, 5.4% of whom had a parenchymal haematoma. Genetic variants yielding a P-value < 0.05 1 × 10-5 were analysed in the validation stage (six cohorts), formed by 580 ischaemic stroke patients with 12.1% haemorrhagic events. All participants received recombinant tissue-plasminogen activator; cases were parenchymal haematoma type 1 or 2 as defined by the European Cooperative Acute Stroke Study (ECASS) criteria. Genome-wide significant findings (P < 5 × 10-8) were characterized by in silico functional annotation, gene expression, and DNA regulatory elements. We analysed 7 989 272 single nucleotide polymorphisms and identified a genome-wide association locus on chromosome 20 in the discovery cohort; functional annotation indicated that the ZBTB46 gene was driving the association for chromosome 20. The top single nucleotide polymorphism was rs76484331 in the ZBTB46 gene [P = 2.49 × 10-8; odds ratio (OR): 11.21; 95% confidence interval (CI): 4.82-26.55]. In the replication cohort (n = 580), the rs76484331 polymorphism was associated with parenchymal haematoma (P = 0.01), and the overall association after meta-analysis increased (P = 1.61 × 10-8; OR: 5.84; 95% CI: 3.16-10.76). ZBTB46 codes the zinc finger and BTB domain-containing protein 46 that acts as a transcription factor. In silico studies indicated that ZBTB46 is expressed in brain tissue by neurons and endothelial cells. Moreover, rs76484331 interacts with the promoter sites located at 20q13. In conclusion, we identified single nucleotide variants in the ZBTB46 gene associated with a higher risk of parenchymal haematoma following recombinant tissue-plasminogen activator treatment.

Accurate Prediction of Persistent Upper Extremity Impairment in Patients With Ischemic Stroke.

OBJECTIVE: To develop a simple and effective risk score for predicting which stroke patients will have persistent impairment of upper extremity motor function at 90 days. DESIGN: Post hoc analysis of clinical trial patients hospitalized with acute ischemic stroke who were followed for 90 days to determine functional outcome. SETTING: Patient were hospitalized at facilities across the United States. PARTICIPANTS: We created a harmonized cohort of individual patients (N=1653) from the NINDS tPA, ALIAS part 2, IMS-III, DEFUSE 3, and FAST-MAG trials. We split the cohort into balanced derivation and validation samples. INTERVENTIONS: Not applicable. MAIN OUTCOME MEASURES: The primary outcome was persistent arm impairment, defined as a National Institutes of Health Stroke Scale (NIHSS) arm domain score of 2 to 4 at 90 days in patients who had a 24-hour NIHSS arm score of 1 or more. We used least absolute shrinkage and selection operator regression to determine the elements of the persistent upper extremity impairment (PUPPI) index, which we validated as a predictive tool. RESULTS: We included 1653 patients (827 derivation, 826 validation), of whom 803 (48.6%) had persistent arm impairment. The PUPPI index gives 1 point each for age 55 years or older and NIHSS values of worse arm (4), worse leg (>2), facial palsy (3), and total NIHSS (≥10). The optimal cutpoint for the PUPPI index was 3 or greater, at which the area under the curve was greater than 0.75 for the derivation and validation cohorts and when using NIHSS values from either 24 hours or in a subacute or discharge time window. Results were similar across different levels of stroke severity. CONCLUSION: The PUPPI index uses readily available information to accurately predict persistent upper extremity motor impairment at 90 days poststroke. The PUPPI index can be administered in minutes and could be used as inclusion criterion in recovery-related clinical trials or, with additional development, as a prognostic tool for patients, caregivers, and clinicians.

Genomic characterization of autozygosity and recent inbreeding trends in all major breeds of US dairy cattle.

Maintaining a genetically diverse dairy cattle population is critical to preserving adaptability to future breeding goals and avoiding declines in fitness. This study characterized the genomic landscape of autozygosity and assessed trends in genetic diversity in 5 breeds of US dairy cattle. We analyzed a sizable genomic data set containing 4,173,679 pedigreed and genotyped animals of the Ayrshire, Brown Swiss, Guernsey, Holstein, and Jersey breeds. Runs of homozygosity (ROH) of 2 Mb or longer in length were identified in each animal. The within-breed means for number and the combined length of ROH were highest in Jerseys (62.66 ± 8.29 ROH and 426.24 ± 83.40 Mb, respectively; mean ± SD) and lowest in Ayrshires (37.24 ± 8.27 ROH and 265.05 ± 85.00 Mb, respectively). Short ROH were the most abundant, but moderate to large ROH made up the largest proportion of genome autozygosity in all breeds. In addition, we identified ROH islands in each breed. This revealed selection patterns for milk production, productive life, health, and reproduction in most breeds and evidence for parallel selective pressure for loci on chromosome 6 between Ayrshire and Brown Swiss and for loci on chromosome 20 between Holstein and Jersey. We calculated inbreeding coefficients using 3 different approaches, pedigree-based (FPED), marker-based using a genomic relationship matrix (FGRM), and segment-based using ROH (FROH). The average inbreeding coefficient ranged from 0.06 in Ayrshires and Brown Swiss to 0.08 in Jerseys and Holsteins using FPED, from 0.22 in Holsteins to 0.29 in Guernsey and Jerseys using FGRM, and from 0.11 in Ayrshires to 0.17 in Jerseys using FROH. In addition, the effective population size at past generations (5-100 generations ago), the yearly rate of inbreeding, and the effective population size in 3 recent periods (2000-2009, 2010-2014, and 2015-2018) were determined in each breed to ascertain current and historical trends of genetic diversity. We found a historical trend of decreasing effective population size in the last 100 generations in all breeds and breed differences in the effect of the recent implementation of genomic selection on inbreeding accumulation.

Truncation of IFT80 causes early embryonic loss in Holstein cattle associated with Holstein haplotype 2.

Recessive alleles represent genetic risk in populations that have undergone bottleneck events. We present a comprehensive framework for identification and validation of these genetic defects, including haplotype-based detection, variant selection from sequence data, and validation using knockout embryos. Holstein haplotype 2 (HH2), which causes embryonic death, was used to demonstrate the approach. Holstein haplotype 2 was identified using a deficiency-of-homozygotes approach and confirmed to negatively affect conception rate and stillbirths. Five carriers were present in a group of 183 sequenced Holstein bulls selected to maximize the coverage of unique haplotypes. Three variants concordant with haplotype calls were found in HH2: a high-priority frameshift mutation resulting, and 2 low-priority variants (1 synonymous variant, 1 premature stop codon). The frameshift in intraflagellar 80 (IFT80) was confirmed in a separate group of Holsteins from the 1000 Bull Genomes Project that shared no animals with the discovery set. IFT80-null embryos were generated by truncating the IFT80 transcript at exon 2 or 11 using a CRISPR-Cas9 system. Abattoir-derived oocytes were fertilized in vitro, and zygotes were injected at the one-cell stage either with a guide RNA and CAS9 mRNA complex (n = 100) or Cas9 mRNA (control, n = 100) before return to culture, and replicated 3 times. IFT80 is activated at the 8-cell stage, and IFT80-null embryos arrested at this stage of development, which is consistent with data from mouse hypomorphs and HH2 carrier-to-carrier matings. This frameshift in IFT80 on chromosome 1 at 107,172,615 bp (p.Leu381fs) disrupts WNT and hedgehog signaling, and is responsible for the death of homozygous embryos.

Single-cell transcriptomic analyses of dairy cattle ruminal epithelial cells during weaning.

Using the 10× Genomics Chromium Controller, we obtained scRNA-seq data of 5064 and 1372 individual cells from two Holstein calf ruminal epithelial tissues before and after weaning, respectively. We detected six distinct cell clusters, designated their cell types, and reported their marker genes. We then examined these clusters' underlining cell types and relationships by performing cell cycle, pseudotime trajectory, regulatory network, weighted gene co-expression network and gene ontology analyses. By integrating these cell marker genes with Holstein GWAS signals, we found they were enriched for animal production and body conformation traits. Finally, we confirmed their cell identities by comparing them with human and mouse stomach epithelial cells. This study presents an initial effort to implement single-cell transcriptomic analysis in cattle, and demonstrates ruminal tissue epithelial cell types and their developments during weaning, opening the door for new discoveries about tissue/cell type roles in complex traits at single-cell resolution.