Exome sequencing identifies novel genes underlying primary congenital glaucoma in the National Birth Defects Prevention Study.

BACKGROUND: Primary congenital glaucoma (PCG) affects approximately 1 in 10,000 live born infants in the United States (U.S.). PCG has a autosomal recessive inheritance pattern, and variable expressivity and reduced penetrance have been reported. Likely causal variants in the most commonly mutated gene, CYP1B1, are less prevalent in the U.S., suggesting that alternative genes may contribute to the condition. This study utilized exome sequencing to investigate the genetic architecture of PCG in the U.S. and to identify novel genes and variants. METHODS: We studied 37 family trios where infants had PCG and were part of the National Birth Defects Prevention Study (births 1997-2011), a U.S. multicenter study of birth defects. Samples underwent exome sequencing and sequence reads were aligned to the human reference sample (NCBI build 37/hg19). Variant filtration was conducted under de novo and Mendelian inheritance models using GEMINI. RESULTS: Among candidate variants, CYP1B1 was most represented (five trios, 13.5%). Twelve probands (32%) had potentially pathogenic variants in other genes not previously linked to PCG but important in eye development and/or to underlie Mendelian conditions with potential phenotypic overlap (e.g., CRYBB2, RXRA, GLI2). CONCLUSION: Variation in the genes identified in this population-based study may help to further explain the genetics of PCG.

ClOneHORT: Approaches for Improved Fidelity in Generative Models of Synthetic Genomes.

Motivation: Deep generative models have the potential to overcome difficulties in sharing individual-level genomic data by producing synthetic genomes that preserve the genomic associations specific to a cohort while not violating the privacy of any individual cohort member. However, there is significant room for improvement in the fidelity and usability of existing synthetic genome approaches. Results: We demonstrate that when combined with plentiful data and with population-specific selection criteria, deep generative models can produce synthetic genomes and cohorts that closely model the original populations. Our methods improve fidelity in the site-frequency spectra and linkage disequilibrium decay and yield synthetic genomes that can be substituted in downstream local ancestry inference analysis, recreating results with .91 to .94 accuracy. Availability: The model described in this paper is freely available at github.com/rlaboulaye/clonehort .

Targeted grazing reduces a widespread wetland plant invader with minimal nutrient impacts, yet native community recovery is limited.

Targeted grazing to control undesirable plant species is increasingly of interest across a diversity of ecosystems, particularly as an alternative or complement to widely used herbicides. However, there are limited comprehensive evaluations of targeted grazing that evaluate both invasive species management effectiveness and potential negative effects on the ecosystem. Phragmites australis, a tall-statured, dense perennial invasive grass from Eurasia, is a pervasive problem in wetlands across the North American continent. As with many invasive species where management has historically relied on herbicides and resistance is a growing concern, land managers seek viable alternatives that have minimal negative ecosystem impacts. Grazing has been used for millennia to manage native Phragmites in Europe. Similarly, in its invasive range within North America, small-scale studies suggest Phragmites may be suppressed by grazers. Yet, the effectiveness of grazing at large scales and its effects on broader ecosystem properties remain largely unknown. We evaluated the influence of targeted grazing on vegetation, soil nutrients, and water nutrients over two years in large plots (∼300x the size of previous studies). We also tested the effects of mowing, a treatment that can be used to facilitate grazer access to large, dense Phragmites stands. In line with our predictions, we found that cattle grazing effectively suppressed invasive Phragmites over two years. Mowing reduced litter, and moderately reduced standing dead Phragmites, both of which suppress native plant germination in this system. However, these reductions in Phragmites were not accompanied by indications of native plant community recovery, as we had optimistically predicted. Despite the potential for grazing to reduce nutrient sequestration by plants and fertilize soils, we were surprised to find no clear negative effects of grazing on nutrient mobilization to groundwater or floodwater. Taken together, our findings indicate that targeted grazing, when implemented at broad scales over short time frames, is effective at achieving invasive plant management goals without sizable nutrient impacts. However, additional steps will be needed to achieve the restoration of diverse, robust native plant communities.

