Gene-based Analysis Research Articles

Revealing the Diversity of the Mycobiome in Different Phases of Ticks: ITS Gene-Based Analysis

Ticks are obligate ectoparasites and vectors of a variety of pathogens in humans and animals. Certain tick-borne pathogens (TBPs) have been identified as the cause of zoonoses, posing potentially significant threats to the human health and livestock industries. Fungi are one of the major TBPs that can affect ticks and cause disease in humans. At present, there are few studies on the diversity of fungal microbial communities carried by Ixodes. Therefore, profiling tick-borne fungi will contribute to understanding the tick-fungal interaction. This study evaluated the community profile and differences in the fungal microbiome in Ixodidae collected on parasitic ticks or nonparasitic ticks in Wuwei, Gansu Province, China. The Shannon index, Simpson index, and Richness index were used to evaluate the diversity of mycobiome. Principle coordinates analysis (PCoA) was conducted to determine patterns of diversity in mycobiome. Using correlation analysis to determine the correlation of mycobiome. The results show that the high-throughput sequencing of the internal transcribed spacer gene generated 3,634,943 raw reads and 7,482 amplicon sequence variants. The dominant tick species in this region was Dermacentor nuttalli (Ixodidae). The mycobiome belonged to four classes—Dothideomycetes, Sordariomycetes, Ustilaginomycetes, and Tremellomycetes—and more than 261 genera, the most abundant genera were Cladosporium, Purpureocillium, Aureobasidium, Tranzscheliella, and Sporormiella. Alpha diversity indicated that the abundance and evenness of mycobiome were marginally higher in nonparasitic ticks than in parasitic ticks. PCoA showed that the community structures of parasitic ticks vary from nonparasitic ticks, samples from nonparasitic ticks tended to cluster more closely than those from the parasitic ticks. Correlation analysis indicated that there was a significant positive correlation or negative correlation between the mycobiome. Our results indicate that the mycobiome carried by Dermacentor nuttalli had rich diversity, and there was a significant difference in mycobiome between parasitic ticks and nonparasitic ticks. These findings may conducive to understand the complex interaction between ticks and commensal fungi and provide help for the further exploration of the behavioral characteristics of ticks and formulation of effective biological control measures.

Genetic interactions of schizophrenia using gene-based statistical epistasis exclusively identify nervous system-related pathways and key hub genes.

Background: The relationship between genotype and phenotype is governed by numerous genetic interactions (GIs), and the mapping of GI networks is of interest for two main reasons: 1) By modelling biological robustness, GIs provide a powerful opportunity to infer compensatory biological mechanisms via the identification of functional relationships between genes, which is of interest for biological discovery and translational research. Biological systems have evolved to compensate for genetic (i.e., variations and mutations) and environmental (i.e., drug efficacy) perturbations by exploiting compensatory relationships between genes, pathways and biological processes; 2) GI facilitates the identification of the direction (alleviating or aggravating interactions) and magnitude of epistatic interactions that influence the phenotypic outcome. The generation of GIs for human diseases is impossible using experimental biology approaches such as systematic deletion analysis. Moreover, the generation of disease-specific GIs has never been undertaken in humans. Methods: We used our Indian schizophrenia case-control (case-816, controls-900) genetic dataset to implement the workflow. Standard GWAS sample quality control procedure was followed. We used the imputed genetic data to increase the SNP coverage to analyse epistatic interactions across the genome comprehensively. Using the odds ratio (OR), we identified the GIs that increase or decrease the risk of a disease phenotype. The SNP-based epistatic results were transformed into gene-based epistatic results. Results: We have developed a novel approach by conducting gene-based statistical epistatic analysis using an Indian schizophrenia case-control genetic dataset and transforming these results to infer GIs that increase the risk of schizophrenia. There were ∼9.5 million GIs with a p-value 1 10-5. Approximately 4.8 million GIs showed an increased risk (OR > 1.0), while ∼4.75 million GIs had a decreased risk (OR <1.0) for schizophrenia. Conclusion: Unlike model organisms, this approach is specifically viable in humans due to the availability of abundant disease-specific genome-wide genotype datasets. The study exclusively identified brain/nervous system-related processes, affirming the findings. This computational approach fills a critical gap by generating practically non-existent heritable disease-specific human GIs from human genetic data. These novel datasets can train innovative deep-learning models, potentially surpassing the limitations of conventional GWAS.

Poultry litter biochar soil amendment affects microbial community structures, promotes phosphorus cycling and growth of barley (Hordeumvulgare)

