Genome-wide Approach Research Articles

Genetic and structural basis of colistin resistance in Klebsiella pneumoniae: Unravelling the molecular mechanisms

ObjectiveAntimicrobial resistance (AMR), together with multidrug resistance (MDR), mainly among Gram-negative bacteria, has been on the rise. Colistin (polymyxin E) remains one of the primary available last resorts to treat infections caused by MDR bacteria during the rapid emergence of global resistance. As the exact mechanism of bacterial resistance to colistin remains undetermined, this study warranted elucidation of the underlying mechanisms of colistin resistance and heteroresistance among carbapenem-resistant Klebsiella pneumoniae isolates. MethodsMolecular analysis was carried out on the resistant isolates using a genome-wide characterisation approach, as well as MALDI-TOF mass spectrometry, to identify lipid A. ResultsAmong the 32 carbapenem-resistant K. pneumoniae isolates, several isolates showed resistance and intermediate resistance to colistin. The seven isolates with intermediate resistance exhibited the “skip-well” phenomenon, attributed to the presence of resistant subpopulations. The three isolates with full resistance to colistin showed ions using MALDI-TOF mass spectrometry at m/z of 1840 and 1824 representing bisphosphorylated and hexa-acylated lipid A, respectively, with or without hydroxylation at position C’-2 of the fatty acyl chain. Studying the genetic environment of mgrB locus revealed the presence of two insertion sequences that disrupted the mgrB locus in the three colistin-resistant isolates: IS1R and IS903B. ConclusionsOur findings show that colistin resistance/heteroresistance was inducible with mutations in chromosomal regulatory networks controlling the lipid A moiety and insertion sequences disrupting the mgrB gene, leading to elevated minimum inhibitory concentration values and treatment failure. Different treatment strategies should be employed to avoid colistin heteroresistance-linked treatment failures, mainly through combination therapy using colistin with carbapenems, aminoglycosides, or tigecycline.

Semaphorin 6D tunes amygdalar circuits for emotional, metabolic, and inflammatory outputs

Regulated neural-metabolic-inflammatory responses are essential for maintaining physiological homeostasis. However, the molecular machinery that coordinates neural, metabolic, and inflammatory responses is largely unknown. Here, we show that semaphorin 6D (SEMA6D) coordinates anxiogenic, metabolic, and inflammatory outputs from the amygdala by maintaining synaptic homeostasis. Using genome-wide approaches, we identify SEMA6D as a pleiotropic gene for both psychiatric and metabolic traits in human. Sema6d deficiency increases anxiety in mice. When fed a high-fat diet, Sema6d-/- mice display attenuated obesity and enhanced myelopoiesis compared with control mice due to higher sympathetic activity via the β3-adrenergic receptor. Genetic manipulation and spatial and single-nucleus transcriptomics reveal that SEMA6D in amygdalar interneurons is responsible for regulating anxiogenic and autonomic responses. Mechanistically, SEMA6D is required for synaptic maturation and γ-aminobutyric acid transmission. These results demonstrate that SEMA6D is important for the normal functioning of the neural circuits in the amygdala, coupling emotional, metabolic, and inflammatory responses.

NucMACC: An MNase-seq pipeline to identify structurally altered nucleosomes in the genome.

Micrococcal nuclease sequencing is the state-of-the-art method for determining chromatin structure and nucleosome positioning. Data analysis is complex due to the AT-dependent sequence bias of the endonuclease and the requirement for high sequencing depth. Here, we present the nucleosome-based MNase accessibility (nucMACC) pipeline unveiling the regulatory chromatin landscape by measuring nucleosome accessibility and stability. The nucMACC pipeline represents a systematic and genome-wide approach for detecting unstable ("fragile") nucleosomes. We have characterized the regulatory nucleosome landscape in Drosophila melanogaster, Saccharomyces cerevisiae, and mammals. Two functionally distinct sets of promoters were characterized, one associated with an unstable nucleosome and the other being nucleosome depleted. We show that unstable nucleosomes present intermediate states of nucleosome remodeling, preparing inducible genes for transcriptional activation in response to stimuli or stress. The presence of unstable nucleosomes correlates with RNA polymerase II proximal pausing. The nucMACC pipeline offers unparalleled precision and depth in nucleosome research and is a valuable tool for future nucleosome studies.

P-284 Large-scale genome-wide association study (GWAS) to determine the genetic underpinnings of female genital tract polyps (FGT)