Genome-wide association study reveals shared and distinct genetic architecture underlying fatty acid and bioactive oxylipin metabolites in the Hispanic Community Health Study/Study of Latinos (HCHS/SOL).

Bioactive fatty acid-derived oxylipin molecules play key roles in mediating inflammation and oxidative stress, which underlie many chronic diseases. Circulating levels of fatty acids and oxylipins are influenced by both environmental and genetic factors; characterizing the genetic architecture of bioactive lipids could yield new insights into underlying biological pathways. Thus, we performed a genome wide association study (GWAS) of n=81 fatty acids and oxylipins in n=11,584 Hispanic Community Health Study/Study of Latinos (HCHS/SOL) participants with genetic and lipidomic data measured at study baseline (58.6% female, mean age = 46.1 years, standard deviation = 13.8 years). Additionally, given the effects of central obesity on inflammation, we examined interactions with waist circumference using two-degree-of-freedom joint tests. Heritability estimates ranged from 0% to 47.9%, and 48 of the 81oxylipins and fatty acids were significantly heritable. Moreover, 40 (49.4%) of the 81 oxylipins and fatty acids had at least one genome-wide significant (p< 6.94E-11) variant resulting in 19 independent genetic loci involved in fatty acid and oxylipin synthesis, as well as downstream pathways. Four loci (lead variant minor allele frequency [MAF] range: 0.08-0.50), including the desaturase-encoding FADS and the OATP1B1 transporter protein-encoding SLCO1B1, exhibited associations with four or more fatty acids and oxylipins. The majority of the 15 remaining loci (87.5%) (lead variant MAF range = 0.03-0.45, mean = 0.23) were only associated with one oxylipin or fatty acid, demonstrating evidence of distinct genetic effects. Finally, while most loci identified in two-degree-of-freedom tests were previously identified in our main effects analyses, we also identified an additional rare variant (MAF = 0.002) near CARS2, a locus previously implicated in inflammation. Our analyses revealed shared and distinct genetic architecture underlying fatty acids and oxylipins, providing insights into genetic factors and motivating future multi-omics work to characterize these compounds and elucidate their roles in disease pathways.

Integrating Genetic and Transcriptomic Data to Identify Genes Underlying Obesity Risk Loci.