Phosphorus (P) is a non-replaceable, finite component of fertilizers. The imbalanced resource distribution and possible depletion of P impose challenges on current crop production worldwide. The aim of this study was to assess the impact of poultry litter biochar on plant growth and P mobilizing capability of the microbiome in comparison to a mineral fertilizer application. Spring barley (Hordeum vulgare) was grown in microcosms using a P-limited soil, fertilized with 0 (control), 50 (fertilizer) kg P ha−1 or a poultry litter biochar amendment (biochar, 2% w/w). Biochar amended rhizospheres had significantly higher phosphonate-utilizing bacteria, phoC and phoD gene (acid and alkaline phosphatase) copy numbers and improved P availability. Spring barley dry matter yields were significantly higher for biochar and fertilizer over control; however, P uptake with biochar was higher than with fertilizer. Furthermore, biochar had higher arbuscular mycorrhizal colonization and significantly raised soil pH. Fingerprint-analysis showed significant differences between all treatments for bacterial and fungal communities. 16S rRNA gene-based sequencing analysis revealed increased relative abundance of the phyla Actinobacteriota and Chloroflexi in biochar compared to control, potentially contributing to the ameliorated plant growth conditions. Pearson correlations of both phyla was positive with a range of P cycling variables as well as Morgan's P but negative with acid phosphatase activity. FAPROTAX analysis revealed positive correlations of aromatic compound degradation with phoC and phoD gene abundance, highlighting a possible link between biodegradation and P release. In conclusion, poultry litter biochar could potentially replace mineral P fertilizer for sustainable plant growth in P depleted soil environments.

Genome-wide characterization of soybean malate dehydrogenase genes reveals a positive role of GmMDH2 in salt stress response

Malate dehydrogenase (MDH) is a widely expressed enzyme that plays a key role in plant growth, development, and the stress response. However, information on MDH genes in the soybean genome is limited. Seventeen members of the soybean MDH family were identified by genome-wide analysis, and the presence of conserved protein motifs was analyzed. The genes were divided into five clusters according to their phylogenetic relationships. The intracellular localizations of six GmMDHs were determined by confocal microscopy on Arabidopsis mesophyll protoplasts. Transcripts of GmMDHs were significantly increased by abiotic stress (drought, salt, and alkalinity) and hormone treatments, as shown by analysis of cis-regulatory elements and quantitative real-time polymerase chain reaction (qRT-PCR). GmMDHs displayed unique expression patterns in diverse soybean tissues. It is noteworthy that under salt stress, the expression levels of a chloroplast isoform (GmMDH2) were unusually high, presumably indicating a critical role in soybean responses to salinity. Expression of GmMDH2 in Escherichia coli showed that the recombinant enzyme had NADP-dependent MDH activity. The redox states of the nicotinamide adenine dinucleotide phosphate (NADPH) pool and antioxidant activities were shown to be modulated by GmMDH2 gene overexpression, which in turn reduced reactive oxygen species (ROS) formation in transgenic soybean, significantly enhancing the salt stress resistance. Gene-based association analysis showed that variations in GmMDH2 were strongly linked to seedling salt tolerance. A polymorphism possibly associated with salt tolerance was discovered in the promoter region of GmMDH2. These findings not only improve our understanding of the stress response mechanism by identifying and characterizing the MDH gene family throughout the soybean genome but it also identified a potential candidate gene for the future enchancement of salt tolerance in the soybean.

Sparse spectral graph analysis and its application to gastric cancer drug resistance-specific molecular interplays identification.

Uncovering acquired drug resistance mechanisms has garnered considerable attention as drug resistance leads to treatment failure and death in patients with cancer. Although several bioinformatics studies developed various computational methodologies to uncover the drug resistance mechanisms in cancer chemotherapy, most studies were based on individual or differential gene expression analysis. However the single gene-based analysis is not enough, because perturbations in complex molecular networks are involved in anti-cancer drug resistance mechanisms. The main goal of this study is to reveal crucial molecular interplay that plays key roles in mechanism underlying acquired gastric cancer drug resistance. To uncover the mechanism and molecular characteristics of drug resistance, we propose a novel computational strategy that identified the differentially regulated gene networks. Our method measures dissimilarity of networks based on the eigenvalues of the Laplacian matrix. Especially, our strategy determined the networks' eigenstructure based on sparse eigen loadings, thus, the only crucial features to describe the graph structure are involved in the eigenanalysis without noise disturbance. We incorporated the network biology knowledge into eigenanalysis based on the network-constrained regularization. Therefore, we can achieve a biologically reliable interpretation of the differentially regulated gene network identification. Monte Carlo simulations show the outstanding performances of the proposed methodology for differentially regulated gene network identification. We applied our strategy to gastric cancer drug-resistant-specific molecular interplays and related markers. The identified drug resistance markers are verified through the literature. Our results suggest that the suppression and/or induction of COL4A1, PXDN and TGFBI and their molecular interplays enriched in the Extracellular-related pathways may provide crucial clues to enhance the chemosensitivity of gastric cancer. The developed strategy will be a useful tool to identify phenotype-specific molecular characteristics that can provide essential clues to uncover the complex cancer mechanism.

