Alleles In Genome Research Articles

Adaptive evolution of A-to-I auto-editing site in Adar of eusocial insects

BackgroundAdenosine-to-inosine (A-to-I) RNA editing is a co-/post-transcriptional modification introducing A-to-G variations in RNAs. There is extensive discussion on whether the flexibility of RNA editing exerts a proteomic diversification role, or it just acts like hardwired mutations to correct the genomic allele. Eusocial insects evolved the ability to generate phenotypically differentiated individuals with the same genome, indicating the involvement of epigenetic/transcriptomic regulation.MethodsWe obtained the genomes of 104 Hymenoptera insects and the transcriptomes of representative species. Comparative genomic analysis was performed to parse the evolutionary trajectory of a regulatory Ile > Met auto-recoding site in Adar gene.ResultsAt genome level, the pre-editing Ile codon is conserved across a node containing all eusocial hymenopterans. At RNA level, the editing events are confirmed in representative species and shows considerable condition-specificity. Compared to random expectation, the editable Ile codon avoids genomic substitutions to Met or to uneditable Ile codons, but does not avoid mutations to other unrelated amino acids.ConclusionsThe flexibility of Adar auto-recoding site in Hymenoptera is selectively maintained, supporting the flexible RNA editing hypothesis. We proposed a new angle to view the adaptation of RNA editing, providing another layer to explain the great phenotypical plasticity of eusocial insects.

Genomic co-localization of variation affecting agronomic and human gut microbiome traits in a meta-analysis of diverse sorghum.

Substantial functional metabolic diversity exists within species of cultivated grain crops that directly or indirectly provide more than half of all calories consumed by humans around the globe. While such diversity is the molecular currency used for improving agronomic traits, diversity is poorly characterized for its effects on human nutrition and utilization by gut microbes. Moreover, we know little about agronomic traits' potential tradeoffs and pleiotropic effects on human nutritional traits. Here, we applied a quantitative genetics approach using a meta-analysis and parallel genome-wide association studies of Sorghum bicolor traits describing changes in the composition and function of human gut microbe communities, and any of 200 sorghum seed and agronomic traits across a diverse sorghum population to identify associated genetic variants. A total of 15 multiple-effect loci (MEL) were initially found where different alleles in the sorghum genome produced changes in seed that affected the abundance of multiple bacterial taxa across 2 human microbiomes in automated in vitro fermentations. Next, parallel genome-wide studies conducted for seed, biochemical, and agronomic traits in the same population identified significant associations within the boundaries of 13/15 MEL for microbiome traits. In several instances, the colocalization of variation affecting gut microbiome and agronomic traits provided hypotheses for causal mechanisms through which variation could affect both agronomic traits and human gut microbes. This work demonstrates that genetic factors affecting agronomic traits in sorghum seed can also drive significant effects on human gut microbes, particularly bacterial taxa considered beneficial. Understanding these pleiotropic relationships will inform future strategies for crop improvement toward yield, sustainability, and human health.

Machine learning-based classification reveals distinct clusters of non-coding genomic allelic variations associated with Erm-mediated antibiotic resistance.

Antibiotic resistance (AR) poses a global health threat as the efficacy of available antibiotics has rapidly eroded due to the widespread transmission of AR genes. Using Erm-dependent MLS resistance as a model, this study highlights the significance of non-coding genomic allelic variations. Through a comprehensive analysis of upstream regulatory elements within the erm family, we elucidated the evolutionary emergence and development of AR mechanisms. Leveraging population-wide machine learning (ML)-based genomic analysis, we transformed substantial non-random allelic variations into discernible clusters of elements, enabling precise prediction of MLS phenotypes from non-coding regions. These findings offer deeper insight into AR evolution and demonstrate the potential of harnessing non-coding genomic allele data for accurately predicting AR phenotypes.

The cacao gene atlas: a transcriptome developmental atlas reveals highly tissue-specific and dynamically-regulated gene networks in Theobroma cacao L