Genome-wide association studies (GWAS) have identified numerous body mass index (BMI) loci. However, most underlying mechanisms from risk locus to BMI remain unknown. Leveraging omics data through integrative analyses could provide more comprehensive views of biological pathways on BMI. We analyzed genotype and blood gene expression data in up to 5,619 samples from the Framingham Heart Study (FHS). Using 3,992 single nucleotide polymorphisms (SNPs) at 97 BMI loci and 20,692 transcripts within 1 Mb, we performed separate association analyses of transcript with BMI and SNP with transcript (PBMI and PSNP, respectively) and then a correlated meta-analysis between the full summary data sets (PMETA). We identified transcripts that met Bonferroni-corrected significance for each omic, were more significant in the correlated meta-analysis than each omic, and were at least nominally associated with BMI in FHS data. Among 308 significant SNP-transcript-BMI associations, we identified seven genes (NT5C2, GSTM3, SNAPC3, SPNS1, TMEM245, YPEL3, and ZNF646) in five association regions. Using an independent sample of blood gene expression data, we validated results for SNAPC3 and YPEL3. We tested for generalization of these associations in hypothalamus, nucleus accumbens, and liver and observed significant (PMETA<0.05 & PMETA

Clonal hematopoiesis driven by mutated DNMT3A promotes inflammatory bone loss.

Clonal hematopoiesis of indeterminate potential (CHIP) arises from aging-associated acquired mutations in hematopoietic progenitors, which display clonal expansion and produce phenotypically altered leukocytes. We associated CHIP-DNMT3A mutations with a higher prevalence of periodontitis and gingival inflammation among 4,946 community-dwelling adults. To model DNMT3A-driven CHIP, we used mice with the heterozygous loss-of-function mutation R878H, equivalent to the human hotspot mutation R882H. Partial transplantation with Dnmt3aR878H/+ bone marrow (BM) cells resulted in clonal expansion of mutant cells into both myeloid and lymphoid lineages and an elevated abundance of osteoclast precursors in the BM and osteoclastogenic macrophages in the periphery. DNMT3A-driven clonal hematopoiesis in recipient mice promoted naturally occurring periodontitis and aggravated experimentally induced periodontitis and arthritis, associated with enhanced osteoclastogenesis, IL-17-dependent inflammation and neutrophil responses, and impaired regulatory T cell immunosuppressive activity. DNMT3A-driven clonal hematopoiesis and, subsequently, periodontitis were suppressed by rapamycin treatment. DNMT3A-driven CHIP represents a treatable state of maladaptive hematopoiesis promoting inflammatory bone loss.

Whole genome sequencing based analysis of inflammation biomarkers in the Trans-Omics for Precision Medicine (TOPMed) consortium.

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38 465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program (with varying sample size by trait, where the minimum sample size was n = 737 for MMP-1). We identified 22 distinct single-variant associations across 6 traits-E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin-that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

Determinants of mosaic chromosomal alteration fitness.

Clonal hematopoiesis (CH) is characterized by the acquisition of a somatic mutation in a hematopoietic stem cell that results in a clonal expansion. These driver mutations can be single nucleotide variants in cancer driver genes or larger structural rearrangements called mosaic chromosomal alterations (mCAs). The factors that influence the variations in mCA fitness and ultimately result in different clonal expansion rates are not well understood. We used the Passenger-Approximated Clonal Expansion Rate (PACER) method to estimate clonal expansion rate as PACER scores for 6,381 individuals in the NHLBI TOPMed cohort with gain, loss, and copy-neutral loss of heterozygosity mCAs. Our mCA fitness estimates, derived by aggregating per-individual PACER scores, were correlated (R2 = 0.49) with an alternative approach that estimated fitness of mCAs in the UK Biobank using population-level distributions of clonal fraction. Among individuals with JAK2 V617F clonal hematopoiesis of indeterminate potential or mCAs affecting the JAK2 gene on chromosome 9, PACER score was strongly correlated with erythrocyte count. In a cross-sectional analysis, genome-wide association study of estimates of mCA expansion rate identified a TCL1A locus variant associated with mCA clonal expansion rate, with suggestive variants in NRIP1 and TERT.

Prairie restoration promotes the abundance and diversity of mutualistic arbuscular mycorrhizal fungi.

Predicting how biological communities assemble in restored ecosystems can assist in conservation efforts, but most research has focused on plants, with relatively little attention paid to soil microbial organisms that plants interact with. Arbuscular mycorrhizal (AM) fungi are an ecologically significant functional group of soil microbes that form mutualistic symbioses with plants and could therefore respond positively to plant community restoration. To evaluate the effects of plant community restoration on AM fungi, we compared AM fungal abundance, species richness, and community composition of five annually cultivated, conventionally managed agricultural fields with paired adjacent retired agricultural fields that had undergone prairie restoration 5-9 years prior to sampling. We hypothesized that restoration stimulates AM fungal abundance and species richness, particularly for disturbance-sensitive taxa, and that gains of new taxa would not displace AM fungal species present prior to restoration due to legacy effects. AM fungal abundance was quantified by measuring soil spore density and root colonization. AM fungal species richness and community composition were determined in soils and plant roots using DNA high-throughput sequencing. Soil spore density was 2.3 times higher in restored prairies compared to agricultural fields, but AM fungal root colonization did not differ between land use types. AM fungal species richness was 2.7 and 1.4 times higher in restored prairies versus agricultural fields for soil and roots, respectively. The abundance of Glomeraceae, a disturbance-tolerant family, decreased by 25% from agricultural to restored prairie soils but did not differ in plant roots. The abundance of Claroideoglomeraceae and Diversisporaceae, both disturbance-sensitive families, was 4.6 and 3.2 times higher in restored prairie versus agricultural soils, respectively. Species turnover was higher than expected relative to a null model, indicating that AM fungal species were gained by replacement. Our findings demonstrate that restoration can promote a relatively rapid increase in the abundance and diversity of soil microbial communities that had been degraded by decades of intensive land use, and community compositional change can be predicted by the disturbance tolerance of soil microbial taxonomic and functional groups.

Sex-specific single cell-level transcriptomic signatures of Rett syndrome disease progression.

Dominant X-linked diseases are uncommon due to female X chromosome inactivation (XCI). While random XCI usually protects females against X-linked mutations, Rett syndrome (RTT) is a female neurodevelopmental disorder caused by heterozygous MECP2 mutation. After 6-18 months of typical neurodevelopment, RTT girls undergo poorly understood regression. We performed longitudinal snRNA-seq on cerebral cortex in a construct-relevant Mecp2e1 mutant mouse model of RTT, revealing transcriptional effects of cell type, mosaicism, and sex on progressive disease phenotypes. Across cell types, we observed sex differences in the number of differentially expressed genes (DEGs) with 6x more DEGs in mutant females than males. Unlike males, female DEGs emerged prior to symptoms, were enriched for homeostatic gene pathways in distinct cell types over time, and correlated with disease phenotypes and human RTT cortical cell transcriptomes. Non-cell-autonomous effects were prominent and dynamic across disease progression of Mecp2e1 mutant females, indicating wild-type-expressing cells normalizing transcriptional homeostasis. These results improve understanding of RTT progression and treatment.

Methylation patterns associated with C-reactive protein in racially and ethnically diverse populations.

Systemic low-grade inflammation is a feature of chronic disease. C-reactive protein (CRP) is a common biomarker of inflammation and used as an indicator of disease risk; however, the role of inflammation in disease is not completely understood. Methylation is an epigenetic modification in the DNA which plays a pivotal role in gene expression. In this study we evaluated differential DNA methylation patterns associated with blood CRP level to elucidate biological pathways and genetic regulatory mechanisms to improve the understanding of chronic inflammation. The racially and ethnically diverse participants in this study were included as 50% White, 41% Black or African American, 7% Hispanic or Latino/a, and 2% Native Hawaiian, Asian American, American Indian, or Alaska Native (total n = 13,433) individuals. We replicated 113 CpG sites from 87 unique loci, of which five were novel (CADM3, NALCN, NLRC5, ZNF792, and cg03282312), across a discovery set of 1,150 CpG sites associated with CRP level (p < 1.2E-7). The downstream pathways affected by DNA methylation included the identification of IFI16 and IRF7 CpG-gene transcript pairs which contributed to the innate immune response gene enrichment pathway along with NLRC5, NOD2, and AIM2. Gene enrichment analysis also identified the nuclear factor-kappaB transcription pathway. Using two-sample Mendelian randomization (MR) we inferred methylation at three CpG sites as causal for CRP levels using both White and Black or African American MR instrument variables. Overall, we identified novel CpG sites and gene transcripts that could be valuable in understanding the specific cellular processes and pathogenic mechanisms involved in inflammation.

Transplanted sagebrush "wildlings" exhibit higher survival than greenhouse-grown tubelings yet both recruit new plants.

BACKGROUND: Land uses such as crop production, livestock grazing, mining, and urban development have contributed to degradation of drylands worldwide. Loss of big sagebrush (Artemisia tridentata) on disturbed drylands across the western U.S. has prompted massive efforts to re-establish this foundational species. There has been growing interest in avoiding the severe limitations experienced by plants at the seed and seedling stages by instead establishing plants from containerized greenhouse seedlings ("tubelings"). In some settings, a potential alternative approach is to transplant larger locally-collected plants ("wildlings"). We compared the establishment of mountain big sagebrush (A. tridentata ssp. vaseyana) from tubelings vs. wildlings in southeastern Idaho. A mix of native and non-native grass and forb species was drill-seeded in a pasture previously dominated by the introduced forage grass, smooth brome (Bromus inermis). We then established 80 m x 80 m treatment plots and planted sagebrush tubelings (n = 12 plots, 1200 plants) and wildlings (n = 12 plots, 1200 plants). We also established seeded plots (n = 12) and untreated control plots (n = 6) for long-term comparison. We tracked project expenses in order to calculate costs of using tubelings vs. wildlings as modified by probability of success. RESULTS: There was high (79%) tubeling and low (10%) wildling mortality within the first year. Three years post-planting, chance of survival for wildlings was significantly higher than that of tubelings (85% and 14% respectively). Despite high up-front costs of planting wildlings, high survival rates resulted in their being < 50% of the cost of tubelings on a per-surviving plant basis. Additionally, by the third year post-planting 34% of surviving tubelings and 95% of surviving wildlings showed evidence of reproduction (presence / absence of flowering stems), and the two types of plantings recruited new plants via seed (3.7 and 2.4 plants, respectively, per surviving tubeling/wildling). CONCLUSIONS: Our results indicate that larger plants with more developed root systems (wildlings) may be a promising avenue for increasing early establishment rates of sagebrush plants in restoration settings. Our results also illustrate the potential for tubelings and wildlings to improve restoration outcomes by "nucleating" the landscape via recruitment of new plants during ideal climate conditions.

Type 1 Diabetes Genetic Risk in 109,954 Veterans With Adult-Onset Diabetes: The Million Veteran Program (MVP).

OBJECTIVE: To characterize high type 1 diabetes (T1D) genetic risk in a population where type 2 diabetes (T2D) predominates. RESEARCH DESIGN AND METHODS: Characteristics typically associated with T1D were assessed in 109,594 Million Veteran Program participants with adult-onset diabetes, 2011-2021, who had T1D genetic risk scores (GRS) defined as low (0 to <45%), medium (45 to <90%), high (90 to <95%), or highest (≥95%). RESULTS: T1D characteristics increased progressively with higher genetic risk (P < 0.001 for trend). A GRS ≥90% was more common with diabetes diagnoses before age 40 years, but 95% of those participants were diagnosed at age ≥40 years, and their characteristics resembled those of individuals with T2D in mean age (64.3 years) and BMI (32.3 kg/m2). Compared with the low-risk group, the highest-risk group was more likely to have diabetic ketoacidosis (low GRS 0.9% vs. highest GRS 3.7%), hypoglycemia prompting emergency visits (3.7% vs. 5.8%), outpatient plasma glucose <50 mg/dL (7.5% vs. 13.4%), a shorter median time to start insulin (3.5 vs. 1.4 years), use of a T1D diagnostic code (16.3% vs. 28.1%), low C-peptide levels if tested (1.8% vs. 32.4%), and glutamic acid decarboxylase antibodies (6.9% vs. 45.2%), all P < 0.001. CONCLUSIONS: Characteristics associated with T1D were increased with higher genetic risk, and especially with the top 10% of risk. However, the age and BMI of those participants resemble those of people with T2D, and a substantial proportion did not have diagnostic testing or use of T1D diagnostic codes. T1D genetic screening could be used to aid identification of adult-onset T1D in settings in which T2D predominates.

Experimental Models to Study Immune Dysfunction in the Pathogenesis of Parkinson's Disease.

Parkinson's disease (PD) is a chronic, age-related, progressive multisystem disease associated with neuroinflammation and immune dysfunction. This review discusses the methodological approaches used to study the changes in central and peripheral immunity in PD, the advantages and limitations of the techniques, and their applicability to humans. Although a single animal model cannot replicate all pathological features of the human disease, neuroinflammation is present in most animal models of PD and plays a critical role in understanding the involvement of the immune system (IS) in the pathogenesis of PD. The IS and its interactions with different cell types in the central nervous system (CNS) play an important role in the pathogenesis of PD. Even though culture models do not fully reflect the complexity of disease progression, they are limited in their ability to mimic long-term effects and need validation through in vivo studies. They are an indispensable tool for understanding the interplay between the IS and the pathogenesis of this disease. Understanding the immune-mediated mechanisms may lead to potential therapeutic targets for the treatment of PD. We believe that the development of methodological guidelines for experiments with animal models and PD patients is crucial to ensure the validity and consistency of the results.

Genetic and transcriptomic landscape of colonic diverticulosis.

OBJECTIVE: Colonic diverticulosis is a prevalent condition among older adults, marked by the presence of thin-walled pockets in the colon wall that can become inflamed, infected, haemorrhage or rupture. We present a case-control genetic and transcriptomic study aimed at identifying the genetic and cellular determinants underlying this condition and the relationship with other gastrointestinal disorders. DESIGN: We conducted DNA and RNA sequencing on colonic tissue from 404 patients with (N=172) and without (N=232) diverticulosis. We investigated variation in the transcriptome associated with diverticulosis and further integrated this variation with single-cell RNA-seq data from the human intestine. We also integrated our expression quantitative trait loci with genome-wide association study using Mendelian randomisation (MR). Furthermore, a Polygenic Risk Score analysis gauged associations between diverticulosis severity and other gastrointestinal disorders. RESULTS: We discerned 38 genes with differential expression and 17 with varied transcript usage linked to diverticulosis, indicating tissue remodelling as a primary diverticula formation mechanism. Diverticula formation was primarily linked to stromal and epithelial cells in the colon including endothelial cells, myofibroblasts, fibroblasts, goblet, tuft, enterocytes, neurons and glia. MR highlighted five genes including CCN3, CRISPLD2, ENTPD7, PHGR1 and TNFSF13, with potential causal effects on diverticulosis. Notably, ENTPD7 upregulation was confirmed in diverticulosis cases. Additionally, diverticulosis severity was positively correlated with genetic predisposition to diverticulitis. CONCLUSION: Our results suggest that tissue remodelling is a primary mechanism for diverticula formation. Individuals with an increased genetic proclivity to diverticulitis exhibit a larger numbers of diverticula on colonoscopy.

Serum Metabolomic Markers of Protein-Rich Foods and Incident CKD: Results From the Atherosclerosis Risk in Communities Study.

Rationale & Objective: While urine excretion of nitrogen estimates the total protein intake, biomarkers of specific dietary protein sources have been sparsely studied. Using untargeted metabolomics, this study aimed to identify serum metabolomic markers of 6 protein-rich foods and to examine whether dietary protein-related metabolites are associated with incident chronic kidney disease (CKD). Study Design: Prospective cohort study. Setting & Participants: A total of 3,726 participants from the Atherosclerosis Risk in Communities study without CKD at baseline. Exposures: Dietary intake of 6 protein-rich foods (fish, nuts, legumes, red and processed meat, eggs, and poultry), serum metabolites. Outcomes: Incident CKD (estimated glomerular filtration rate < 60 mL/min/1.73 m2 with ≥25% estimated glomerular filtration rate decline relative to visit 1, hospitalization or death related to CKD, or end-stage kidney disease). Analytical Approach: Multivariable linear regression models estimated cross-sectional associations between protein-rich foods and serum metabolites. C statistics assessed the ability of the metabolites to improve the discrimination of highest versus lower 3 quartiles of intake of protein-rich foods beyond covariates (demographics, clinical factors, health behaviors, and the intake of nonprotein food groups). Cox regression models identified prospective associations between protein-related metabolites and incident CKD. Results: Thirty significant associations were identified between protein-rich foods and serum metabolites (fish, n = 8; nuts, n = 5; legumes, n = 0; red and processed meat, n = 5; eggs, n = 3; and poultry, n = 9). Metabolites collectively and significantly improved the discrimination of high intake of protein-rich foods compared with covariates alone (difference in C statistics = 0.033, 0.051, 0.003, 0.024, and 0.025 for fish, nuts, red and processed meat, eggs, and poultry-related metabolites, respectively; P < 1.00 × 10-16 for all). Dietary intake of fish was positively associated with 1-docosahexaenoylglycerophosphocholine (22:6n3), which was inversely associated with incident CKD (HR, 0.82; 95% CI, 0.75-0.89; P = 7.81 × 10-6). Limitations: Residual confounding and sample-storage duration. Conclusions: We identified candidate biomarkers of fish, nuts, red and processed meat, eggs, and poultry. A fish-related metabolite, 1-docosahexaenoylglycerophosphocholine (22:6n3), was associated with a lower risk of CKD.

In this study, we aimed to identify associations between protein-rich foods (fish, nuts, legumes, red and processed meat, eggs, and poultry) and serum metabolites, which are small biological molecules involved in metabolism. Metabolites significantly associated with a protein-rich food individually and collectively improved the discrimination of the respective protein-rich food, suggesting that these metabolites should be prioritized in future diet biomarker research. We also studied associations between significant diet-related metabolites and incident kidney disease. One fish-related metabolite was associated with a lower kidney disease risk. This finding supports the recent nutritional guidelines recommending a Mediterranean diet, which includes fish as the main dietary protein source.

Metabolomic Association and Risk Prediction With Heart Failure in Older Adults.

BACKGROUND: Older adults have markedly increased risks of heart failure (HF), specifically HF with preserved ejection fraction (HFpEF). Identifying novel biomarkers can help in understanding HF pathogenesis and improve at-risk population identification. This study aimed to identify metabolites associated with incident HF, HFpEF, and HF with reduced ejection fraction and examine risk prediction in older adults. METHODS: Untargeted metabolomic profiling was performed in Black and White adults from the ARIC study (Atherosclerosis Risk in Communities) visit 5 (n=3719; mean age, 75 years). We applied Cox regressions to identify metabolites associated with incident HF and its subtypes. The metabolite risk score (MRS) was constructed and examined for associations with HF, echocardiographic measures, and HF risk prediction. Independent samples from visit 3 (n=1929; mean age, 58 years) were used for replication. RESULTS: Sixty metabolites (hazard ratios range, 0.79-1.49; false discovery rate, <0.05) were associated with incident HF after adjusting for clinical risk factors, eGFR, and NT-proBNP (N-terminal pro-B-type natriuretic peptide). Mannonate, a hydroxy acid, was replicated (hazard ratio, 1.36 [95% CI, 1.19-1.56]) with full adjustments. MRS was associated with an 80% increased risk of HF per SD increment, and the highest MRS quartile had 8.7× the risk of developing HFpEF than the lowest quartile. High MRS was also associated with unfavorable values of cardiac structure and function. Adding MRS over clinical risk factors and NT-proBNP improved 5-year HF risk prediction C statistics from 0.817 to 0.850 (∆C, 0.033 [95% CI, 0.017-0.047]). The association between MRS and incident HF was replicated after accounting for clinical risk factors (P<0.05). CONCLUSIONS: Novel metabolites associated with HF risk were identified, elucidating disease pathways, specifically HFpEF. An MRS was associated with HF risk and improved 5-year risk prediction in older adults, which may assist at at-risk population identification.

Exploiting meta-analysis of genome-wide interaction with serum 25-hydroxyvitamin D to identify novel genetic loci associated with pulmonary function.

BACKGROUND: Higher 25-hydroxyvitamin D (25(OH)D) concentrations in serum has a positive association with pulmonary function. Investigating genome-wide interactions with 25(OH)D may reveal new biological insights into pulmonary function. OBJECTIVES: We aimed to identify novel genetic variants associated with pulmonary function by accounting for 25(OH)D interactions. METHODS: We included 211,264 participants from the observational United Kingdom Biobank study with pulmonary function tests (PFTs), genome-wide genotypes, and 25(OH)D concentrations from 4 ancestral backgrounds-European, African, East Asian, and South Asian. Among PFTs, we focused on forced expiratory volume in the first second (FEV1) and forced vital capacity (FVC) because both were previously associated with 25(OH)D. We performed genome-wide association study (GWAS) analyses that accounted for variant×25(OH)D interaction using the joint 2 degree-of-freedom (2df) method, stratified by participants' smoking history and ancestry, and meta-analyzed results. We evaluated interaction effects to determine how variant-PFT associations were modified by 25(OH)D concentrations and conducted pathway enrichment analysis to examine the biological relevance of our findings. RESULTS: Our GWAS meta-analyses, accounting for interaction with 25(OH)D, revealed 30 genetic variants significantly associated with FEV1 or FVC (P2df <5.00×10-8) that were not previously reported for PFT-related traits. These novel variant signals were enriched in lung function-relevant pathways, including the p38 MAPK pathway. Among variants with genome-wide-significant 2df results, smoking-stratified meta-analyses identified 5 variants with 25(OH)D interactions that influenced FEV1 in both smoking groups (never smokers P1df interaction<2.65×10-4; ever smokers P1df interaction<1.71×10-5); rs3130553, rs2894186, rs79277477, and rs3130929 associations were only evident in never smokers, and the rs4678408 association was only found in ever smokers. CONCLUSION: Genetic variant associations with lung function can be modified by 25(OH)D, and smoking history can further modify variant×25(OH)D interactions. These results expand the known genetic architecture of pulmonary function and add evidence that gene-environment interactions, including with 25(OH)D and smoking, influence lung function.

Differential structure and function of phosphorus-mineralizing microbial communities in organic and upper mineral soil horizons across a temperate rainforest chronosequence.

Microbial community structure and function were assessed in the organic and upper mineral soil across a ~4000-year dune-based chronosequence at Big Bay, New Zealand, where total P declined and the proportional contribution of organic soil in the profile increased with time. We hypothesized that the organic and mineral soils would show divergent community evolution over time with a greater dependency on the functionality of phosphatase genes in the organic soil layer as it developed. The structure of bacterial, fungal, and phosphatase-harbouring communities was examined in both horizons across 3 dunes using amplicon sequencing, network analysis, and qPCR. The soils showed a decline in pH and total phosphorus (P) over time with an increase in phosphatase activity. The organic horizon had a wider diversity of Class A (phoN/phoC) and phoD-harbouring communities and a more complex microbiome, with hub taxa that correlated with P. Bacterial diversity declined in both horizons over time, with enrichment of Planctomycetes and Acidobacteria. More complex fungal communities were evident in the youngest dune, transitioning to a dominance of Ascomycota in both soil horizons. Higher phosphatase activity in older dunes was driven by less diverse P-mineralizing communities, especially in the organic horizon.

Acute Adenoviral Infection Elicits an Arrhythmogenic Substrate Prior to Myocarditis.

BACKGROUND: Viral cardiac infection represents a significant clinical challenge encompassing several etiological agents, disease stages, complex presentation, and a resulting lack of mechanistic understanding. Myocarditis is a major cause of sudden cardiac death in young adults, where current knowledge in the field is dominated by later disease phases and pathological immune responses. However, little is known regarding how infection can acutely induce an arrhythmogenic substrate before significant immune responses. Adenovirus is a leading cause of myocarditis, but due to species specificity, models of infection are lacking, and it is not understood how adenoviral infection may underlie sudden cardiac arrest. Mouse adenovirus type-3 was previously reported as cardiotropic, yet it has not been utilized to understand the mechanisms of cardiac infection and pathology. METHODS: We have developed mouse adenovirus type-3 infection as a model to investigate acute cardiac infection and molecular alterations to the infected heart before an appreciable immune response or gross cardiomyopathy. RESULTS: Optical mapping of infected hearts exposes decreases in conduction velocity concomitant with increased Cx43Ser368 phosphorylation, a residue known to regulate gap junction function. Hearts from animals harboring a phospho-null mutation at Cx43Ser368 are protected against mouse adenovirus type-3-induced conduction velocity slowing. Additional to gap junction alterations, patch clamping of mouse adenovirus type-3-infected adult mouse ventricular cardiomyocytes reveals prolonged action potential duration as a result of decreased IK1 and IKs current density. Turning to human systems, we find human adenovirus type-5 increases phosphorylation of Cx43Ser368 and disrupts synchrony in human induced pluripotent stem cell-derived cardiomyocytes, indicating common mechanisms with our mouse whole heart and adult cardiomyocyte data. CONCLUSIONS: Together, these findings demonstrate that adenoviral infection creates an arrhythmogenic substrate through direct targeting of gap junction and ion channel function in the heart. Such alterations are known to precipitate arrhythmias and likely contribute to sudden cardiac death in acutely infected patients.