Oxomer- and Reporter Gene-Based Analysis of FIH Activity in Cells.

Cellular and tissue adaptations to oxygen deprivation (hypoxia) are necessary for both normal physiology and disease. Responses to hypoxia are initiated by the cellular oxygen sensors prolyl-4-hydroxylase domain (PHD) proteins 1-3 and factor inhibiting HIF (FIH). These enzymes regulate the transcription factor hypoxia-inducible factor (HIF) in a hypoxia-sensitive manner. FIH also regulates proteins outside the HIF pathway, including the deubiquitinase OTUB1. Numerous preclinical analyses have demonstrated thattreatment with HIF hydroxylase inhibitors is beneficial and protective in many hypoxia-associated diseases. However, clinically available HIF hydroxylase inhibitors increase erythropoietin (EPO) gene expression and red blood cell production, which can be detrimental in hypoxia-associated conditions, such as ischemia/reperfusion injury of the heart or chronic inflammation. Our understanding of the relevance of FIH in (patho)physiology is only in its infancy, but FIH activity does not govern erythropoietin expression. Therefore, it is of prime interest to assess the relevance of FIH activity in (patho)physiology in detail, as it may contribute to developing novel therapeutic options for treating hypoxia-associated diseases that do not affect Epo regulation. Here, we describe specific protocols for two different methods to assess FIH enzymatic activity within cells, using a HIF-dependent firefly luciferase-reporter gene and an oxomer-dependent assay. Oxomers are oxygen-dependent stable protein oligomers formed by FIH, for example, with the deubiquitinase OTUB1. Oxomer formation directly depends on FIH activity, providing a suitable cellular readout for an easy-to-use analysis of FIH enzymatic activity in cellulo. These techniques permit an analysis of FIH activity toward HIF and outside the HIF pathway, allowing the investigation of FIH activity under different (patho)physiological conditions and assessment of novel (putative) inhibitors.

Dissecting the shared genetic architecture between Alzheimer's disease and frailty: a cross-trait meta-analyses of genome-wide association studies.

Introduction: Frailty is the most common medical condition affecting the aging population, and its prevalence increases in the population aged 65 or more. Frailty is commonly diagnosed using the frailty index (FI) or frailty phenotype (FP) assessments. Observational studies have indicated the association of frailty with Alzheimer's disease (AD). However, the shared genetic and biological mechanism of these comorbidity has not been studied. Methods: To assess the genetic relationship between AD and frailty, we examined it at single nucleotide polymorphism (SNP), gene, and pathway levels. Results: Overall, 16 genome-wide significant loci (15 unique loci) (p meta-analysis < 5 × 10-8) and 22 genes (21 unique genes) were identified between AD and frailty using cross-trait meta-analysis. The 8 shared loci implicated 11 genes: CLRN1-AS1, CRHR1, FERMT2, GRK4, LINC01929, LRFN2, MADD, RP11-368P15.1, RP11-166N6.2, RNA5SP459, and ZNF652 between AD and FI, and 8 shared loci between AD and FFS implicated 11 genes: AFF3, C1QTNF4, CLEC16A, FAM180B, FBXL19, GRK4, LINC01104, MAD1L1, RGS12, ZDHHC5, and ZNF521. The loci 4p16.3 (GRK4) was identified in both meta-analyses. The colocalization analysis supported the results of our meta-analysis in these loci. The gene-based analysis revealed 80 genes between AD and frailty, and 4 genes were initially identified in our meta-analyses: C1QTNF4, CRHR1, MAD1L1, and RGS12. The pathway analysis showed enrichment for lipoprotein particle plasma, amyloid fibril formation, protein kinase regulator, and tau protein binding. Conclusion: Overall, our results provide new insights into the genetics of AD and frailty, suggesting the existence of non-causal shared genetic mechanisms between these conditions.

Side-Stream Enhanced Biological Phosphorus Removal (S2EBPR) enables effective phosphorus removal in a pilot-scale A-B stage shortcut nitrogen removal system for mainstream municipal wastewater treatment