BackgroundTheobroma cacao, the cocoa tree, is a tropical crop grown for its highly valuable cocoa solids and fat which are the basis of a 200-billion-dollar annual chocolate industry. However, the long generation time and difficulties associated with breeding a tropical tree crop have limited the progress of breeders to develop high-yielding disease-resistant varieties. Development of marker-assisted breeding methods for cacao requires discovery of genomic regions and specific alleles of genes encoding important traits of interest. To accelerate gene discovery, we developed a gene atlas composed of a large dataset of replicated transcriptomes with the long-term goal of progressing breeding towards developing high-yielding elite varieties of cacao.ResultsWe describe the creation of the Cacao Transcriptome Atlas, its global characterization and define sets of genes co-regulated in highly organ- and temporally-specific manners. RNAs were extracted and transcriptomes sequenced from 123 different tissues and stages of development representing major organs and developmental stages of the cacao lifecycle. In addition, several experimental treatments and time courses were performed to measure gene expression in tissues responding to biotic and abiotic stressors. Samples were collected in replicates (3–5) to enable statistical analysis of gene expression levels for a total of 390 transcriptomes. To promote wide use of these data, all raw sequencing data, expression read mapping matrices, scripts, and other information used to create the resource are freely available online. We verified our atlas by analyzing the expression of genes with known functions and expression patterns in Arabidopsis (ACT7, LEA19, AGL16, TIP13, LHY, MYB2) and found their expression profiles to be generally similar between both species. We also successfully identified tissue-specific genes at two thresholds in many tissue types represented and a set of genes highly conserved across all tissues.ConclusionThe Cacao Gene Atlas consists of a gene expression browser with graphical user interface and open access to raw sequencing data files as well as the unnormalized and CPM normalized read count data mapped to several cacao genomes. The gene atlas is a publicly available resource to allow rapid mining of cacao gene expression profiles. We hope this resource will be used to help accelerate the discovery of important genes for key cacao traits such as disease resistance and contribute to the breeding of elite varieties to help farmers increase yields.

Observed germline BRCA1/2 and TP53 tumor genomic profiling allele frequencies.

e22533 Background: Tumor genomic profiling (TGP) is increasingly used to assist with cancer treatment planning. While these tests are not intended to find germline variants, they are a possible incidental outcome. Currently there are no clear guidelines on when to suspect if a somatic variant identified on TGP may need confirmatory germline genetic testing. Variant allele frequency (VAF) is one tool commonly used to assess if a variant may be germline. While some studies have demonstrated germline VAF ranges on TGP, the data is limited. The purpose of this study was to review VAF ranges of BRCA1/2 and TP53 variants at a large reference laboratory to better understand VAF ranges of confirmed germline variants. Methods: A retrospective review of test results from 184 patients who underwent TGP and hereditary cancer germline testing at a large reference laboratory was performed. Patients were included in the study if they had both TGP and germline testing performed, regardless of the order of testing. Results were analyzed to determine variants identified and VAF. Analysis of BRCA1/2 and TP53 variants included patients who had variants identified on TGP and had germline testing that would have confirmed the variants if present in the germline. The VAFs were obtained, and simple statistics were used to describe the findings. Results: Twenty-one patients were found to have 23 BRCA1/2 variants identified on TGP (2 patients had 2 BRCA1/2 TGP variants). Of these TGP variants, 16 were confirmed on germline testing (69.6%). The TGP VAF for these patients ranged from 38.2% to 85.1%. The mean was 61.0% and the median was 61.8%. The remaining 7 TGP variants that were not confirmed on germline testing ranged from 7.2% to 63.4%. The mean was 31.1% and the median was 28%. Thirty-eight patients were found to have 41 TP53 variants identified on TGP (3 patients had 2 TP53 TGP variants). Of these patients, 1 had a TP53 variant confirmed on germline testing with a VAF of 44.7%. The 40 remaining variants in this dataset were not confirmed on germline testing. The VAF for these 40 variants ranged from 1.3% to 95.3%. The mean was 38.4% and the median was 34.5%. Conclusions: In this study, 69.6% of BRCA1/2 TGP variants were confirmed on germline testing compared to 2.4% of TP53 TGP variants. Confirmed germline variants identified on TGP had VAFs that were > 38%. It’s important to note that TP53 variants may have very high VAFs (up to 95.3% in this study) and not be.

Polyploidisation pleiotropically buffers ageing in hepatocytes

