Allelic Sequences Research Articles

Identification and molecular marker development for peel color gene in melon (Cucumis melo L.)

Peel color is an important appearance quality of melons that greatly affects consumer preferences. In this study, a near-isogenic line NIL-G (dark green peel) was generated from B8 (grey-green peel) and B15 (white peel). The F2 population constructed by crossing NIL-G and B15 was used to study the inheritance pattern of peel color, and bulked-segregant analysis sequencing (BSA-seq) was employed to identify the interval in which the target gene was located. Genetic analysis results showed that the dark green peel trait at maturity is controlled by a dominant gene. BSA-seq and molecular markers were used to localize the candidate gene in a 263.7 kb interval of chromosome 4, which contained the CmAPRR2 gene with known functions. Moreover, allelic sequence analysis revealed four SNP variations of the CmAPRR2 gene in B15, of which SNP.G614331A was located at the junction of the 6th exon and 6th intron. The G-to-A mutation caused alternative splicing of the transcript of CmAPRR2 in B15, generating two transcripts (CmAPRR2-A and CmAPRR2-B) with premature termination codons. Furthermore, the Kompetitive Allele Specific PCR (KASP) marker, APRR2-G/A, was developed based on this SNP and shown to co-segregate with the peel color phenotype in the F2 population. Compared to white-peel B15, the expression level of CmAPRR2 in dark green peel NIL-G was higher at each growth stage. Therefore, CmAPRR2 may be the key gene controlling the fruit color of melons. Overall, this study identified a novel allelic variant of CmAPRR2 that leads to white peel formation in mature melons. The study also provides a theoretical basis for further research on the gene regulatory mechanism of melon peel colors and may promote the future use of molecular marker-assisted selection to modify melon peel colors.

Sequence analysis and identification of a novel HLA-DPB1*02:01:69 allele by third-generation sequencing

To analyze the sequence of a novel HLA-DPB1 allele in an individual. A individual identified from the database of blood donors for matched platelet transfusion at the Blood Center of Zhejiang Province in May 2022 was selected as the study subject. HLA genotype of the individual was determined by next-generation sequencing (NGS) on an Ion Torrent S5 platform. The sequence of the HLA-DPB1 locus was also determined by NGS on an Illumina Miseq platform and third-generation sequencing using Oxford Nanopore MinION. This study was approved by the Blood Center of Zhejiang Province (Ethics No. 2021-001). A novel HLA-DPB1*02 allele was identified in the specimen, for which the closest genotype was HLA-DPB1*02:new,17:01:01G, with the variant located in exon 3. Meanwhile, the NGS also revealed a novel HLA-DPB1*17 allele, with the closest genotype being HLA-DPB1*02:01,17:new. Both the HLA-DPB1*17:01:01:01 and HLA-DPB1*02 alleles were identified by third-generation sequencing. Compared with the HLA-DPB1*02:01:02:01 allele, the novel allele had a G>A variation at position 369 in the exon 3, which however did not result in amino acid change. A novel HLA-DPB1 allele has been identified and validated by both NGS and TGS, which has been named as HLA-DPB1*02:01:69 by the World Health Organization Committee on Nomenclature of Factors of the HLA System.

Conserved allomorphs of MR1 drive the specificity of MR1-restricted TCRs.