Abstract Study question Can a large-scale GWAS meta-analysis identify the genomic risk loci and associated candidate genes for FGT polyps? Summary answer We identified ten genomic risk loci, with two exonic missense variants in PRIM1 and COL17A1 genes, highlighting the shared mechanisms with cancerous processes. What is known already Small-scale candidate gene studies have focused on biological processes such as oestrogen stimulation and inflammation to clarify the biology behind FGT polyps. However, the exact mechanism for the development of polyps is still elusive. At the same time, a genome-wide approach, which has become the gold standard in complex disease genetics, has never been used to uncover the genetics of the FGT polyps. Study design, size, duration We performed a genome wide association study (GWAS) meta-analysis including a total of 25,100 women with FGT polyps (International Classification of Disease, ICD-10 diagnosis code N84) and 207,193 female controls (without N84 code) of European ancestry from the FinnGen study (11,092 cases and 94,394 controls) and the Estonian Biobank (EstBB, 14,008 cases and 112,799 controls). Participants/materials, setting, methods A meta-analysis and functional annotation of GWAS signals were performed using the Functional mapping and annotation of genetic associations (FUMA) tool. Gene prioritization for the identified variants was conducted using the variant-to-gene (V2G) score from the OpenTarget genetics portal, followed by FUMA analysis. To determine associations with other phenotypes, we performed a look-up of associated variants across multiple traits, a genetic correlation analysis, and a phenome-wide association study (PheWAS) with ICD10 diagnosis codes. Main results and the role of chance Our GWAS meta-analysis revealed ten significant (P < 5 x 10-8) genomic risk loci. Two signals, rs2277339 (P = 7.6 x 10-10) and rs1265005 (P = 1.1 x 10-9) (in linkage disequilibrium (LD) with rs805698 r2 = 0.75), are exonic missense variants in PRIM1, and COL17A1 genes, respectively. Based on the literature, these genes may play a role in cellular proliferation. Several of the identified genomic loci had previously been linked to endometrial cancer and/or uterine fibroids. Thus, highlighting the shared mechanisms underlying tissue overgrowth and cancerous processes, which may be relevant to the development of polyps. Genetic correlation analysis revealed a negative correlation between sex hormone-binding globulin (SHBG) and the risk of FGT polyps (rg = -0,21, se = 0.04, P = 2.9 x 10-6), and on the phenotypic level (PheWAS), the strongest associations were observed with endometriosis, leiomyoma of the uterus and excessive, frequent and irregular menstruation. Limitations, reasons for caution In this study, we focused broadly on FGT polyps and did not differentiate between the polyp subtypes. The prevalence of FGT polyps led us to assume that most women included in the study had endometrial polyps. Wider implications of the findings The study findings have the potential to significantly enhance our understanding of the genetic mechanisms involved, paving the way for future functional follow-up, which in turn could improve the diagnosis, risk assessment, and targeted treatment options, since surgery is the only line of treatment available for diagnosed polyps. Trial registration number not applicable

Characterization of the GGP gene family in potato (Solanum tuberosum L.) and Pepper (Capsicum annuum L.) and its expression analysis under hormonal and abiotic stresses

GDP-l-galactose phosphorylase (GGP) is a key rate-limiting enzyme in plant ascorbic acid synthesis, which plays an important role in plant growth and development as well as stress response. However, the presence of GGP and its function in potato and pepper are not known. In this study, we first identified two GGP genes in each potato and pepper genomes using a genome-wide search approach. We then analyzed their physicochemical properties, conserved domains, protein structures and phylogenetic relationships. Phylogenetic tree analysis revealed that members of the potato and pepper GGP gene families are related to eggplant (Solanum melongena L.), Arabidopsis (Arabidopsis thaliana L.), tobacco (Nicotiana tabacum L.) and tomato (Solanum lycopersicum L.), with tomato being the most closely related. The promoter sequences mainly contain homeopathic elements such as light-responsive, hormone-responsive and stress-responsive, with light-responsive elements being the most abundant. By analyzing the structure of the genes, it was found that there is no transmembrane structure or signal peptide in the GGP gene family of potatoes and peppers, and that all of its members are hydrophilic proteins. The expression profiles of different tissues show that StGGP1 has the highest expression levels in leaves, StGGP2 has the highest expression levels in stamens, and CaGGPs have the highest expression levels in the early stages of fruit development (Dev1). It was found that StGGPs and CaGGPs genes showed different response to phytohormones and abiotic stresses. Abscisic acid (ABA) treatment induced the most significant change in the expression of StGGPs, while the expression of CaGGPs showed the most pronounced change under methyl jasmonate (MeJA) treatment. StGGPs responded mainly to dark treatment, whereas CaGGPs responded mainly to NaCl stress. These results provide an important basis for a detailed study about the functions of GGP homologous genes in potato and pepper in response to abiotic stresses.