While the adsorption/bio-oxidation (A/B) process has been widely studied for carbon capture and shortcut nitrogen (N) removal, its integration with enhanced biological phosphorus (P) removal (EBPR) has been considered challenging and thus unexplored. Here, full-scale pilot testing with an integrated system combining A-stage high-rate activated sludge with B-stage partial (de)nitrification/anammox and side-stream EBPR (HRAS-P(D)N/A-S2EBPR) was conducted treating real municipal wastewater. The results demonstrated that, despite the relatively low influent carbon load, the B-stage P(D)N-S2EBPR system could achieve effective P removal performance, with the carbon supplement and redirection of the A-stage sludge fermentate to the S2EBPR. The novel process configuration design enabled a system shift in carbon flux and distribution for efficient EBPR, and provided unique selective factors for ecological niche partitioning among different key functionally relevant microorganisms including polyphosphate accumulating organisms (PAOs) and glycogen-accumulating organisms (GAOs). The combined nitrite from B-stage to S2EBPR and aerobic-anoxic conditions in our HRAS-P(D)N/A-S2EBPR system promoted DPAOs for simultaneous internal carbon-driven denitrification via nitrite and P removal. 16S rRNA gene-based oligotyping analysis revealed high phylogenetic microdiversity within the Accumulibacter population and discovered coexistence of certain oligotypes of Accumulibacter and Competibacter correlated with efficient P removal. Single-cell Raman micro-spectroscopy-based phenotypic profiling showed high phenotypic microdiversity in the active PAO community and the involvement of unidentified PAOs and internal carbon-accumulating organisms that potentially played an important role in system performance. This is the first pilot study to demonstrate that the P(D)N-S2EBPR system could achieve shortcut N removal and influent carbon-independent EBPR simultaneously, and the results provided insights into the effects of incorporating S2EBPR into A/B process on metabolic activities, microbial ecology, and resulted system performance.

Toll receptor ligand Spätzle 4 responses to the highly pathogenic Enterococcus faecalis from Varroa mites in honeybees.

Honeybees play a major role in crop pollination, which supports the agricultural economy and international food supply. The colony health of honeybees is threatened by the parasitic mite Varroa destructor, which inflicts physical injury on the hosts and serves as the vector for variable viruses. Recently, it shows that V. destructor may also transmit bacteria through the feeding wound, yet it remains unclear whether the invading bacteria can exhibit pathogenicity to the honeybees. Here, we incidentally isolate Enterococcus faecalis, one of the most abundant bacteria in Varroa mites, from dead bees during our routine generation of microbiota-free bees in the lab. In vivo tests show that E. faecalis is only pathogenic in Apis mellifera but not in Apis cerana. The expression of antimicrobial peptide genes is elevated following infection in A. cerana. The gene-based molecular evolution analysis identifies positive selection of genes encoding Späetzle 4 (Spz4) in A. cerana, a signaling protein in the Toll pathway. The amino acid sites under positive selection are related to structural changes in Spz4 protein, suggesting improvement of immunity in A. cerana. The knock-down of Spz4 in A. cerana significantly reduces the survival rates under E. faecalis challenge and the expression of antimicrobial peptide genes. Our results indicate that bacteria associated with Varroa mites are pathogenic to adult bees, and the positively selected gene Spz4 in A. cerana is crucial in response to this mite-related pathogen.

Not all parents are the same: Diverse strategies of symbiont transmission in seaweeds.

Different marine seaweed species have been shown to harbour specific bacterial communities, however, the extent to which vertical symbiont transmission from parents to offspring contributes to host-specificity is unclear. Here we use fluorescence and electron microscopy as well as 16S rRNA gene-based community analysis to investigate symbiont transmission in members of the three major seaweed groups (green Chlorophyta, red Rhodophyta and brown Phaeophyceae). We found seaweeds employ diverse strategies to transfer symbionts to their progeny. For instance, the green Ulva australis does not appear to have the capacity for vertical transmission. In contrast, the brown Phyllospora comosa adopts a non-selective vertical transmission. The red Delisea pulchra demonstrates weak selectivity in symbiont transmission, while the brown Hormosira banksii exhibits a strongly selective symbiont transfer. Mucilage on the gametes appears to facilitate vertical transmission and transferred bacteria have predicted properties that could support early development of the seaweeds. Previous meta-analysis has indicated that vertical transmission is rare in aquatic compared to terrestrial environments, however, our results contribute to the growing evidence that this might not be the case and that instead vertical transmission with various degrees of symbiont selection occurs in the ecologically important group of seaweeds.

Identification of novel rare variants for anxiety: an exome-wide association study in the UK Biobank