Major histocompatibility complex class-1-related protein (MR1), unlike human leukocyte antigen (HLA) class-1, was until recently considered to be monomorphic. MR1 presents metabolites in the context of host responses to bacterial infection. MR1-restricted TCRs specific to tumor cells have been described, raising interest in their potential therapeutic application for cancer treatment. The diversity of MR1-ligand biology has broadened with the observation that single nucleotide variants (SNVs) exist within MR1 and that allelic variants can impact host immunity. The TCR from a MR1-restricted T-cell clone, MC.7.G5, with reported cancer specificity and pan-cancer activity, was cloned and expressed in Jurkat E6.1 TCRαβ- β2M- CD8+ NF-κB:CFP NFAT:eGFP AP-1:mCherry cells or in human donor T cells. Functional activity of 7G5.TCR-T was demonstrated using cytotoxicity assays and by measuring cytokine release after co-culture with cancer cell lines with or without loading of previously described MR1 ligands. MR1 allele sequencing was undertaken after the amplification of the MR1 gene region by PCR. In vivo studies were undertaken at Labcorp Drug Development (Ann Arbor, MI, USA) or Epistem Ltd (Manchester, UK). The TCR cloned from MC.7.G5 retained MR1-restricted functional cytotoxicity as 7G5.TCR-T. However, activity was not pan-cancer, as initially reported with the clone MC.7.G5. Recognition was restricted to cells expressing a SNV of MR1 (MR1*04) and was not cancer-specific. 7G5.TCR-T and 7G5-like TCR-T cells reacted to both cancer and healthy cells endogenously expressing MR1*04 SNVs, which encode R9H and H17R substitutions. This allelic specificity could be overcome by expressing supraphysiological levels of the wild-type MR1 (MR1*01) in cell lines. Healthy individuals harbor T cells reactive to MR1 variants displaying self-ligands expressed in cancer and benign tissues. Described "cancer-specific" MR1-restricted TCRs need further validation, covering conserved allomorphs of MR1. Ligands require identification to ensure targeting MR1 is restricted to those specific to cancer and not normal tissues. For the wider field of immunology and transplant biology, the observation that MR1*04 may behave as an alloantigen warrants further study. .

NASTRA: accurate analysis of short tandem repeat markers by nanopore sequencing with repeat-structure-aware algorithm.

Short-tandem repeats (STRs) are the type of genetic markers extensively utilized in biomedical and forensic applications. Due to sequencing noise in nanopore sequencing, accurate analysis methods are lacking. We developed NASTRA, an innovative tool for Nanopore Autosomal Short Tandem Repeat Analysis, which overcomes traditional database-based methods' limitations and provides a precise germline analysis of STR genetic markers without the need for allele sequence reference. Demonstrating high accuracy in cell line authentication testing and paternity testing, NASTRA significantly surpasses existing methods in both speed and accuracy. This advancement makes it a promising solution for rapid cell line authentication and kinship testing, highlighting the potential of nanopore sequencing for in-field applications.

Targeted Sequencing of CTX-M Alleles Encoding for Extended-Spectrum Beta-Lactamases in Seattle Area Wastewater

Targeted Sequencing of CTX-M Alleles Encoding for Extended-Spectrum Beta-Lactamases in Seattle Area Wastewater

CRISPR-CLEAR: Nucleotide-Resolution Mapping of Regulatory Elements via Allelic Readout of Tiled Base Editing.

CRISPR tiling screens have advanced the identification and characterization of regulatory sequences but are limited by low resolution arising from the indirect readout of editing via guide RNA sequencing. This study introduces CRISPR-CLEAR , an end-to-end experimental assay and computational pipeline, which leverages targeted sequencing of CRISPR-introduced alleles at the endogenous target locus following dense base-editing mutagenesis. This approach enables the dissection of regulatory elements at nucleotide resolution, facilitating a direct assessment of genotype-phenotype effects.

Extending HLA allele sequences using next-generation sequencing technology.

Extended sequences for 13 HLA alleles were found which had limited coverage previously.

A Newly Incompatibility F Replicon Allele (FIB81) in Extensively Drug-Resistant Escherichia coli Isolated from Diseased Broilers.

Multiple drug resistance (MDR) has gained pronounced attention among Enterobacterales. The transfer of multiple antimicrobial resistance genes, frequently carried on conjugative incompatibility F (IncF) plasmids and facilitating interspecies resistance transmission, has been linked to Salmonella spp. and E. coli in broilers. In Egypt, the growing resistance is exacerbated by the limited clinical efficacy of many antimicrobials. In this study, IncF groups were screened and characterized in drug-resistant Salmonella spp. and E. coli isolated from broilers. The antimicrobial resistance profile, PCR-based replicon typing of bacterial isolates pre- and post-plasmid curing, and IncF replicon allele sequence typing were investigated. Five isolates of E. coli (5/31; 16.13%) and Salmonella spp. (5/36; 13.89%) were pan-susceptible to the examined antimicrobial agents, and 85.07% of tested isolates were MDR and extensively drug-resistant (XDR). Twelve MDR and XDR E. coli and Salmonella spp. isolates were examined for the existence of IncF replicons (FII, FIA, and FIB). They shared resistance to ampicillin, ampicillin/sulbactam, amoxicillin/clavulanate, doxycycline, cefotaxime, and colistin. All isolates carried from one to two IncF replicons. The FII-FIA-FIB+ and FII-FIA+FIB- were the predominant replicon patterns. FIB was the most frequently detected replicon after plasmid curing. Three XDR E. coli isolates that were resistant to 12-14 antimicrobials carried a newly FIB replicon allele with four nucleotide substitutions: C99→A, G112→T, C113→T, and G114→A. These findings suggest that broilers are a significant reservoir of IncF replicons with highly divergent IncF-FIB plasmid incompatibility groups circulating among XDR Enterobacterales. Supporting these data with additional comprehensive epidemiological studies involving replicons other than the IncF can provide insights for implementing efficient policies to prevent the spreading of new replicons to humans.