Polyploidy in hepatocytes has been proposed as a genetic mechanism to buffer against transcriptional dysregulation. Here, we aim to demonstrate the role of polyploidy in modulating gene regulatory networks in hepatocytes during ageing. We performed single-nucleus RNA sequencing in hepatocyte nuclei of different ploidy levels isolated from young and old wild-type mice. Changes in the gene expression and regulatory network were compared to three independent strains that were haploinsufficient for HNF4A, CEBPA or CTCF, representing non-deleterious perturbations. Phenotypic characteristics of the liver section were additionally evaluated histologically, whereas the genomic allele composition of hepatocytes was analysed by BaseScope. We observed that ageing in wild-type mice results in nuclei polyploidy and a marked increase in steatosis. Haploinsufficiency of liver-specific master regulators (HFN4A or CEBPA) results in the enrichment of hepatocytes with tetraploid nuclei at a young age, affecting the genomic regulatory network, and dramatically suppressing ageing-related steatosis tissue wide. Notably, these phenotypes are not the result of subtle disruption to liver-specific transcriptional networks, since haploinsufficiency in the CTCF insulator protein resulted in the same phenotype. Further quantification of genotypes of tetraploid hepatocytes in young and old HFN4A-haploinsufficient mice revealed that during ageing, tetraploid hepatocytes lead to the selection of wild-type alleles, restoring non-deleterious genetic perturbations. Our results suggest a model whereby polyploidisation leads to fundamentally different cell states. Polyploid conversion enables pleiotropic buffering against age-related decline via non-random allelic segregation to restore a wild-type genome. The functional role of hepatocyte polyploidisation during ageing is poorly understood. Using single-nucleus RNA sequencing and BaseScope approaches, we have studied ploidy dynamics during ageing in murine livers with non-deleterious genetic perturbations. We have identified that hepatocytes present different cellular states and the ability to buffer ageing-associated dysfunctions. Tetraploid nuclei exhibit robust transcriptional networks and are better adapted to genomically overcome perturbations. Novel therapeutic interventions aimed at attenuating age-related changes in tissue function could be exploited by manipulation of ploidy dynamics during chronic liver conditions.

A Hypermutable Region in the DISP2 Gene Links to Natural Selection and Late-Onset Neurocognitive Disorders in Humans.

(CCG) short tandem repeats (STRs) are predominantly enriched in genic regions, mutation hotspots for C to T truncating substitutions, and involved in various neurological and neurodevelopmental disorders. However, intact blocks of this class of STRs are widely overlooked with respect to their link with natural selection. The human neuron-specific gene, DISP2 (dispatched RND transporter family member 2), contains a (CCG) repeat in its 5' untranslated region. Here, we sequenced this STR in a sample of 448 Iranian individuals, consisting of late-onset neurocognitive disorder (NCD) (N = 203) and controls (N = 245). We found that the region spanning the (CCG) repeat was highly mutated, resulting in several flanking (CCG) residues. However, an 8-repeat of the (CCG) repeat was predominantly abundant (frequency = 0.92) across the two groups. While the overall distribution of genotypes was not different between the two groups (p > 0.05), we detected four genotypes in the NCD group only (2% of the NCD genotypes, Mid-p = 0.02), consisting of extreme short alleles, 5- and 6-repeats, that were not detected in the control group. The patients harboring those genotypes received the diagnoses of probable Alzheimer's disease and vascular dementia. We also found six genotypes in the control group only (2.5% of the control genotypes, Mid-p = 0.01) that consisted of the 8-repeat and extreme long alleles, 9- and 10-repeats, of which the 10-repeat was not detected in the NCD group. The (CCG) repeat specifically expanded in primates. In conclusion, we report an indication of natural selection at a novel hypermutable region in the human genome and divergent alleles and genotypes in late-onset NhCDs and controls. These findings reinforce the hypothesis that a collection of rare alleles and genotypes in a number of genes may unambiguously contribute to the cognition impairment component of late-onset NCDs.

Identification of germline population variants misclassified as cancer-associated somatic variants.

Databases used for clinical interpretation in oncology rely on genetic data derived primarily from patients of European ancestry, leading to biases in cancer genetics research and clinical practice. One practical issue that arises in this context is the potential misclassification of multi-ancestral population variants as tumor-associated because they are not represented in reference genomes against which tumor sequencing data is aligned. To systematically find misclassified variants, we compared somatic variants in census genes from the Catalogue of Somatic Mutations in Cancer (COSMIC) V99 with multi-ancestral population variants from the Genome Aggregation Databases' Linkage Disequilibrium (GnomAD). By comparing genomic coordinates, reference, and alternate alleles, we could identify misclassified variants in genes associated with cancer. We found 192 of 208 genes in COSMIC's cancer-associated census genes (92.31%) to be associated with variant misclassifications. Among the 1,906,732 variants in COSMIC, 6,957 variants (0.36%) aligned with normal population variants in GnomAD, concerning for misclassification. The African / African American ancestral population included the greatest number of misclassified variants and also had the greatest number of unique misclassified variants. The direct, systematic comparison of variants from COSMIC for co-occurrence in GnomAD supports a more accurate interpretation of tumor sequencing data and reduces bias related to genomic ancestry.