BackgroundRare variants are believed to play a substantial role in the genetic architecture of mental disorders, particularly in coding regions. However, limited evidence supports the impact of rare variants on anxiety. MethodsUsing whole-exome sequencing data from 200,643 participants in the UK Biobank, we investigated the contribution of rare variants to anxiety. Firstly, we computed genetic risk score (GRS) of anxiety utilizing genotype data and summary data from a genome-wide association study (GWAS) on anxiety disorder. Subsequently, we identified individuals within the lowest 50% GRS, a subgroup more likely to carry pathogenic rare variants. Within this subgroup, we classified individuals with the highest 10% 7-item Generalized Anxiety Disorder scale (GAD-7) score as cases (N = 1869), and those with the lowest 10% GAD-7 score were designated as controls (N = 1869). Finally, we conducted gene-based burden tests and single-variant association analyses to assess the relationship between rare variants and anxiety. ResultsTotally, 47,800 variants with MAF ≤0.01 were annotated as non-benign coding variants, consisting of 42,698 nonsynonymous SNVs, 489 nonframeshift substitution, 236 frameshift substitution, 617 stop-gain and 40 stop-loss variants. After variation aggregation, 5066 genes were included in gene-based association analysis. Totally, 11 candidate genes were detected in burden test, such as RNF123 (PBonferroni adjusted = 3.40 × 10−6), MOAP1(PBonferroni adjusted = 4.35 × 10−4), CCDC110 (PBonferroni adjusted = 5.83 × 10−4). Single-variant test detected 9 rare variants, such as rs35726701(RNF123)(PBonferroni adjusted = 3.16 × 10−10) and rs16942615(CAMTA2) (PBonferroni adjusted = 4.04 × 10−4). Notably, RNF123, CCDC110, DNAH2, and CSKMT gene were identified in both tests. ConclusionsOur study identified novel candidate genes for anxiety in protein-coding regions, revealing the contribution of rare variants to anxiety.

Commensal Fecal Microbiota Profiles Associated with Initial Stages of Intestinal Mucosa Damage: A Pilot Study.

Progressive intestinal mucosal damage occurs over years prior to colorectal cancer (CRC) development. The endoscopic screening of polyps and histopathological examination are used clinically to determine the risk and progression of mucosal lesions. We analyzed fecal microbiota compositions using 16S rRNA gene-based metataxonomic analyses and the levels of short-chain fatty acids (SCFAs) using gas chromatography in volunteers undergoing colonoscopy and histopathological analyses to determine the microbiota shifts occurring at the early stages of intestinal mucosa alterations. The results were compared between diagnosis groups (nonpathological controls and polyps), between samples from individuals with hyperplastic polyps or conventional adenomas, and between grades of dysplasia in conventional adenomas. Some microbial taxa from the Bacillota and Euryarchaeota phyla were the most affected when comparing the diagnosis and histopathological groups. Deeper microbiota alterations were found in the conventional adenomas than in the hyperplastic polyps. The Ruminococcus torques group was enriched in both the hyperplastic polyps and conventional adenomas, whereas the family Eggerthellaceae was enriched only in the hyperplastic polyps. The abundance of Prevotellaceae, Oscillospiraceae, Methanobacteriaceae, Streptococcaceae, Christensenellaceae, Erysipelotrichaceae, and Clostridiaceae shifted in conventional adenomas depending on the grade of dysplasia, without affecting the major SCFAs. Our results suggest a reorganization of microbial consortia involved in gut fermentative processes.

UspA gene-based phylogenetic analysis and antigenic epitope prediction for Escherichia coli strains of avian origin.

Pathogenic Escherichia coli (E. coli) is responsible for various local and systemic infections in animal and human populations. Conventional methods for the detection and identification of E. coli are time-consuming and less reliable for atypical strains. The uspA gene has been widely used as a target for the detection of E. coli. The present study was aimed at phylogenetic analysis of the uspA gene sequences to determine the evolutionary relationships between the strains and other members of the Enterobacteriaceae family. In addition, the unique differences in the sequences of the current study with Salmonella and Shigella species were tested using Tajima's molecular clock test. Antigenic epitope prediction was performed to locate the B-cell epitope region of the UspA protein. Two E. coli isolates of avian origin and strains from the National Center for Biotechnology Information (NCBI) database were used for prediction. The Immune Epitope Database (IEDB) server, Bepitope, ABCpred, SVMTrip, and ElliPro server were used to identify B-cell epitopes. The 3D structure was predicted using SWISS-MODEL. Phylogenetic analysis of the isolates from the current study revealed that both OM837340 and OM837341 sequences from the current study had maximum nucleotide homology (nt) of 99.87%-100% with E. coli isolates and minimum nt homology of 84.08% with Salmonella enteritidis and S. Hissar. The isolates in the current study had a homology of 98.87%, while the homology with Shigella species was 99.25%. Seven silent mutations were observed in the coding region of the UspA protein of ECO9LTBW (current study). Modeling of the UspA protein revealed a maximum homology of 67.86% with the Protein Data Bank in Europe (PDBe), also validated by the Ramachandran plot. No significant differences were found in the coding regions of uspA of Salmonella, Shigella, and E. coli with Tajima's test. For the E. coli isolates, a total of 24 linear B-cell and seven discontinuous epitopes were predicted using in-silico analysis. When the results of the predicted peptides were compared, two peptides, namely ARPYNA and YSDLYTGLIDVNLGDMQKRISEE, were found suitable candidates. In conclusion, the uspA gene appears to be conserved among E. coli isolates and can be used for molecular detection.