P-521 Inter-individual variation in gonadotoxicity in female childhood cancer survivors – a genome-wide association study: results from the PanCareLIFE study

Abstract Study question To identify new variants associated with alkylating agent (AA) induced gonadotoxicity among adult female childhood cancer survivors (CCS) using a genome-wide association study approach (GWAS). Summary answer Upon meta-analysis of two independent cohorts, the variant rs78861946 was associated with lower AMH levels in CCS, but not separately in the smaller replication cohort. What is known already Female childhood cancer survivors are at risk of therapy-related gonadal impairment. Alkylating agents are well-established risk factors, and the inter-individual variability in gonadotoxicity may be explained by genetic polymorphisms. Study design, size, duration A genome-wide association study (GWAS) approach was used to discover new variants associated with AA-induced gonadotoxicity among adult female CCS. The discovery cohort included adult female CCS from the pan-European PanCareLIFE cohort (n = 743; median age: 25.8 years), excluding those who received bilateral ovarian irradiation, bilateral oophorectomy, CNS, total body irradiation, or stem cell transplantation. Results were replicated in the USA-based St. Jude Lifetime Cohort (n = 391; median age: 31.3 years) and combined in a meta-analysis. Participants/materials, setting, methods Anti-Müllerian hormone (AMH) levels served as a proxy for ovarian function in both the discovery and replication cohort. A GWAS was performed using a linear regression model adjusted for cyclophosphamide equivalent dose of alkylating agents, age at diagnosis, and age at study. Main results and the role of chance Three genome-wide significant (<5.0x10E-8 ) and 16 genome-wide suggestive (<5.0x10E-6 ) loci were associated with log-transformed AMH levels, adjusted for cyclophosphamide equivalent dose of alkylating agents, age at diagnosis, and age at study in the PanCareLIFE cohort. Based on effect allele frequency (EAF) (>0.01 if not genome-wide significant), p-value (<5.0 × 10E−6 ), and biological relevance, 15 SNPs were selected for replication in the independent replication cohort of the St. Jude Lifetime Cohort. None of the SNPs were statistically significantly associated with AMH levels in this smaller cohort. A meta-analysis indicated thatrs78861946 was associated at borderline genome-wide statistical significance (reference/effect allele: C/T;EAF: 0.04, Beta (SE): −0.484 (0.091), p-value= 9.39 × 10E−8). While the sample size of both the discovery and replication cohorts are large for CCS cohorts, they are very small compared to GWASs in the general population. The limited power is likely one of the reasons the identified SNPs were not statistically significant in our replication cohort of SJLIFE. Limitations, reasons for caution While AMH is a valuable surrogate marker for ovarian reserve, it remains a proxy for ovarian reserve and low AMH levels may still coincide with fertility potential and pregnancy. Wider implications of the findings Identification of risk-contributing variants such as rs78861946 may prove useful to better identify patients who are specifically susceptible for alkylating agent-induced gonadal toxicity. Future polygenic risk scores may aid individualized counselling regarding the treatment-related risk of gonadotoxicity and fertility preservation options for female childhood cancer survivors. Trial registration number not applicable

P-551 Application of simultaneous PGT-M/PGT-A using whole genome sequencing and haplotyping for monogenetic disorders, multiple variants, low penetrance, and de novo disorders