Genome wide association study reveals key genic regions controlling oil palm vegetative characters in wild population of Cameroon

Cameroon’s oil palm germplasm is contributing a lot towards the improvement of genetic variation in the base oil palm population used in breeding programs around the world. There is a new oil palm germplasm collection at Dibamba constituted of 169 wild accessions but whose genomic study was not implemented before now. Genome-wide-association study is considered as a powerful tool to identify genomic specific allele variants controlling the expression of important agronomic traits in crops. This approach could potentially accelerate varietal improvement in plant breeding programs. The objective of this study was to identify genic regions across the investigated oil palm genomes which are conferred to control natural variation for oil palm vegetative traits such as leaflet width and leaf area. The allele and genotype frequencies were significantly detected to be in Hardy-Weinberg equilibrium (P<0.01) for these vegetative traits. Across the genome, important numbers of single nucleotide polymorphisms were associated with the oil palm leaflet width and leaf area characters. This suggests that the considered regions may contain genes controlling the phenotype variation expression of the trait of interest and should be useful under positive selection in subsequent breeding of the oil palm.

Identification of eight genomic protective alleles for mitochondrial diabetes by Kinship-graph convolutional network.

Nearly 85% of maternally inherited diabetes and deafness (MIDD) are caused by the m.3243A>G mutation in the mitochondrial DNA. However, the clinical phenotypes of MIDD may also be influenced by the nuclear genome, this study aimed to investigate nuclear genome variants that influence clinical phenotypes associated with m.3243A>G mutation in MIDD based on whole-genome sequencing of the patients belonging to pedigrees. We analyzed a whole-genome sequencing (WGS) dataset from blood samples of 38 MIDD patients with the m.3243A > G mutation belonging to 10 pedigrees, by developing a Kinship-graph convolutional network approach, called Ki-GCN, integrated with the conventional genome-wide association study (GWAS) methods. We identified eight protective alleles in the nuclear genome that have protective effects against the onset of MIDD, related deafness, and also type 2 diabetes. Based on these eight protective alleles, we constructed an effective logistic regression model to predict the early or late onset of MIDD patients. There are protective alleles in the nuclear genome that are associated with the onset-age of MIDD patients and might also provide protective effects on the deafness derived from MIDD patients.

Genomic prediction and allele mining of agronomic and morphological traits in pea (Pisum sativum) germplasm collections.

Well-performing genomic prediction (GP) models for polygenic traits and molecular marker sets for oligogenic traits could be useful for identifying promising genetic resources in germplasm collections, setting core collections, and establishing molecular variety distinction. This study aimed at (i) defining GP models and key marker sets for predicting 15 agronomic or morphological traits in germplasm collections, (ii) verifying the GP model usefulness also for selection in breeding programs, (iii) investigating the consistency between molecular and phenotypic diversity patterns, and (iv) identifying genomic regions associated with to the target traits. The study was based on phenotyping data and over 41,000 genotyping-by-sequencing-generated SNP markers of 220 landraces or old cultivars belonging to a world germplasm collection and 11 modern cultivars. Non-metric multi-dimensional scaling (NMDS) and an analysis of population genetic structure indicated a high level of genetic differentiation of material from Western Asia, a major West-East diversity gradient, and quite limited genetic diversity of the improved germplasm. Mantel's test revealed a low correlation (r = 0.12) between phenotypic and molecular diversity, which increased (r = 0.45) when considering only the molecular diversity relative to significant SNPs from genome-wide association analyses. These analyses identified, inter alia, several areas of chromosome 6 involved in a largely pleiotropic control of vegetative or reproductive organ pigmentation. We found various significant SNPs for grain and straw yield under severe drought and onset of flowering, and one SNP on chromosome 5 for grain protein content. GP models displayed moderately high predictive ability (0.43 to 0.61) for protein content, grain and straw yield, and onset of flowering, and high predictive ability (0.76) for individual seed weight, based on intra-population, intra-environment cross-validations. The inter-population, inter-environment assessment of the models trained on the germplasm collection for breeding material of three recombinant inbred line (RIL) populations, which was challenged by much narrower diversity of the material, over eight-fold less available markers and quite different test environments, led to an overall loss of predictive ability of about 40% for seed weight, 50% for protein content and straw yield, and 60% for onset of flowering, and no prediction for grain yield. Within-RIL population predictive ability differed among populations.