A Chinese individual with DEL phenotype caused by a novel RHD allele

BackgroundRhD variants are categorized into partial D, weak D, and DEL. The detection of DEL can only be achieved through the adsorption and elution method or molecular techniques. Here, we report a case of DEL phenotypes associated with a novel allele in a Chinese individual. Study design and methodsWe used serological methods such as saline, indirect anti-human globulin, and adsorption-elution. The RHD genotype was determined by the PCR-sequence specific primer (PCR-SSP) method as well as the Sanger dideoxy sequencing. ResultsRBCs of the sample were found to be DEL phenotype by serological testing, with negative reactions in the saline and indirect anti-human globulin tests while positive reactions by the absorption–elution method. The genotyping results revealed a hemizygous (RHDc .1127 T>G/RHD-). The novel allele sequence has been submitted to GenBank (Accession number: OR608456). ConclusionOur study demonstrates a case of a Chinese individual with DEL phenotype caused by a novel allele RHD c .1127 T > G. It expands the database of the DEL variant.

Molecular mechanisms underlying gene regulatory variation of maize metabolic traits.

Variation in gene expression levels is pervasive among individuals and races or varieties, and has substantial agronomic consequences, for example, by contributing to hybrid vigor. Gene expression level variation results from mutations in regulatory sequences (cis) and/or transcription factor (TF) activity (trans), but the mechanisms underlying cis- and/or trans-regulatory variation of complex phenotypes remain largely unknown. Here, we investigated gene expression variation mechanisms underlying the differential accumulation of the insecticidal compounds maysin and chlorogenic acid in silks of widely used maize (Zea mays) inbreds, B73 and A632. By combining transcriptomics and cistromics, we identified 1,338 silk direct targets of the maize R2R3-MYB TF Pericarp color1 (P1), consistent with it being a regulator of maysin and chlorogenic acid biosynthesis. Among these P1 targets, 464 showed allele-specific expression (ASE) between B73 and A632 silks. Allelic DNA-affinity purification sequencing identified 34 examples in which P1 allelic specific binding (ASB) correlated with cis-expression variation. From previous yeast one-hybrid studies, we identified 9 TFs potentially implicated in the control of P1 targets, with ASB to 83 out of 464 ASE genes (cis) and differential expression of 4 out of 9 TFs between B73 and A632 silks (trans). These results provide a molecular framework for understanding universal mechanisms underlying natural variation of gene expression levels, and how the regulation of metabolic diversity is established.

Sequencing-induced artefacts in NGS STR data

Significant progress has been made in recent years in the development of techniques for Next Generation Sequencing (NGS), or Massively Parallel Sequencing (MPS), of forensically relevant short tandem repeat (STR) loci. However, as these technologies are investigated and adopted by forensic laboratories, new challenges unfold that require further scrutiny. In the analysis of DNA profiles generated using the MiSeq FGx sequencing system, we have observed noise sequences with relatively high readcounts that are challenging to distinguish from genuine alleles. These high read count noise sequences appear as allele sequences with one or a few substituted bases compared to a known allele sequence within the profile.An examination of ForenSeq DNA Signature Prep Kit STR noise sequences revealed that the substituted base of a parent allele can align to the same position on the sequence across noise sequences. This suggests that these substitution events occur at specific positions within the amplicon, resulting in multiple noise reads with substitutions at the same position. Mapping of the noise events onto the original raw read positions revealed a high number of events, or “noise spikes”, occurring atspecific positions within a given sequencing run. These noise spikes affected reads across the entire run, agnostic of locus or sample, while the position, occurrence, and amplitude of the spikes differed across runs. The majority of noise sequences with high read counts in a DNA profile were generated from base changes at these spike positions, and could be classified as “noise spike artefacts”.In this paper we present evidence of the noise spike artefacts and their genesis during the sequencing process in the sequencing-by-synthesis (SBS) cycles, as well as the methods developed to detect them. The information and methods will assist laboratories with detecting noise spikes in MiSeq FGx sequencing runs, differentiating authentic allele sequences from noise spike artefacts, and developing protocols for analyst review and handling of MiSeq FGx data.