Abstract Study question Is the concurrent use of high-resolution PGT-A and genome-wide haplotyping with disease-specific assays for PGT-M effective in an unselected patient population at-risk for monogenetic disorders? Summary answer PGT-A/PGT-M is an effective treatment for couples/patients burdened with genetic diseases, that provides optimal family planning, prevents vertical transmission, and reduces personal and societal encumbrance What is known already Societal awareness of genetic testing expanded the use of carrier screening, diagnostic genetic tests, and higher in vitro fertilization (IVF) utilization has increased demand for preimplantation genetic testing for monogenetic disorders (PGT-M), as an alternative to invasive prenatal tests. The scope of PGT-M testing expands beyond common severe childhood disorders, to include serious and mild late onset disorders, cancer predisposition, variants of unknown significance, de novo mutations, and variants in multiple genes. This highlights novel ethical dilemmas, necessitating adjustment of generic genome wide approaches to complement, or eventually replace, targeted disease-specific PGT-M assays. Study design, size, duration This is a retrospective study of 359 couples at-risk for monogenetic disorders referred to the CReATe Fertility Centre, Toronto, Canada. Our objective was to evaluate the efficiency of concurrent use of genome wide haplotyping and disease specific assays for PGT-M and high resolution PGT-A in an unselected patient population at risk for monogenetic disorders, multiple and de novo variants, low penetrance and variants of unknown significance. Participants/materials, setting, methods All patients received comprehensive genetic counselling. For each couple/patient, DNA from gamete providers, related family member/s and whole-genome amplified (WGA)DNA from a single trophectoderm(TE) biopsy were analyzed. Whole-genome SNP-array (Karyomapping, Illumina, CA) and STR-analysis were used to obtain parental haplotypes and Sanger-sequencing was used for direct variant analysis. PGT-A was performed using Illumina platform(NextSeq 550). NxClinical (BioDiscovery) and BlueFuse Multi(Illumina) software were used for data analysis. Main results and the role of chance 1826 embryos, from 235 IVF cycles were analyzed for 88 unique genes (14 cancer predisposing, 42 autosomal recessive and 32 dominant genes) and a total of 282 unique variants. The most frequent genes tested were BRCA1/2 -21% of all cases, followed by HBB-8.7%, CFTR-8.7% and CYP21A2-8.7%. PGT-M for 11 cases was done for a de novo variant detected in a sibling and in 2 cases for a parental mosaic variant. Multiple genes and VUS were analyzed in 12 and 6 cases, respectively. PGT-A analysis showed 31.6% aneuploidy, 9.6% mosaicism and 58.8% euploidy. Genome wide haplotyping was used for euploid and mosaic embryos with an average case call rate of 95%, and average case Allele Drop Out (ADO) rate of 2.7%. Recombination events in 3% of cases prevented assignment of haplotype. Direct mutation testing was obtained for all samples and was 100% concordant with the predicted haplotypes. In cases with no predicted haplotype, direct testing established the PGT-M diagnosis. The distribution for both recessive and dominant disorders was in Mendelian order. To date 66 live births confirmed the results from PGT-M diagnosis. Limitations, reasons for caution Ethical considerations and extensive genetic counselling are prerequisites for PGT-A/M testing implementation. Simultaneous screening for embryo aneuploidy and monogenetic disorders optimizes patient treatment and family planning, but may not be available in all jurisdictions. Although the data was collected prospectively, retrospective design is a limitation. Wider implications of the findings Simultaneous PGT-A and PGT-M with direct variant testing provides optimal embryo prioritization for transfer and significantly reduces time to a healthy pregnancy in couples at risk for inherited disorders. Karyomapping is effective for generic haplotyping, with a low risk of undetected recombination events in the genomic region of interest. Trial registration number not applicable

Mono-methylation of lysine 27 at histone 3 confers lifelong susceptibility to stress

Histone post-translational modifications are critical for mediating persistent alterations in gene expression. By combining unbiased proteomics profiling and genome-wide approaches, we uncovered a role for mono-methylation of lysine 27 at histone H3 (H3K27me1) in the enduring effects of stress. Specifically, mice susceptible to early life stress (ELS) or chronic social defeat stress (CSDS) displayed increased H3K27me1 enrichment in the nucleus accumbens (NAc), a key brain-reward region. Stress-induced H3K27me1 accumulation occurred at genes that control neuronal excitability and was mediated by the VEFS domain of SUZ12, a core subunit of the polycomb repressive complex-2, which controls H3K27 methylation patterns. Viral VEFS expression changed the transcriptional profile of the NAc, led to social, emotional, and cognitive abnormalities, and altered excitability and synaptic transmission of NAc D1-medium spiny neurons. Together, we describe a novel function of H3K27me1 in the brain and demonstrate its role as a "chromatin scar" that mediates lifelong stress susceptibility.

Mapping of flumioxazin tolerance in a snap bean diversity panel leads to the discovery of a master genomic region controlling multiple stress resistance genes.