Opposing, spatially-determined epigenetic forces impose restrictions on stochastic olfactory receptor choice.

Olfactory receptor (OR) choice represents an example of genetically hardwired stochasticity, where every olfactory neuron expresses one out of ~2000 OR alleles in the mouse genome in a probabilistic, yet stereotypic fashion. Here, we propose that topographic restrictions in OR expression are established in neuronal progenitors by two opposing forces: polygenic transcription and genomic silencing, both of which are influenced by dorsoventral gradients of transcription factors NFIA, B, and X. Polygenic transcription of OR genes may define spatially constrained OR repertoires, among which one OR allele is selected for singular expression later in development. Heterochromatin assembly and genomic compartmentalization of OR alleles also vary across the axes of the olfactory epithelium and may preferentially eliminate ectopically expressed ORs with more dorsal expression destinations from this 'privileged' repertoire. Our experiments identify early transcription as a potential 'epigenetic' contributor to future developmental patterning and reveal how two spatially responsive probabilistic processes may act in concert to establish deterministic, precise, and reproducible territories of stochastic gene expression.

Opposing, spatially-determined epigenetic forces impose restrictions on stochastic olfactory receptor choice

Olfactory receptor (OR) choice represents an example of genetically hardwired stochasticity, where every olfactory neuron expresses one out of ~2000 OR alleles in the mouse genome in a probabilistic, yet stereotypic fashion. Here, we propose that topographic restrictions in OR expression are established in neuronal progenitors by two opposing forces: polygenic transcription and genomic silencing, both of which are influenced by dorsoventral gradients of transcription factors NFIA, B, and X. Polygenic transcription of OR genes may define spatially constrained OR repertoires, among which one OR allele is selected for singular expression later in development. Heterochromatin assembly and genomic compartmentalization of OR alleles also vary across the axes of the olfactory epithelium and may preferentially eliminate ectopically expressed ORs with more dorsal expression destinations from this ‘privileged’ repertoire. Our experiments identify early transcription as a potential ‘epigenetic’ contributor to future developmental patterning and reveal how two spatially responsive probabilistic processes may act in concert to establish deterministic, precise, and reproducible territories of stochastic gene expression.

A genomic mutational constraint map using variation in 76,156 human genomes.

The depletion of disruptive variation caused by purifying natural selection (constraint) has been widely used to investigate protein-coding genes underlying human disorders1-4, but attempts to assess constraint for non-protein-coding regions have proved more difficult. Here we aggregate, process and release a dataset of 76,156 human genomes from the Genome Aggregation Database (gnomAD)-the largest public open-access human genome allele frequency reference dataset-and use it to build a genomic constraint map for the whole genome (genomic non-coding constraint of haploinsufficient variation (Gnocchi)). We present a refined mutational model that incorporates local sequence context and regional genomic features to detect depletions of variation. As expected, the average constraint for protein-coding sequences is stronger than thatfor non-coding regions. Within the non-coding genome, constrained regions are enriched for known regulatory elements and variants that are implicated in complex human diseases and traits, facilitating the triangulation of biological annotation, disease association and natural selection to non-coding DNA analysis. More constrained regulatory elements tend to regulate more constrained protein-coding genes, which in turn suggests that non-coding constraint can aid the identification of constrained genes that are as yet unrecognized by current gene constraint metrics. We demonstrate that this genome-wide constraint map improves the identification and interpretation of functional human genetic variation.

Feasibility for non-invasive prenatal fetal blood group and platelet genotyping by massively parallel sequencing: A single test system for multiple atypical red cell, platelet and quality control markers.