25 years of the IPD-IMGT/HLA Database.

Twenty-five years ago, in 1998, the HLA Informatics Group of the Anthony Nolan Research Institute released the IMGT/HLA Database. Since this time, this online resource has acted as the repository for the numerous variant sequences of HLA alleles named by the WHO Nomenclature Committee for Factors of the HLA System. The IPD-IMGT/HLA Database has provided a stable, highly accessible, user-friendly repository for this work. During this time, the technology underlying HLA typing has undergone significant changes. Next generation sequencing (NGS) has superseded previous methodologies of HLA typing and can generate large amounts of high-resolution sequencing data. This has resulted in a drastic increase in the number and complexity of sequences submitted to the database. The challenge for the IPD-IMGT/HLA Database has been to maintain the highest standards of curation, while supporting the core set of tools and functionality to our users with increased numbers of submissions and sequences. Traditional methods of accessing and presenting data have been challenged and new methods utilising new computing technologies have had to be developed to keep pace and support a shifting user demographic.

Variability and Number of Circulating Complementary Sex Determiner (Csd) Alleles in a Breeding Population of Italian Honeybees under Controlled Mating.

In Apis mellifera, csd is the primary gene involved in sex determination: haploid hemizygous eggs develop as drones, while females develop from eggs heterozygous for the csd gene. If diploid eggs are homozygous for the csd gene, diploid drones will develop, but will be eaten by worker bees before they are born. Therefore, high csd allelic diversity is a priority for colony survival and breeding. This study aims to investigate the variability of the hypervariable region (HVR) of the csd gene in bees sampled in an apiary under a selection scheme. To this end, an existing dataset of 100 whole-genome sequences was analyzed with a validated pipeline based on de novo assembly of sequences within the HVR region. In total, 102 allelic sequences were reconstructed and translated into amino acid sequences. Among these, 47 different alleles were identified, 44 of which had previously been observed, while 3 are novel alleles. The results show a high variability in the csd region in this breeding population of honeybees.

Human restricted CHRFAM7A gene increases brain efficiency.

CHRFAM7A, a uniquely human fusion gene, has been associated with neuropsychiatric disorders including Alzheimer's disease, schizophrenia, anxiety, and attention deficit disorder. Understanding the physiological function of CHRFAM7A in the human brain is the first step to uncovering its role in disease. CHRFAM7A was identified as a potent modulator of intracellular calcium and an upstream regulator of Rac1 leading to actin cytoskeleton reorganization and a switch from filopodia to lamellipodia implicating a more efficient neuronal structure. We performed a neurocognitive-MRI correlation exploratory study on 46 normal human subjects to explore the effect of CHRFAM7A on human brain. Dual locus specific genotyping of CHRFAM7A was performed on genomic DNA to determine copy number (TaqMan assay) and orientation (capillary sequencing) of the CHRFAM7A alleles. As only the direct allele is expressed at the protein level and affects α7 nAChR function, direct allele carriers and non-carriers are compared for neuropsychological and MRI measures. Subjects underwent neuropsychological testing to measure motor (Timed 25-foot walk test, 9-hole peg test), cognitive processing speed (Symbol Digit Modalities Test), Learning and memory (California Verbal Learning Test immediate and delayed recall, Brief Visuospatial Memory Test-Revised immediate and delayed recall) and Beck Depression Inventory-Fast Screen, Fatigue Severity Scale. All subjects underwent MRI scanning on the same 3 T GE scanner using the same protocol. Global and tissue-specific volumes were determined using validated cross-sectional algorithms including FSL's Structural Image Evaluation, using Normalization, of Atrophy (SIENAX) and FSL's Integrated Registration and Segmentation Tool (FIRST) on lesion-inpainted images. The cognitive tests were age and years of education-adjusted using analysis of covariance (ANCOVA). Age-adjusted analysis of covariance (ANCOVA) was performed on the MRI data. CHRFAM7A direct allele carrier and non-carrier groups included 33 and 13 individuals, respectively. Demographic variables (age and years of education) were comparable. CHRFAM7A direct allele carriers demonstrated an upward shift in cognitive performance including cognitive processing speed, learning and memory, reaching statistical significance in visual immediate recall (FDR corrected p = 0.018). The shift in cognitive performance was associated with smaller whole brain volume (uncorrected p = 0.046) and lower connectivity by resting state functional MRI in the visual network (FDR corrected p = 0.027) accentuating the cognitive findings. These data suggest that direct allele carriers harbor a more efficient brain consistent with the cellular biology of actin cytoskeleton and synaptic gain of function. Further larger human studies of cognitive measures correlated with MRI and functional imaging are needed to decipher the impact of CHRFAM7A on brain function.