The first exome wide association study in Tunisia: identification of candidate loci and pathways with biological relevance for type 2 diabetes.

Type 2 diabetes (T2D) is a multifactorial disease involving genetic and environmental components. Several genome-wide association studies (GWAS) have been conducted to decipher potential genetic aberrations promoting the onset of this metabolic disorder. These GWAS have identified over 400 associated variants, mostly in the intronic or intergenic regions. Recently, a growing number of exome genotyping or exome sequencing experiments have identified coding variants associated with T2D. Such studies were mainly conducted in European populations, and the few candidate-gene replication studies in North African populations revealed inconsistent results. In the present study, we aimed to discover the coding genetic etiology of T2D in the Tunisian population. We carried out a pilot Exome Wide Association Study (EWAS) on 50 Tunisian individuals. Single variant analysis was performed as implemented in PLINK on potentially deleterious coding variants. Subsequently, we applied gene-based and gene-set analyses using MAGMA software to identify genes and pathways associated with T2D. Potential signals were further replicated in an existing large in-silico dataset, involving up to 177116 European individuals. Our analysis revealed, for the first time, promising associations between T2D and variations in MYORG gene, implicated in the skeletal muscle fiber development. Gene-set analysis identified two candidate pathways having nominal associations with T2D in our study samples, namely the positive regulation of neuron apoptotic process and the regulation of mucus secretion. These two pathways are implicated in the neurogenerative alterations and in the inflammatory mechanisms of metabolic diseases. In addition, replication analysis revealed nominal associations of the regulation of beta-cell development and the regulation of peptidase activity pathways with T2D, both in the Tunisian subjects and in the European in-silico dataset. The present study is the first EWAS to investigate the impact of single genetic variants and their aggregate effects on T2D risk in Africa. The promising disease markers, revealed by our pilot EWAS, will promote the understanding of the T2D pathophysiology in North Africa as well as the discovery of potential treatments.

Whole genome sequence analysis of apparent treatment resistant hypertension status in participants from the Trans-Omics for Precision Medicine program.

Introduction: Apparent treatment-resistant hypertension (aTRH) is characterized by the use of four or more antihypertensive (AHT) classes to achieve blood pressure (BP) control. In the current study, we conducted single-variant and gene-based analyses of aTRH among individuals from 12 Trans-Omics for Precision Medicine cohorts with whole-genome sequencing data. Methods: Cases were defined as individuals treated for hypertension (HTN) taking three different AHT classes, with average systolic BP ≥ 140 or diastolic BP ≥ 90mmHg, or four or more medications regardless of BP (n = 1,705). A normotensive control group was defined as individuals with BP < 140/90mmHg (n = 22,079), not on AHT medication. A second control group comprised individuals who were treatment responsive on one AHT medication with BP < 140/ 90mmHg (n = 5,424). Logistic regression with kinship adjustment using the Scalable and Accurate Implementation of Generalized mixed models (SAIGE) was performed, adjusting for age, sex, and genetic ancestry. We assessed variants using SKAT-O in rare-variant analyses. Single-variant and gene-based tests were conducted in a pooled multi-ethnicity stratum, as well as self-reported ethnic/racial strata (European and African American). Results: One variant in the known HTN locus, KCNK3, was a top finding in the multi-ethnic analysis (p = 8.23E-07) for the normotensive control group [rs12476527, odds ratio (95% confidence interval) = 0.80 (0.74-0.88)]. This variant was replicated in the Vanderbilt University Medical Center's DNA repository data. Aggregate gene-based signals included the genes AGTPBP, MYL4, PDCD4, BBS9, ERG, and IER3. Discussion: Additional work validating these loci in larger, more diverse populations, is warranted to determine whether these regions influence the pathobiology of aTRH.