Non-invasive prenatal tests (NIPT) to predict fetal red cell or platelet antigen status for alloimmunised women are provided for select antigens. This study reports on massively parallel sequencing (MPS) using a red cell and platelet probe panel targeting multiple nucleotide variants, plus individual identification single nucleotide polymorphisms (IISNPs). Maternal blood samples were provided from 33 alloimmunised cases, including seven with two red cell antibodies. Cell-free and genomic DNA was sequenced using targeted MPS and bioinformatically analysed using low-frequency variant detection. The resulting maternal genomic DNA allele frequency was subtracted from the cell-free DNA counterpart. Outcomes were matched against validated phenotyping/genotyping methods, where available. A 2.5% subtractive allele frequency threshold was set after comparing MPS predictions for K, RhC/c, RhE/e and Fya /Fyb against expected outcomes. This threshold was used for subsequent predictions, including HPA-15a, Jka /Jkb , Kpa /Kpb and Lua . MPS outcomes were 97.2% concordant with validated methods; one RhC case was discordantly negative and lacked IISNPs. IISNPs were informative for 30/33 cases as controls. NIPT MPS is feasible for fetal blood group genotyping and covers multiple blood groups and control targets in a single test. Noting caution for the Rh system, this has the potential to provide a personalised service for alloimmunised women.

Chemical Expansion of the Methyltransferase Reaction: Tools for DNA Labeling and Epigenome Analysis.

DNA is the genetic matter of life composed of four major nucleotides which can be further furnished with biologically important covalent modifications. Among the variety of enzymes involved in DNA metabolism, AdoMet-dependent methyltransferases (MTases) combine the recognition of specific sequences and covalent methylation of a target nucleotide. The naturally transferred methyl groups play important roles in biological signaling, but they are poor physical reporters and largely resistant to chemical derivatization. Therefore, an obvious strategy to unlock the practical utility of the methyltransferase reactions is to enable the transfer of "prederivatized" (extended) versions of the methyl group.However, previous enzymatic studies of extended AdoMet analogs indicated that the transalkylation reactions are drastically impaired as the size of the carbon chain increases. In collaborative efforts, we proposed that, akin to enhanced SN2 reactivity of allylic and propargylic systems, addition of a π orbital next to the transferable carbon atom might confer the needed activation of the reaction. Indeed, we found that MTase-catalyzed transalkylations of DNA with cofactors containing a double or a triple C-C bond in the β position occurred in a robust and sequence-specific manner. Altogether, this breakthrough approach named mTAG (methyltransferase-directed transfer of activated groups) has proven instrumental for targeted labeling of DNA and other types of biomolecules (using appropriate MTases) including RNA and proteins.Our further work focused on the propargylic cofactors and their reactions with DNA cytosine-5 MTases, a class of MTases common for both prokaryotes and eukaryotes. Here, we learned that the 4-X-but-2-yn-1-yl (X = polar group) cofactors suffered from a rapid loss of activity in aqueous buffers due to susceptibility of the triple bond to hydration. This problem was remedied by synthetically increasing the separation between X and the triple bond from one to three carbon units (6-X-hex-2-ynyl cofactors). To further optimize the transfer of the bulkier groups, we performed structure-guided engineering of the MTase cofactor pocket. Alanine replacements of two conserved residues conferred substantial improvements of the transalkylation activity with M.HhaI and three other engineered bacterial C5-MTases. Of particular interest were CpG-specific DNA MTases (M.SssI), which proved valuable tools for studies of mammalian methylomes and chemical probing of DNA function.Inspired by the successful repurposing of bacterial enzymes, we turned to more complex mammalian C5-MTases (Dnmt1, Dnmt3A, and Dnmt3B) and asked if they could ultimately lead to mTAG labeling inside mammalian cells. Our efforts to engineer mouse Dnmt1 produced a variant (Dnmt1*) that enabled efficient Dnmt1-directed deposition of 6-azide-hexynyl groups on DNA in vitro. CRISPR-Cas9 editing of the corresponding codons in the genomic Dnmt1 alleles established endogenous expression of Dnmt1* in mouse embryonic stem cells. To circumvent the poor cellular uptake of AdoMet and its analogs, we elaborated their efficient internalization by electroporation, which has finally enabled selective catalysis-dependent azide tagging of natural Dnmt1 targets in live mammalian cells. The deposited chemical groups were then exploited as "click" handles for reading adjoining sequences and precise genomic mapping of the methylation sites. These findings offer unprecedented inroads into studies of DNA methylation in a wide range of eukaryotic model systems.

Characterizing proteomic and transcriptomic features of missense variants in amyotrophic lateral sclerosis genes.