Effective weed management tools are crucial for maintaining the profitable production of snap bean (Phaseolus vulgaris L.). Preemergence herbicides help the crop to gain a size advantage over the weeds, but the few preemergence herbicides registered in snap bean have poor waterhemp (Amaranthus tuberculatus) control, a major pest in snap bean production. Waterhemp and other difficult-to-control weeds can be managed by flumioxazin, an herbicide that inhibits protoporphyrinogen oxidase (PPO). However, there is limited knowledge about crop tolerance to this herbicide. We aimed to quantify the degree of snap bean tolerance to flumioxazin and explore the underlying mechanisms. We investigated the genetic basis of herbicide tolerance using genome-wide association mapping approach utilizing field-collected data from a snap bean diversity panel, combined with gene expression data of cultivars with contrasting response. The response to a preemergence application of flumioxazin was measured by assessing plant population density and shoot biomass variables. Snap bean tolerance to flumioxazin is associated with a single genomic location in chromosome 02. Tolerance is influenced by several factors, including those that are indirectly affected by seed size/weight and those that directly impact the herbicide's metabolism and protect the cell from reactive oxygen species-induced damage. Transcriptional profiling and co-expression network analysis identified biological pathways likely involved in flumioxazin tolerance, including oxidoreductase processes and programmed cell death. Transcriptional regulation of genes involved in those processes is possibly orchestrated by a transcription factor located in the region identified in the GWAS analysis. Several entries belonging to the Romano class, including Bush Romano 350, Roma II, and Romano Purpiat presented high levels of tolerance in this study. The alleles identified in the diversity panel that condition snap bean tolerance to flumioxazin shed light on a novel mechanism of herbicide tolerance and can be used in crop improvement.

Tracking-seq reveals the heterogeneity of off-target effects in CRISPR-Cas9-mediated genome editing.

The continued development of novel genome editors calls for a universal method to analyze their off-target effects. Here we describe a versatile method, called Tracking-seq, for in situ identification of off-target effects that is broadly applicable to common genome-editing tools, including Cas9, base editors and prime editors. Through tracking replication protein A (RPA)-bound single-stranded DNA followed by strand-specific library construction, Tracking-seq requires a low cell input and is suitable for in vitro, ex vivo and in vivo genome editing, providing a sensitive and practical genome-wide approach for off-target detection in various scenarios. We show, using the same guide RNA, that Tracking-seq detects heterogeneity in off-target effects between different editor modalities and between different cell types, underscoring the necessity of direct measurement in the original system.

Genome sequence analyses identify novel risk loci for multiple system atrophy.

Multiple system atrophy (MSA) is an adult-onset, sporadic synucleinopathy characterized by parkinsonism, cerebellar ataxia, and dysautonomia. The genetic architecture of MSA is poorly understood, and treatments are limited to supportive measures. Here, we performed a comprehensive analysis of whole genome sequence data from 888 European-ancestry MSA cases and 7,128 controls to systematically investigate the genetic underpinnings of this understudied neurodegenerative disease. We identified four significantly associated risk loci using a genome-wide association study approach. Transcriptome-wide association analyses prioritized USP38-DT, KCTD7, and lnc-KCTD7-2 as novel susceptibility genes for MSA within these loci, and single-nucleus RNA sequence analysis found that the associated variants acted as cis-expression quantitative trait loci for multiple genes across neuronal and glial cell types. In conclusion, this study highlights the role of genetic determinants in the pathogenesis of MSA, and the publicly available data from this study represent a valuable resource for investigating synucleinopathies.

Discovery of runs-of-homozygosity diplotype clusters and their associations with diseases in UK Biobank.

Runs-of-homozygosity (ROH) segments, contiguous homozygous regions in a genome were traditionally linked to families and inbred populations. However, a growing literature suggests that ROHs are ubiquitous in outbred populations. Still, most existing genetic studies of ROH in populations are limited to aggregated ROH content across the genome, which does not offer the resolution for mapping causal loci. This limitation is mainly due to a lack of methods for the efficient identification of shared ROH diplotypes. Here, we present a new method, ROH-DICE (runs-of-homozygous diplotype cluster enumerator), to find large ROH diplotype clusters, sufficiently long ROHs shared by a sufficient number of individuals, in large cohorts. ROH-DICE identified over 1 million ROH diplotypes that span over 100 single nucleotide polymorphisms (SNPs) and are shared by more than 100 UK Biobank participants. Moreover, we found significant associations of clustered ROH diplotypes across the genome with various self-reported diseases, with the strongest associations found between the extended human leukocyte antigen (HLA) region and autoimmune disorders. We found an association between a diplotype covering the homeostatic iron regulator (HFE) gene and hemochromatosis, even though the well-known causal SNP was not directly genotyped or imputed. Using a genome-wide scan, we identified a putative association between carriers of an ROH diplotype in chromosome 4 and an increase in mortality among COVID-19 patients (p-value = 1.82 × 10-11). In summary, our ROH-DICE method, by calling out large ROH diplotypes in a large outbred population, enables further population genetics into the demographic history of large populations. More importantly, our method enables a new genome-wide mapping approach for finding disease-causing loci with multi-marker recessive effects at a population scale.