Identification of Superior Haplotypes and Haplotype Combinations for Grain Size- and Weight-Related Genes for Breeding Applications in Rice (Oryza sativa L.).

The identification of superior haplotypes and haplotype combinations is essential for haplotype-based breeding (HBB), which provides selection targets for genomics-assisted breeding. In this study, genotypes of 42 functional genes in rice were analyzed by targeted capture sequencing in a panel of 180 Indica rice accessions. In total, 69 SNPs/Indels in seven genes were detected to be associated with grain length (GL), grain width (GW), ratio of grain length-width (L/W) and thousand-grain weight (TGW) using candidate gene-based association analysis, including BG1 and GS3 for GL, GW5 for GW, BG1 and GW5 for L/W, and AET1, SNAC1, qTGW3, DHD1 and GW5 for TGW. Furthermore, two haplotypes were identified for each of the seven genes according to these associated SNPs/Indels, and the amount of genetic variation explained by different haplotypes ranged from 3.24% to 27.66%. Additionally, three, three and eight haplotype combinations for GL, L/W and TGW explained 25.38%, 5.5% and 22.49% of the total genetic variation for each trait, respectively. Further analysis showed that Minghui63 had the superior haplotype combination Haplotype Combination 4 (HC4) for TGW. The most interesting finding was that some widely used restorer lines derived from Minghui63 also have the superior haplotype combination HC4, and our breeding varieties and lines using the haplotype-specific marker panel also confirmed that the TGW of the lines was much higher than that of their sister lines without HC4, suggesting that TGW-HC4 is the superior haplotype combination for TGW and can be utilized in rice breeding.

Sex bias in celiac disease: XWAS and monocyte eQTLs in women identify TMEM187 as a functional candidate gene

BackgroundCeliac disease (CeD) is an immune-mediated disorder that develops in genetically predisposed individuals upon gluten consumption. HLA risk alleles explain 40% of the genetic component of CeD, so there have been continuing efforts to uncover non-HLA loci that can explain the remaining heritability. As in most autoimmune disorders, the prevalence of CeD is significantly higher in women. Here, we investigated the possible involvement of the X chromosome on the sex bias of CeD.MethodsWe performed a X chromosome-wide association study (XWAS) and a gene-based association study in women from the CeD Immunochip (7062 cases, 5446 controls). We also constructed a database of X chromosome cis-expression quantitative trait loci (eQTLs) in monocytes from unstimulated (n = 226) and lipopolysaccharide (LPS)-stimulated (n = 130) female donors and performed a Summary-data-based MR (SMR) analysis to integrate XWAS and eQTL information. We interrogated the expression of the potentially causal gene (TMEM187) in peripheral blood mononuclear cells (PBMCs) from celiac patients at onset, on a gluten-free diet, potential celiac patients and non-celiac controls.ResultsThe XWAS and gene-based analyses identified 13 SNPs and 25 genes, respectively, 22 of which had not been previously associated with CeD. The X chromosome cis-eQTL analysis found 18 genes with at least one cis-eQTL in naïve female monocytes and 8 genes in LPS-stimulated female monocytes, 2 of which were common to both situations and 6 were unique to LPS stimulation. SMR identified a potentially causal association of TMEM187 expression in naïve monocytes with CeD in women, regulated by CeD-associated, eQTL-SNPs rs7350355 and rs5945386. The CeD-risk alleles were correlated with lower TMEM187 expression. These results were replicated using eQTLs from LPS-stimulated monocytes. We observed higher levels of TMEM187 expression in PBMCs from female CeD patients at onset compared to female non-celiac controls, but not in male CeD individuals.ConclusionUsing X chromosome genotypes and gene expression data from female monocytes, SMR has identified TMEM187 as a potentially causal candidate in CeD. Further studies are needed to understand the implication of the X chromosome in the higher prevalence of CeD in women.

Abstract B069: Lung cancer risk in women and racial minorities and genetic variants of UGT genes