Within recent years, there has been a growing number of genes associated with amyotrophic lateral sclerosis (ALS), resulting in an increasing number of novel variants, particularly missense variants, many of which are of unknown clinical significance. Here, we leverage the sequencing efforts of the ALS Knowledge Portal (3,864 individuals with ALS and 7,839 controls) and Project MinE ALS Sequencing Consortium (4,366 individuals with ALS and 1,832 controls) to perform proteomic and transcriptomic characterization of missense variants in 24 ALS-associated genes. The two sequencing datasets were interrogated for missense variants in the 24 genes, and variants were annotated with genomic database minor allele frequencies, ClinVar pathogenicity classifications, protein sequence features including Uniprot functional site annotations and PhosphoSitePlus post-translational modification (PTM) site annotations, structural features from AlphaFold predicted monomeric 3D structures, and transcriptomic expression levels from Genotype-Tissue Expression (GTEx). We then applied missense variant enrichment and gene-burden testing following binning of variation based on the selected proteomic and transcriptomic features to identify those most relevant to pathogenicity in ALS-associated genes. Using predicted human protein structures from AlphaFold, we determined that missense variants carried by individuals with ALS were significantly enriched in β-sheets and α-helices, as well as in core, buried, or moderately buried regions. At the same time, we identified that hydrophobic amino acid residues, compositionally biased protein regions and protein-protein interaction regions are predominantly enriched in missense variants carried by individuals with ALS. Assessment of expression level based on transcriptomics also revealed enrichment of variants of high and medium expression across all tissues and within the brain. We further explored enriched features of interest using burden analyses and identified individual genes were indeed driving certain enrichment signals. A case study is presented for SOD1 to demonstrate proof of concept of how enriched features may aid in defining variant pathogenicity. Our results present proteomic and transcriptomic features that are important indicators of missense variant pathogenicity in ALS and are distinct from features associated with neurodevelopmental disorders.

Abstract PD17-01: Genomic analysis of circulating tumor DNA (ctDNA) from patients with HR+, HER2-mutant metastatic breast cancer (MBC) enrolled in SUMMIT: mechanisms of acquired resistance to neratinib + fulvestrant + trastuzumab (N+F+T)

Abstract Background: In the SUMMIT trial, original cohorts of patients with HR+, HER2-negative (locally assessed), HER2-mutant MBC received N alone or N+F, with promising clinical response rates but abbreviated duration. Clinical progression coincided with emergence of additional HER2 mutations and/or amplification of the mutant allele [Smyth et al. Cancer Discov 2020;10:198–213]. Addition of T to the combination was postulated to prolong response; the combination of N+F+T in heavily pretreated patients with HR+, HER2-mutant MBC who had received CDK4/6 inhibitors (n=51) yielded a confirmed overall response rate (ORR) of 35.3%, median duration of response (DOR) of 14.3 months, clinical benefit rate (CBR) of 47.1%, and median progression-free survival (PFS) of 8.2 months [Jhaveri et al. J Clin Oncol 2022;40:1028]. Seven of these 51 patients were part of a randomized (1:1:1) cohort who received N+F+T, F+T, or F alone. Patients randomized to F+T or F could crossover to N+F+T upon progression. No patients responded to F or F+T; however, one of four patients who crossed over to N+F+T upon progression on F+T responded to the triplet, as did two of six who crossed over upon progression on F. We undertook ctDNA sequencing in patients who responded to N+F+T upfront and after crossover. Methods: Patients with HR+, HER2-negative MBC with activating HER2 mutation(s) and prior CDK4/6i received N+F+T (oral N 240 mg/d with loperamide prophylaxis, im F 500 mg d1, d15, and d29 of cycle 1 then q4w, iv T 8 mg/kg initially then 6 mg/kg q3w) or F+T or F alone. Efficacy endpoints included investigator-assessed ORR and CBR (RECIST v1.1), DOR, and best overall response. ctDNA was collected at baseline, every third cycle during treatment, and at the end of treatment and analyzed by next-generation sequencing. Samples were analyzed using the TEMPUS xF assay. Results: Sequencing results from ctDNA analysis are pending at the time of this abstract submission. Baseline HER2 mutations and co-alterations will be reported and compared with those found in tissue samples. Genomic spectrum and variant allele frequencies in samples taken at baseline, on treatment, and at the end of treatment from patients who experienced complete or partial response to N+F+T and then progressed (n=25) will be sequenced and mechanisms of acquired resistance will be posited. ctDNA genomic profiles of serial time points from patients randomized to F or F+T before and after crossover to N+F+T (n=10) will also be evaluated. Conclusions: Similarities and differences between the mechanisms of acquired resistance to N+F+T, and those previously reported to be associated with progression on N or N+F, will be discussed. Citation Format: Cynthia Ma, James Waisman, Adam M. Brufsky, Eddy S. Yang, Hans Wildiers, John P. Crown, Sarina A. Piha-Paul, Jennifer M. Suga, José Ángel García-Sáenz, Valentina Gambardella, Angel Guerrero, Salomon Stemmer, Ron Bose, Tonya Novara-Demgen, Daniel DiPrimeo, Lisa D. Eli, Komal Jhaveri. Genomic analysis of circulating tumor DNA (ctDNA) from patients with HR+, HER2-mutant metastatic breast cancer (MBC) enrolled in SUMMIT: mechanisms of acquired resistance to neratinib + fulvestrant + trastuzumab (N+F+T) [abstract]. In: Proceedings of the 2022 San Antonio Breast Cancer Symposium; 2022 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2023;83(5 Suppl):Abstract nr PD17-01.