Genetic Foundation of Male Spur Length and Its Correlation with Female Egg Production in Chickens.

Spurs, which mainly appear in roosters, are protrusions near the tarsometatarsus on both sides of the calves of chickens, and are connected to the tarsometatarsus by a bony core. As a male-biased morphological characteristic, the diameter and length of spurs vary significantly between different individuals, mainly related to genetics and age. As a specific behavior of hens, egg-laying also varies greatly between individuals in terms of traits such as age at first egg (AFE), egg weight (EW), and so on. At present, there are few studies on chicken spurs. In this study, we investigated the inheritance pattern of the spur trait in roosters with different phenotypes and the correlations between spur length, body weight at 18 weeks of age (BW18), shank length at 18 weeks of age (SL18), and the egg-laying trait in hens (both hens and roosters were from the same population and were grouped according to their family). These traits related to egg production included AFE, body weight at first egg (BWA), and first egg weight (FEW). We estimated genetic parameters based on pedigree and phenotype data, and used variance analysis to calculate broad-sense heritability for correcting the parameter estimation results. The results showed that the heritability of male left and right spurs ranged from 0.6 to 0.7. There were significant positive correlations between left and right spur length, BW18, SL18, and BWA, as well as between left and right spur length and AFE. We selected 35 males with the longest spurs and 35 males with the shortest spurs in the population, and pooled them into two sets to obtain the pooled genome sequencing data. After genome-wide association and genome divergency analysis by FST, allele frequency differences (AFDs), and XPEHH methods, we identified 7 overlapping genes (CENPE, FAT1, FAM149A, MANBA, NFKB1, SORBS2, UBE2D3) and 14 peak genes (SAMD12, TSPAN5, ENSGALG00000050071, ENSGALG00000053133, ENSGALG00000050348, CNTN5, TRPC6, ENSGALG00000047655,TMSB4X, LIX1, CKB, NEBL, PRTFDC1, MLLT10) related to left and right spur length through genome-wide selection signature analysis and a genome-wide association approach. Our results identified candidate genes associated with chicken spurs, which helps to understand the genetic mechanism of this trait and carry out subsequent research around it.

Ribosomal DNA copy number is associated with body mass in humans and other mammals

Body mass results from a complex interplay between genetics and environment. Previous studies of the genetic contribution to body mass have excluded repetitive regions due to the technical limitations of platforms used for population scale studies. Here we apply genome-wide approaches, identifying an association between adult body mass and the copy number (CN) of 47S-ribosomal DNA (rDNA). rDNA codes for the 18 S, 5.8 S and 28 S ribosomal RNA (rRNA) components of the ribosome. In mammals, there are hundreds of copies of these genes. Inter-individual variation in the rDNA CN has not previously been associated with a mammalian phenotype. Here, we show that rDNA CN variation associates with post-pubertal growth rate in rats and body mass index in adult humans. rDNA CN is not associated with rRNA transcription rates in adult tissues, suggesting the mechanistic link occurs earlier in development. This aligns with the observation that the association emerges by early adulthood.

A high-density genome-wide approach reveals novel genetic markers linked to small ruminant lentivirus susceptibility in sheep.