Nanopore sequencing in distinguishing between wild-type and vaccine strains of Varicella-Zoster virus

BackgroundThe introduction of varicella vaccines into routine pediatric immunization programs has led to a considerable reduction in varicella incidence. However, there have been reports of varicella, herpes zoster, and meningitis caused by the vaccine strain of varicella-zoster virus (VZV), raising concerns. Establishing the relationship between the wild-type and vaccine strains in VZV infections among previously vaccinated individuals is crucial. Differences in the single nucleotide polymorphisms (SNPs) among vaccine strains can be utilized to identify the strain. In this study, we employed nanopore sequencing to identify VZV strains and analyzed clinical samples. MethodsWe retrospectively examined vesicle and cerebrospinal fluid samples from patients with VZV infections. One sample each of the wild-type and vaccine strains, previously identified using allelic discrimination real-time PCR and direct sequencing, served as controls. Ten samples with undetermined VZV strains were included. After DNA extraction, a long PCR targeting the VZV ORF62 region was executed. Nanopore sequencing identified SNPs, allowing discrimination between the vaccine and wild-type strains. ResultsNanopore sequencing confirmed SNPs at previously reported sites (105,705, 106,262, 107,136, and 107,252), aiding in distinguishing between wild-type and vaccine strains. Among the ten unknown samples, nine were characterized as wild strains and one as a vaccine strain. Even in samples with low VZV DNA levels, nanopore sequencing was effective in strain identification. ConclusionThis study validates that nanopore sequencing is a reliable method for differentiating between the wild-type and vaccine strains of VZV. Its ability to produce long-read sequences is remarkable, allowing simultaneous confirmation of known SNPs and the detection of new mutations. Nanopore sequencing can serve as a valuable tool for the swift and precise identification of wild-type and vaccine strains and has potential applications in future VZV surveillance.

Abstract 1612: Timing ecDNA acquisition across cancer

Abstract Extrachromosomal DNA (ecDNA) are acentric, circular DNA fragments which form after DNA double strand breaks. They segregate randomly during mitosis, promoting high-level amplifications of positively selected oncogenes such as EGFR and MYC. The prevalence of ecDNA across several cancer types, such as glioblastoma and sarcoma, has made promising therapeutic targets. However, while DNA damage is a hallmark of nearly all cancers, the formation of ecDNA is not, and the determinants for ecDNA acquisition are unclear. Furthermore, whether the timing of tumors’ acquisition of ecDNA (i.e., early or late clonal or sub-clonal) affects patient survival and drug resistance is unclear. Here, we leverage orthoganol information from ecDNA amplicons, such as mutation allele frequency and sequence synteny, to develop computational methods for elucidating the evolutionary history of ecDNA. Applying these to several cancer subtypes in TCGA, we find that most ecDNAs arise early in tumor evolution. We expect these results will help identify determinants of ecDNA acquisition and its temporal influence on patient survival. Citation Format: Andrew R. Lynch, Peter Van Loo. Timing ecDNA acquisition across cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1612.