Abstract The major research focus of this project is the identification of genetic variants involved in disparities in carcinogen metabolism and risk for cancer. Specifically, this research focuses on the identification of genetic variants in carcinogen metabolism genes, that increase lung cancer risk, particularity when the risk differs between men and women or by racial groups. Although cigarette smoking causes the vast majority of lung cancers, most lifetime smokers do not develop this disease. Many studies suggest that women have greater susceptibility than men for lung cancer, and that risk varies by racial groups. Genetic variation in tobacco carcinogen-metabolizing enzymes could give details as to why only some smokers develop lung cancer and why women or certain racial groups may be at a greater risk. Genetic changes have been shown to vary significantly between individuals of different races, which could contribute to racial disparities in cancer risk. Uridine glucuronosyltransferase (UGT) genes detoxify many carcinogens found in tobacco, so variants in UGT genes could explain differences in lung cancer risk among smokers. Many genetic association studies of UGT genes and lung cancer risk have been performed to identify alleles (SNPs and Copy Number Variants) in UGT genes that increase risk for lung cancer, and to determine if some variants display a gender or racial-specific association with this cancer. cBioPortal provides gene-based visualizations and analyses, so that users can find altered genes and/or networks within a study of interest or across all cancer studies. cBioPortal is a web application platform that facilitates data exploration and analysis through the use of a variety of visualization and analytical tools using populations of cancer cases and controls. cBioPortal has been utilized in this study to visualize variations in UGT genes across many lung cancer data sets. A copy number variant that deletes the entire UGT2B17 gene (which metabolizes the lung cancer carcinogen, NNAL) is a particular focus in this study using data from cBioPortal, as well as common copy number alterations (amplifications) in UGT3A1 and UGT3A2. Citation Format: James W. Gallagher, Talayshia Vicks, Carla J. Gallagher. Lung cancer risk in women and racial minorities and genetic variants of UGT genes [abstract]. In: Proceedings of the 16th AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2023 Sep 29-Oct 2;Orlando, FL. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2023;32(12 Suppl):Abstract nr B069.

Type 2 Diabetes Modifies the Association of CAD Genomic Risk Variants With Subclinical Atherosclerosis.

Individuals with type 2 diabetes (T2D) have an increased risk of coronary artery disease (CAD), but questions remain about the underlying pathology. Identifying which CAD loci are modified by T2D in the development of subclinical atherosclerosis (coronary artery calcification [CAC], carotid intima-media thickness, or carotid plaque) may improve our understanding of the mechanisms leading to the increased CAD in T2D. We compared the common and rare variant associations of known CAD loci from the literature on CAC, carotid intima-media thickness, and carotid plaque in up to 29 670 participants, including up to 24 157 normoglycemic controls and 5513 T2D cases leveraging whole-genome sequencing data from the Trans-Omics for Precision Medicine program. We included first-order T2D interaction terms in each model to determine whether CAD loci were modified by T2D. The genetic main and interaction effects were assessed using a joint test to determine whether a CAD variant, or gene-based rare variant set, was associated with the respective subclinical atherosclerosis measures and then further determined whether these loci had a significant interaction test. Using a Bonferroni-corrected significance threshold of P<1.6×10-4, we identified 3 genes (ATP1B1, ARVCF, and LIPG) associated with CAC and 2 genes (ABCG8 and EIF2B2) associated with carotid intima-media thickness and carotid plaque, respectively, through gene-based rare variant set analysis. Both ATP1B1 and ARVCF also had significantly different associations for CAC in T2D cases versus controls. No significant interaction tests were identified through the candidate single-variant analysis. These results highlight T2D as an important modifier of rare variant associations in CAD loci with CAC.

Polyphasic Characterization and Genomic Insights into an Aerobic Denitrifying Bacterium, Shewanella zhuhaiensis sp. nov., Isolated from a Tidal Flat Sediment.

A new, facultatively anaerobic, light-yellow, and rod-shaped bacterium designated as 3B26T isolated from Qi'ao Island's tidal flat sediment was identified. Strain 3B26T can hydrolyze gelatin, aesculin, and skim milk. The major cellular fatty acids were identified as iso-C15:0, referred to as summed feature 3, and C16:0; the polar lipids consisted of phosphatidylethanolamine, phosphatidylglycerol, and phospholipid; and the quinones contained Q-7, Q-8, MK-7, and MMK7. The genomic size of strain 3B26T was 4,682,650 bp, and its genomic DNA G + C content was 54.8%. While a 16S rRNA gene-based phylogenetic analysis confirmed that strain 3B26T belongs to the genus Shewanella, both phylogenomic inference and genomic comparison revealed that strain 3B26T is distinguishable from its relatives, and digital DNA-DNA hybridization (dDDH) values of 24.4-62.6% and average nucleotide identities (ANIs) of 83.5-95.6% between them were below the 70% dDDH and 96% ANI thresholds for bacterial species delineation. Genomic functional analysis demonstrated that strain 3B26T possesses complete gene clusters of eicosapentaenoic acid biosynthesis and denitrification. Based on the evidence above, strain 3B26T is considered to represent a novel species of the genus Shewanella, and the name Shewanella zhuhaiensis sp. nov. (type strain 3B26T = GDMCC 1.2057T = KCTC 82339T) is proposed.