Visna/Maedi virus (VMV) is lentiviral disease of sheep responsible for severe production losses. Multiple genomic regions associated with infection were reported indicating genetic complexity. In this study, a combined genome-wide approach using a high-density SNP array has been performed, comparing VMV-infected (n = 78) and non-infected (n = 66) individuals of the Valle del Belice breed. The serological tests showed a seroprevalence of 26%. The comparison among results from different approaches (GWAS, Fisher's exact test and the FST analysis) revealed two association signals: on OAR03 close to the GRIN2B gene and on OAR05 close to the TMEM232 gene. To the best of our knowledge, there has been no previous association between these genes and lentiviral infection in any species. The GRIN2B gene plays a role in pain response, synaptic transmission, and receptor clustering, while TMEM232 is involved in the development of immune-related disorders. The results highlighted new aspects of the genetic complexity related to the resistance/susceptibility to VMV in sheep, confirming that studies on different breeds can lead to different results. The ideal approach for validation of the markers identified in our study is to use samples from a population independent from the discovery population with the same phenotype used in the discovery stage.

Unveiling genetic signatures associated with resilience to neonatal diarrhea in lambs through two GWAS approaches

Neonatal diarrhea presents a significant global challenge due to its multifactorial etiology, resulting in high morbidity and mortality rates, and substantial economic losses. While molecular-level studies on genetic resilience/susceptibility to neonatal diarrhea in farm animals are scarce, prior observations indicate promising research directions. Thus, the present study utilizes two genome-wide association approaches, pKWmEB and MLM, to explore potential links between genetic variations in innate immunity and neonatal diarrhea in Karacabey Merino lambs. Analyzing 707 lambs, including 180 cases and 527 controls, revealed an overall prevalence rate of 25.5%. The pKWmEB analysis identified 13 significant SNPs exceeding the threshold of ≥ LOD 3. Moreover, MLM detected one SNP (s61781.1) in the SLC22A8 gene (p-value, 1.85eE-7), which was co-detected by both methods. A McNemar’s test was conducted as the final assessment to identify whether there are any major effective markers among the detected SNPs. Results indicate that four markers—oar3_OAR1_122352257, OAR17_77709936.1, oar3_OAR18_17278638, and s61781.1—have a substantial impact on neonatal diarrhea prevalence (odds ratio: 2.03 to 3.10; statistical power: 0.88 to 0.99). Therefore, we propose the annotated genes harboring three of the associated markers, TIAM1, YDJC, and SLC22A8, as candidate major genes for selective breeding against neonatal diarrhea.

Whole-exome sequencing for genetic diagnosis of idiopathic liver injury in children.

Genome-wide approaches, such as whole-exome sequencing (WES), are widely used to decipher the genetic mechanisms underlying inter-individual variability in disease susceptibility. We aimed to dissect inborn monogenic determinants of idiopathic liver injury in otherwise healthy children. We thus performed WES for 20 patients presented with paediatric-onset recurrent elevated transaminases (rELT) or acute liver failure (ALF) of unknown aetiology. A stringent variant screening was undertaken on a manually-curated panel of 380 genes predisposing to inherited human diseases with hepatobiliary involvement in the OMIM database. We identified rare nonsynonymous variants in nine genes in six patients (five rELT and one ALF). We next performed a case-level evaluation to assess the causal concordance between the gene mutated and clinical symptoms of the affected patient. A genetic diagnosis was confirmed in four rELT patients (40%), among whom two carried novel mutations in ACOX2 or PYGL, and two had previously-reported morbid variants in ABCB4 or PHKA2. We also detected rare variants with uncertain clinical significance in CDAN1, JAG1, PCK2, SLC27A5 or VPS33B in rELT or ALF patients. In conclusion, implementation of WES improves diagnostic yield and enables precision management in paediatric cases of liver injury with unknown aetiology, in particular recurrent hypertransaminasemia.

Classification of HER2 status across multiple cancers using epigenomic profiles from a novel liquid biopsy assay.

3047 Background: The efficacy of trastuzumab deruxtecan in multiple cancers led to a recent FDA granting of priority review for the treatment of adults with HER2-positive solid tumors. This underscores the critical need for reliable assessment of HER2 expression across cancers. The current standard using IHC/ISH for HER2 scoring is fraught with challenges including scoring discordance and tumor heterogeneity. Here we describe the ability to classify HER2 status from 1 ml of plasma, using epigenomic signatures from a novel multi-analyte liquid biopsy (LBx) platform, offering a minimally invasive approach for patient (pt) selection across cancers. Methods: We selected 179 samples from 172 pts with advanced breast (BC), gastro-esophageal (GEA) and ovarian (OV) cancers who had associated HER2 status scored from tissue-based IHC/ISH according to ASCO/CAP guidelines (table). Samples were taken at baseline or at progression and those with detectable cell free DNA, as assessed by iChorCNA, were profiled for genome-wide epigenomic signals across histone modifications associated with active enhancers, promoters and DNA methylation. A regularized regression model developed to classify HER2 status in BC cell lines and was refined and validated for HER2 status prediction for each pt cohort (HER2+ = 3+, 2+/ISH+; HER2- = 2+/ISH-, 1+, 0). The BC cell-line derived HER2 classifier was refined to incorporate GEA and OV cancer specific features before being applied to those samples. Performance was assessed via AUC in a leave-one-out cross-validation schema. We also evaluated HER2 status at progression in a subset of pts with benchmarked HER2 IHC to assess dynamic changes in receptor status by LBx. Results: The epigenomic HER2 classifier was applied to all samples with ctDNA detectable by ichorCNA (90 of 179; 50%). HER2 classification of BC pts by epigenomic liquid biopsy was concordant with standard tissue-based IHC for 64/72 (89%) BC samples (AUC 0.9, table). HER2 classifier predictions for all longitudinally collected samples were concordant with IHC-based HER2 status, including the one patient whose status switched from HER2+ to HER2- at progression. Accurate classification of 11/14 (79%, GEA) and 4/4 (100%, OV) pts was achieved using the indication-refined HER2 classifier. Conclusions: We demonstrate proof of concept for a HER2 classification approach using comprehensive epigenomic signals from 1 ml of plasma that could be applied across multiple cancers. With further development, our genome-wide profiling approach could alleviate clinical constraints associated with multiple tissue-based HER2 scoring assays (IHC/ISH) and enable longitudinal monitoring of HER2 status on therapy. [Table: see text]

Genome-wide CRISPR screening identifies the pivotal role of ANKRD42 in colorectal cancer metastasis through EMT regulation.

Colorectal cancer (CRC), a pervasive and lethal malignancy of gastrointestinal cancer, imposes significant challenges due to the occurrence of distant metastasis in advanced stages. Understanding the intricate regulatory mechanisms driving CRC distant metastasis is of paramount importance. CRISPR-Cas9 screening has emerged as a powerful tool for investigating tumor initiation and progression. However, its application in studying CRC distant metastasis remains largely unexplored. To establish a model that faithfully recapitulates CRC liver metastasis in patients, we developed an in vivo genome-wide CRISPR-Cas9 screening approach using a spleen-injected liver metastasis mouse model. Through comprehensive screening of a whole-genome sgRNA library, we identified ANKRD42 as a pivotal regulatory gene facilitating CRC liver metastasis. Analysis of the TCGA database and our clinical cohorts unveiled heightened ANKRD42 expression in metastases. At the cellular level, the attenuation of ANKRD42 impaired the migration and invasion processes of tumor cells. In vivo experiments further validated these observations, highlighting the diminished liver metastatic capacity of tumor cells upon ANKRD42 knockdown. To unravel the specific mechanisms by which ANKRD42 regulates CRC distant metastasis, we leveraged patient-derived organoid (PDO) models. Depleting ANKRD42 in PDOs sourced from liver metastases precipitated the downregulation of pivotal genes linked to epithelial-mesenchymal transition (EMT), including CDH2 and SNAI2, thereby effectively suppressing tumor metastasis. This study not only establishes a conceptual framework but also identifies potential therapeutic avenues for advanced-stage distant metastasis in CRC patients.

RadioGWAS links radiome to genome to discover driver genes with somatic mutations for heterogeneous tumor image phenotype in pancreatic cancer

Addressing the significant level of variability exhibited by pancreatic cancer necessitates the adoption of a systems biology approach that integrates molecular data, biological properties of the tumors, medical images, and clinical features of the patients. In this study, a comprehensive multi-omics methodology was employed to examine a distinctive collection of patient dataset containing rapid autopsy tumor and normal tissue samples as well as longitudinal imaging with a focus on pancreatic cancer. By performing a whole exome sequencing analysis on tumor and normal tissues to identify somatic gene variants and a radiomic feature analysis to tumor CT images, the genome-wide association approach established a connection between pancreatic cancer driver genes and relevant radiomic features, enabling a thorough and quantitative assessment of the heterogeneity of pancreatic tumors. The significant association between sets of genes and radiomic features revealed the involvement of genes in shaping tumor morphological heterogeneity. Some results of the association established a connection between the molecular level mechanism and their outcomes at the level of tumor structural heterogeneity. Because tumor structure and tumor structural heterogeneity are related to the patients’ overall survival, patients who had pancreatic cancer driver gene mutations with an association to a certain radiomic feature have been observed to experience worse survival rates than cases without these somatic mutations. Furthermore, the association analysis has revealed potential gene mutations and radiomic feature candidates that warrant further investigation in future research endeavors.