Sequence analysis and secondary structure prediction of autosomal STR alleles using next generation sequencing (NGS) data

The inception of NGS technology in forensic DNA analysis explores the sequence-based alleles in STR markers. This allows the in-depth analysis of observed STR allelic sequences to understand their nature, stability, inheritance, and possible generation of artifacts during their analysis. In the current study, 100 allelic sequences at 20 STR markers generated using the NGS technique were analyzed for their sequence topographies and prediction of secondary structures. The G + C content of the alleles observed in 20 STR markers used in this study varied from 58.65 ± 1.367% (D2S1338) to 7.62 ± 0.844% (D12ATA63). The average exact mass of one stand and for the opposite stand mass was found to be highest in FGA (25,963.248 ± 1623.271; 27,501.720 ± 1712.691), whereas, D4S2408 generated the lowest exact mass of one stand and for the opposite stand mass (11,136.354 ± 1521.757; 11,582.021 ± 1585.486). As expected, none of the STR markers showed the presence of open reading frames, Codons, and CRISPR sequences. Three STR markers viz. D2S1338, TH01, and D5S2800 showed the presence of restriction sites for Cac8I, TspGWI, TspDTI, AccI, and Hpy8I enzymes. Phylogenetic analysis reveals the close association of alleles between the D12ATA63 and D19S433 markers. Stable pseudoknots were predicted at alleles of D2S1338 showing an average energy of −0.76 with the highest number of nucleotides present in the pseudoknots i.e., 21.33, suggesting this marker is more prone to generate amplification artifacts.

Full-length sequence of the novel allele, HLA-B*58:01:40, identified in a Chinese individual.

HLA-B*58:01:40 differs from HLA-B*58:01:01 by a single nucleotide change in exon 3, 507 C- > T (codon 145.3 CGC- > CGT).

Use of a Candida albicans SC5314 PacBio HiFi reads dataset to close gaps in the reference genome assembly, reveal a subtelomeric gene family, and produce accurate phased allelic sequences.

Candida albicans SC5314 is the most-often used strain for molecular manipulation of the species. The SC5314 reference genome sequence is the result of considerable effort from many scientists and has advanced research into fungal biology and pathogenesis. Although the resource is highly developed and presented in a phased diploid format, the sequence includes gaps and does not extend to the telomeres on its eight chromosome pairs. Accurate SC5314 genome assembly is complicated by the presence of extensive repeated sequences and considerable allelic length variation at some loci. Advances in genome sequencing technology provide the tools to obtain highly accurate long-read data that span even the most-difficult-to-assemble genome regions. Here, we describe derivation of a PacBio HiFi data set and creation of a collapsed haploid telomere-to-telomere assembly of the SC5314 genome (ASM3268872v1) that revealed previously unknown features of the strain. ASM3268872v1 subtelomeric distances were up to 19 kb larger than in the reference genome and revealed a family of highly conserved DNA helicase-encoding genes at 10 of the 16 chromosome ends. We also describe alignments of individual HiFi reads to deduce accurate diploid sequences for the most notoriously difficult-to-assemble C. albicans genes: the agglutinin-like sequence (ALS) gene family. We provide a tutorial that demonstrates how the HiFi reads can be visualized to explore any region of interest. Availability of the HiFi reads data set and the ASM3268872v1 comparative guide assembly will streamline research efforts because accurate diploid sequences can be derived using simple in silico methods rather than time-consuming laboratory-bench approaches.

Ancient tree genomes for old questions.

Most foundational work on the evolution and migration of plant species relies on genomic data from contemporary samples. Ancient plant samples can give us access to allele sequences and distributions on the landscape dating back to the mid Holocene or earlier (Gugerli etal.,2005). Nuclear DNA from ancient wood, however, has been mostly inaccessible until now. In a From the Cover article in this issue of Molecular Ecology, Wagner etal. (2023) present the first resequenced nuclear genomes from ancient oak wood, including two samples dated to the 15th century and one that dates to more than 3500 years ago. These ancient tree genomes open the possibility for investigating species adaptation, migration, divergence, and hybridisation in the deep past. They pave the way for what we hope will be a new era in the use of paleogenomics to study Holocene tree histories.