Single-cell whole-genome sequencing, haplotype analysis in prenatal diagnosis of monogenic diseases.

Monogenic inherited diseases are common causes of congenital disabilities, leading to severe economic and mental burdens on affected families. In our previous study, we demonstrated the validity of cell-based noninvasive prenatal testing (cbNIPT) in prenatal diagnosis by single-cell targeted sequencing. The present research further explored the feasibility of single-cell whole-genome sequencing (WGS) and haplotype analysis of various monogenic diseases with cbNIPT. Four families were recruited: one with inherited deafness, one with hemophilia, one with large vestibular aqueduct syndrome (LVAS), and one with no disease. Circulating trophoblast cells (cTBs) were obtained from maternal blood and analyzed by single-cell 15X WGS. Haplotype analysis showed that CFC178 (deafness family), CFC616 (hemophilia family), and CFC111 (LVAS family) inherited haplotypes from paternal and/or maternal pathogenic loci. Amniotic fluid or fetal villi samples from the deafness and hemophilia families confirmed these results. WGS performed better than targeted sequencing in genome coverage, allele dropout (ADO), and false-positive (FP) ratios. Our findings suggest that cbNIPT by WGS and haplotype analysis have great potential for use in prenatally diagnosing various monogenic diseases.

Genome-wide analysis of NBS-LRR genes revealed contribution of disease resistance from Saccharum spontaneum to modern sugarcane cultivar.

During plant evolution, nucleotide-binding sites (NBS) and leucine-rich repeat (LRR) genes have made significant contributions to plant disease resistance. With many high-quality plant genomes sequenced, identification and comprehensive analyses of NBS-LRR genes at whole genome level are of great importance to understand and utilize them. In this study, we identified the NBS-LRR genes of 23 representative species at whole genome level, and researches on NBS-LRR genes of four monocotyledonous grass species, Saccharum spontaneum, Saccharum officinarum, Sorghum bicolor and Miscanthus sinensis, were focused. We found that whole genome duplication, gene expansion, and allele loss could be factors affecting the number of NBS-LRR genes in the species, and whole genome duplication is likely to be the main cause of the number of NBS-LRR genes in sugarcane. Meanwhile, we also found a progressive trend of positive selection on NBS-LRR genes. These studies further elucidated the evolutionary pattern of NBS-LRR genes in plants. Transcriptome data from multiple sugarcane diseases revealed that more differentially expressed NBS-LRR genes were derived from S. spontaneum than from S. officinarum in modern sugarcane cultivars, and the proportion was significantly higher than the expected. This finding reveals that S. spontaneum has a greater contribution to disease resistance for modern sugarcane cultivars. In addition, we observed allelespecific expression of seven NBS-LRR genes under leaf scald, and 125 NBS-LRR genes responding to multiple diseases were identified. Finally, we built a plant NBS-LRR gene database to facilitate subsequent analysis and use of NBSLRR genes obtained here. In conclusion, this study complemented and completed the research of plant NBS-LRR genes, and discussed how NBS-LRR genes responding to sugarcane diseases, which provided a guide and genetic resources for further research and utilization of NBS-LRR genes.