Pangenome graph construction from genome alignments with Minigraph-Cactus

Pangenome references address biases of reference genomes by storing a representative set of diverse haplotypes and their alignment, usually as a graph. Alternate alleles determined by variant callers can be used to construct pangenome graphs, but advances in long-read sequencing are leading to widely available, high-quality phased assemblies. Constructing a pangenome graph directly from assemblies, as opposed to variant calls, leverages the graph's ability to represent variation at different scales. Here we present the Minigraph-Cactus pangenome pipeline, which creates pangenomes directly from whole-genome alignments, and demonstrate its ability to scale to 90 human haplotypes from the Human Pangenome Reference Consortium. The method builds graphs containing all forms of genetic variation while allowing use of current mapping and genotyping tools. We measure the effect of the quality and completeness of reference genomes used for analysis within the pangenomes and show that using the CHM13 reference from the Telomere-to-Telomere Consortium improves the accuracy of our methods. We also demonstrate construction of a Drosophila melanogaster pangenome.

Expanding deep phenotypic spectrum associated with atypical pathogenic structural variations overlapping 15q11-q13 imprinting region

BACKGROUND

The 15q11-q13 region is a genetic locus with genes subject to genomic imprinting, significantly influencing neurodevelopment. Genomic imprinting is an epigenetic phenomenon that causes differential gene expression based on the parent of origin. In most diploid organisms, gene expression typically involves an equal contribution from both maternal and paternal alleles, shaping the phenotype. Nevertheless, in mammals, including humans, mice, and marsupials, the functional equivalence of parental alleles is not universally maintained. Notably, during male and female gametogenesis, parental alleles may undergo differential marking or imprinting, thereby modifying gene expression without altering the underlying DNA sequence. Neurodevelopmental disorders, such as Prader-Willi syndrome (PWS) (resulting from the absence of paternally expressed genes in this region), Angelman syndrome (AS) (associated with the absence of the maternally expressed UBE3A gene), and 15q11-q13 duplication syndrome (resulting from the two common forms of duplications-either an extra isodicentric 15 chromosome or an interstitial 15 duplication), are the outcomes of genetic variations in this imprinting region.

METHODS

Conducted a genomic study to identify the frequency of pathogenic variants impacting the 15q11-q13 region in an ethnically homogenous population from Bangladesh. Screened all known disorders from the DECIPHER database and identified variant enrichment within this cohort. Using the Horizon analysis platform, performed enrichment analysis, requiring at least >60% overlap between a copy number variation and a disorder breakpoint. Deep clinical phenotyping was carried out through multiple examination sessions to evaluate a range of clinical symptoms.

RESULTS

This study included eight individuals with clinically suspected PWS/AS, all previously confirmed through chromosomal microarray analysis, which revealed chromosomal breakpoints within the 15q11-q13 region. Among this cohort, six cases (75%) exhibited variable lengths of deletions, whereas two cases (25%) showed duplications. These included one type 2 duplication, one larger atypical duplication, one shorter type 2 deletion, one larger type 1 deletion, and four cases with atypical deletions. Furthermore, thorough clinical assessments led to the diagnosis of four PWS patients, two AS patients, and two individuals with 15q11-q13 duplication syndrome.

CONCLUSION

Our deep phenotypic observations identified a spectrum of clinical features that overlap and are unique to PWS, AS, and Dup15q syndromes. Our findings establish genotype-phenotype correlation for patients impacted by variable structural variations within the 15q11-q13 region.

Ancestral allele of DNA polymerase gamma modifies antiviral tolerance

Mitochondria are critical modulators of antiviral tolerance through the release of mitochondrial RNA and DNA (mtDNA and mtRNA) fragments into the cytoplasm after infection, activating virus sensors and type-I interferon (IFN-I) response1-4. The relevance of these mechanisms for mitochondrial diseases remains understudied. Here we investigated mitochondrial recessive ataxia syndrome (MIRAS), which is caused by a common European founder mutation in DNA polymerase gamma (POLG1)5. Patients homozygous for the MIRAS variant p.W748S show exceptionally variable ages of onset and symptoms5, indicating that unknown modifying factors contribute to disease manifestation. We report that the mtDNA replicase POLG1 has a role in antiviral defence mechanisms to double-stranded DNA and positive-strand RNA virus infections (HSV-1, TBEV and SARS-CoV-2), and its p.W748S variant dampens innate immune responses. Our patient and knock-in mouse data show that p.W748S compromises mtDNA replisome stability, causing mtDNA depletion, aggravated by virus infection. Low mtDNA and mtRNA release into the cytoplasm and a slow IFN response in MIRAS offer viruses an early replicative advantage, leading to an augmented pro-inflammatory response, a subacute loss of GABAergic neurons and liver inflammation and necrosis. A population databank of around 300,000 Finnish individuals6 demonstrates enrichment of immunodeficient traits in carriers of the POLG1 p.W748S mutation. Our evidence suggests that POLG1 defects compromise antiviral tolerance, triggering epilepsy and liver disease. The finding has important implications for the mitochondrial disease spectrum, including epilepsy, ataxia and parkinsonism.

Detection of the HLA-C*07:04:29 allele in a Chinese individual

HLA-C*07:04:29 differs from HLA-C*07:04:01:01 by a single substitution in exon 4.

Characterization of the novel HLA-C*06:376N allele by Pacific Biosciences HiFi sequencing in a Chinese individual

HLA-C*06:376N differs from HLA-C*06:02:01:01 by seven nucleotide changes in exon 2, intron 2, and exon 3.

Identification of the novel HLA-B*13:191 allele by next-generation sequencing

The novel HLA-B*13:191 allele was detected during the HLA typing for kidney transplantation.

The novel HLA-DRB1*11:323 allele characterized by two different sequencing-based typing techniques

The novel allele HLA-DRB1*11:323 differs from HLA-DRB1*11:01:02:01 by one non-synonymous nucleotide substitution in exon 2.

The pseudoentropy of allele frequency trajectories, the persistence of variation, and the effective population size

To concisely describe how genetic variation, at individual loci or across whole genomes, changes over time, and to follow transitory allelic changes, we introduce a quantity related to entropy, that we term pseudoentropy. This quantity emerges in a diffusion analysis of the mean time a mutation segregates in a population. For a neutral locus with an arbitrary number of alleles, the mean time of segregation is generally proportional to the pseudoentropy of initial allele frequencies. After the initial time point, pseudoentropy generally decreases, but other behaviours are possible, depending on the genetic diversity and selective forces present. For a biallelic locus, pseudoentropy and entropy coincide, but they are distinct quantities with more than two alleles. Thus for populations with multiple biallelic loci, the language of entropy suffices. Then entropy, combined across loci, serves as a concise description of genetic variation. We used individual based simulations to explore how this entropy behaves under different evolutionary scenarios. In agreement with predictions, the entropy associated with unlinked neutral loci decreases over time. However, deviations from free recombination and neutrality have clear and informative effects on the entropy's behaviour over time. Analysis of publicly available data of a natural D. melanogaster population, that had been sampled over seven years, using a sliding-window approach, yielded considerable variation in entropy trajectories of different genomic regions. These mostly follow a pattern that suggests a substantial effective population size and a limited effect of positive selection on genome-wide diversity over short time scales.

Characterisation of the novel HLA-B*58:02:04 allele by next-generation sequencing

HLA-B*58:02:04 differs from HLA-B*58:02:01 by one synonymous nucleotide in codon 215 in exon 4.

The Association Between the 5-Hydroxytryptamine Receptor 2A Gene Variants rs6311 and rs6313 and Obstructive Sleep Apnea in the Iranian Kurdish Population

Introduction: Sleep is one of the most significant parts of everyone's life. Most people sleep for about one-third of their lives. Sleep disorders negatively impact the quality of life. Obstructive sleep apnea (OSA) is a severe sleep disorder that significantly impacts the patient's life and their family members. This study aimed to investigate the relationship between rs6313 and rs6311 polymorphisms in the serotonin receptor type 2A gene and OSA in the Kurdish population. Materials and Methods: The study's population comprises 100 OSA sufferers and 100 healthy people. Polysomnography diagnostic tests were done on both the patient and control groups. The polymerase chain reaction (PCR)-restriction fragment length polymorphism (RFLP) was used to investigate the relationship between OSA and LEPR gene polymorphisms. Results: Statistical analysis showed a significant relationship between genotype frequencies of patient and control groups of rs6311 with OSA in dominant [odds ratio (OR) = 5.203, p < 0.001) and codominant models (OR = 9.7, p < 0.001). Also, there was a significant relationship between genotype frequencies of patient and control groups of rs6313 with OSA in dominant (OR = 10.565, p < 0.001) and codominant models (OR = 5.938, p < 0.001). Conclusions: Findings from the study demonstrated that the two polymorphisms rs6311 and rs6313 could be effective at causing OSA; however, there was no correlation between the severity of the disease and either of the two polymorphisms.

The novel HLA-DPB1*1500:01N allele, identified by Sanger dideoxy nucleotide sequencing in a Chinese individual

HLA-DPB1*1500:01N differs from HLA-DPB1*05:01:01:01 by one nucleotide in exon 3.

Identification of the novel HLA-C allele, HLA-C*07:02:150

A novel HLA-C*07 allele, now officially designated HLA-C*07:02:150, was identified by next-generation sequencing.

Identification of two novel HLA-DRB1 alleles, HLA-DRB1*08:03:13 and HLA-DRB1*08:119

Compared with HLA-DRB1*08:03:02:01, the alleles HLA-DRB1*08:03:13 and HLA-DRB1*08:119 each show one nucleotide substitution, respectively.

The association of CYP2C19 LoF alleles with adverse clinical outcomes in stroke patients taking clopidogrel: An updated meta-analysis

The aggregated risk of recurrent stroke in stroke/transient ischemic attack (TIA) patients carrying CYP2C19 LoF alleles who take clopidogrel has not been investigated recently, and the available research is limited. This study aimed to perform an updated meta-analysis to assess the association between CYP2C19 LoF alleles and the risk of recurrent stroke in patients taking clopidogrel. Databases were searched for the literature on eligible studies. The end points were recurrent stroke, composite vascular events, and bleeding events. Odds ratios (ORs) were calculated using RevMan software, where p < 0.05 was considered statistically significant. Patients carrying CYP2C19 LoF alleles who were treated with clopidogrel had a significantly increased risk of recurrent ischemic stroke compared with non-carriers (OR 2.18, 96% CI 1.80-2.63; p < 0.00001). The risk of recurrent stroke was only significantly different in Asian patients (OR 2.29, 96% CI 1.88-2.80; p < 0.00001) but not in patients of other ethnicities; however, there were a limited number of studies in other ethnic groups. Both observational studies (OR 2.83, 96% CI 2.20-3.65; p < 0.00001) and RCTs (OR 1.48, 96% CI 1.10-1.98; p = 0.009) found associations with a significantly increased risk of recurrent ischemic stroke. Asian stroke patients or TIA patients carrying CYP2C19 LoF alleles and taking clopidogrel were at a significantly higher risk of recurrent ischemic stroke than non-carriers. Significantly increased risk of recurrent ischemic stroke was found in both observational studies and RCTs.

The novel HLA-DPA1*02:117 allele, identified by Sanger dideoxy nucleotide sequencing in a Chinese individual

HLA-DPA1*02:117 differs from HLA-DPA1*02:02:02:01 by one nucleotide in exon 2.

[Genetic Diversity of Blastocystis in Diarrheal Cases: Identification of Subtypes and Alleles]

Blastocystis spp. are the most common intestinal protozoan parasites detected in human stool samples. While identified long before today, its pathogenicity remains controversial. It is generally asymptomatic but in symptomatic cases, many gastrointestinal symptoms, especially diarrhea, have been associated with Blastocystis infection. In recent years, the relationship between the symptoms observed in cases and Blastocystis subtypes (ST) has been reported. The aim of this study was to detect Blastocystis in diarrheal cases admitted to the Aydın Adnan Menderes University Faculty of Medicine, Department of Parasitology Laboratory, to determine subtypes and allele diversity and to investigate its relationship with clinical symptoms. For this purpose, diarrheal stool samples of 200 cases were included in the study and their demographic characteristics (age, gender, residence) and clinical findings (abdominal pain, dyspepsia, nausea-vomiting, weakness, weight loss, anal itching, rash, urticaria) were recorded. Blastocystis was detected by direct microscope method (DM) and by molecular analyses which were performed with polymerase chain reaction (PCR). Subtype diversity was determined based on DNA sequence analysis by PCR targeting the Blastocystis ribosomal ribonucleic acid small subunit (SSU rRNA) gene. In addition, alleles related to Blastocystis subtypes were determined and statistically compared between all data and clinical findings. In the current study, Blastocystis was detected in 31 (15.5%) samples by DM and in 35 (17.5%) samples by PCR specific to the Blastocystis SSU rRNA gene among 200 diarrheal stool samples. No statistical difference was detected between Blastocystis and demographic characteristics. Dyspepsia and nausea-vomiting symptoms differed significantly in cases with Blastocystis compared to negative ones (p= 0.0025, p= 0.0498). Blastocystis subtype was detected in 33 samples by SSU rRNA sequence analysis, and the subtype distribution was ST1 (n= 10, 30.3%), ST2 (n= 4, 12.1%) and ST3 (n= 19, 57.6%). In the statistical evaluation between clinical findings and Blastocystis subtypes, a relationship was found between dyspepsia and Blastocystis ST3 (p= 0.0039). The allele diversity of Blastocystis subtypes was determined as allele 4 (10/10) in all ST1, allele 11 (2/4) and 12 (2/4) in ST2, allele 34 (14/19), 36 (4/19), and 38 (1/19) in ST3. In conclusion, our study provides important data on the molecular epidemiological characteristics of the Blastocystis by determining positivity, subtypes and alleles in diarrheal cases. Therefore, within the scope of the one health approach, comprehensive molecular epidemiological studies are required to determine the presence and genotypes of Blastocystis in human, animal and environmental samples.

High-resolution KIR and HLA genotyping in three Chinese ethnic minorities reveals distinct origins

Polymorphism of killer-cell immunoglobulin-like receptors (KIRs) and their HLA class I ligands impacts the effector activity of cytotoxic NK cell and T cell subsets. Therefore, understanding the extent and implications of KIR and HLA class I genetic polymorphism across various populations is important for immunological and medical research. In this study, we conducted a high-resolution investigation of KIR and HLA class I diversity in three distinct Chinese ethnic minority populations. We studied the She, Yugur, and Tajik, and compared them with the Zhejiang Han population (Zhe), which represents the majority Southern Han ethnicity. Our findings revealed that the Tajik population exhibited the most diverse KIR copy number, allele, and haplotype diversity among the four populations. This diversity aligns with their proposed ancestral origin, closely resembling that of Iranian populations, with a relatively higher presence of KIR-B genes, alleles, and haplotypes compared with the other Chinese populations. The Yugur population displayed KIR distributions similar to those of the Tibetans and Southeast Asians, whereas the She population resembled the Zhe and other East Asians, as confirmed by genetic distance analysis of KIR. Additionally, we identified 12.9% of individuals across the three minority populations as having KIR haplotypes characterized by specific gene block insertions or deletions. Genetic analysis based on HLA alleles yielded consistent results, even though there were extensive variations in HLA alleles. The observed variations in KIR interactions, such as higher numbers of 2DL1-C2 interactions in Tajik and Yugur populations and of 2DL3-C1 interactions in the She population, are likely shaped by demographic and evolutionary mechanisms specific to their local environments. Overall, our findings offer valuable insights into the distribution of KIR and HLA diversity among three distinct Chinese ethnic minority populations, which can inform future clinical and population studies.

Tissue-specific enhancer-gene maps from multimodal single-cell data identify causal disease alleles

Translating genome-wide association study (GWAS) loci into causal variants and genes requires accurate cell-type-specific enhancer-gene maps from disease-relevant tissues. Building enhancer-gene maps is essential but challenging with current experimental methods in primary human tissues. Here we developed a nonparametric statistical method, SCENT (single-cell enhancer target gene mapping), that models association between enhancer chromatin accessibility and gene expression in single-cell or nucleus multimodal RNA sequencing and ATAC sequencing data. We applied SCENT to 9 multimodal datasets including >120,000 single cells or nuclei and created 23 cell-type-specific enhancer-gene maps. These maps were highly enriched for causal variants in expression quantitative loci and GWAS for 1,143 diseases and traits. We identified likely causal genes for both common and rare diseases and linked somatic mutation hotspots to target genes. We demonstrate that application of SCENT to multimodal data from disease-relevant human tissue enables the scalable construction of accurate cell-type-specific enhancer-gene maps, essential for defining noncoding variant function.

Selfish conflict underlies RNA-mediated parent-of-origin effects

Genomic imprinting-the non-equivalence of maternal and paternal genomes-is a critical process that has evolved independently in many plant and mammalian species1,2. According to kinship theory, imprinting is the inevitable consequence of conflictive selective forces acting on differentially expressed parental alleles3,4. Yet, how these epigenetic differences evolve in the first place is poorly understood3,5,6. Here we report the identification and molecular dissection of a parent-of-origin effect on gene expression that might help to clarify this fundamental question. Toxin-antidote elements (TAs) are selfish elements that spread in populations by poisoning non-carrier individuals7-9. In reciprocal crosses between two Caenorhabditis tropicalis wild isolates, we found that the slow-1/grow-1 TA is specifically inactive when paternally inherited. This parent-of-origin effect stems from transcriptional repression of the slow-1 toxin by the PIWI-interacting RNA (piRNA) host defence pathway. The repression requires PIWI Argonaute and SET-32 histone methyltransferase activities and is transgenerationally inherited via small RNAs. Remarkably, when slow-1/grow-1 is maternally inherited, slow-1 repression is halted by a translation-independent role of its maternal mRNA. That is, slow-1 transcripts loaded into eggs-but not SLOW-1 protein-are necessary and sufficient to counteract piRNA-mediated repression. Our findings show that parent-of-origin effects can evolve by co-option of the piRNA pathway and hinder the spread of selfish genes that require sex for their propagation.

Detection of the novel HLA-B*55:01:31, HLA-C*07:1113 alleles and confirmation of HLA-C*12:392

Novel HLA-B*55:01:31, HLA-C*07:1113 alleles and confirmatory HLA-C*12:392 allele were detected during the HLA typing process.

The novel HLA-DRB1*03:210 allele characterised by two different sequencing-based typing techniques

The novel allele HLA-DRB1*03:210 differs from HLA-DRB1*03:01:01:01 by one non-synonymous nucleotide substitution in exon 3.

Identification of the novel HLA-DPA1*02:02:15 allele by next-generation sequencing

The novel HLA-DPA1*02:02:15 allele differs from HLA-DPA1*02:02:02:01 by one nucleotide substitution in exon 1.

Interpreting population- and family-based genome-wide association studies in the presence of confounding

A central aim of genome-wide association studies (GWASs) is to estimate direct genetic effects: the causal effects on an individual's phenotype of the alleles that they carry. However, estimates of direct effects can be subject to genetic and environmental confounding and can also absorb the "indirect" genetic effects of relatives' genotypes. Recently, an important development in controlling for these confounds has been the use of within-family GWASs, which, because of the randomness of mendelian segregation within pedigrees, are often interpreted as producing unbiased estimates of direct effects. Here, we present a general theoretical analysis of the influence of confounding in standard population-based and within-family GWASs. We show that, contrary to common interpretation, family-based estimates of direct effects can be biased by genetic confounding. In humans, such biases will often be small per-locus, but can be compounded when effect-size estimates are used in polygenic scores (PGSs). We illustrate the influence of genetic confounding on population- and family-based estimates of direct effects using models of assortative mating, population stratification, and stabilizing selection on GWAS traits. We further show how family-based estimates of indirect genetic effects, based on comparisons of parentally transmitted and untransmitted alleles, can suffer substantial genetic confounding. We conclude that, while family-based studies have placed GWAS estimation on a more rigorous footing, they carry subtle issues of interpretation that arise from confounding.

Novel regulatory variant in ABO intronic RUNX1 binding site inducing A3 phenotype

BACKGROUND AND OBJECTIVES

Mixed-field agglutination in ABO phenotyping (A3, B3) has been linked to genetically different blood cell populations such as in chimerism, or to rare variants in either ABO exon 7 or regulatory regions. Clarification of such cases is challenging and would greatly benefit from sequencing technologies that allow resolving full-gene haplotypes at high resolution.

MATERIALS AND METHODS

We used long-read sequencing by Oxford Nanopore Technologies to sequence the entire ABO gene, amplified in two overlapping long-range PCR fragments, in a blood donor presented with A3B phenotype. Confirmation analyses were carried out by Sanger sequencing and included samples from other family members.

RESULTS

Our data revealed a novel heterozygous g.10924C>A variant on the ABO*A allele located in the transcription factor binding site for RUNX1 in intron 1 (+5.8 kb site). Inheritance was shown by the results of the donor's mother, who shared the novel variant and the anti-A specific mixed-field agglutination.

CONCLUSION

We discovered a regulatory variant in the 8-bp RUNX1 motif of ABO, which extends current knowledge of three other variants affecting the same motif and also leading to A3 or B3 phenotypes. Overall, long-range PCR combined with nanopore sequencing proved powerful and showed great potential as an emerging strategy for resolving cases with cryptic ABO phenotypes.

Identification of the novel HLA-C allele, HLA-C*12:351Q, using next-generation sequencing

The novel allele HLA-C*12:351Q differs from HLA-C*12:02:02:01 by a single nucleotide deletion in exon 5.

Saliva direct PCR protocol for HLA-DQB1*02 genotyping

Celiac disease (CD) is an immune disorder, that is triggered by gluten ingestion in genetically predisposed individuals. The HLA-DQB1*02 allele is the main predisposing genetic factor and a candidate for first-line genotyping screening. We designed and validated a simple, DNA purification-free PCR protocol directly from crude saliva, enabling the detection of the DQB1*02 allele. This assay also distinguishes homozygous from heterozygous carriers. We propose this method for use in mass screening and/or epidemiological studies.

Network-based analysis predicts interacting genetic modifiers from a meta-mapping study of spike-wave discharge in mice

Absence seizures are characterized by brief lapses in awareness accompanied by a hallmark spike-and-wave discharge (SWD) electroencephalographic pattern and are common to genetic generalized epilepsies (GGEs). While numerous genes have been associated with increased risk, including some Mendelian forms with a single causal allele, most cases of GGE are idiopathic and there are many unknown genetic modifiers of GGE influencing risk and severity. In a previous meta-mapping study, crosses between transgenic C57BL/6 and C3HeB/FeJ strains, each carrying one of three SWD-causing mutations (Gabrg2tm1Spet(R43Q) , Scn8a8j or Gria4spkw1 ), demonstrated an antagonistic epistatic interaction between loci on mouse chromosomes 2 and 7 influencing SWD. These results implicate universal modifiers in the B6 background that mitigate SWD severity through a common pathway, independent of the causal mutation. In this study, we prioritized candidate modifiers in these interacting loci. Our approach integrated human genome-wide association results with gene interaction networks and mouse brain gene expression to prioritize candidate genes and pathways driving variation in SWD outcomes. We considered candidate genes that are functionally associated with human GGE risk genes and genes with evidence for coding or non-coding allele effects between the B6 and C3H backgrounds. Our analyses output a summary ranking of gene pairs, one gene from each locus, as candidates for explaining the epistatic interaction. Our top-ranking gene pairs implicate microtubule function, cytoskeletal stability and cell cycle regulation as novel hypotheses about the source of SWD variation across strain backgrounds, which could clarify underlying mechanisms driving differences in GGE severity in humans.

Determination for KIR genotype and allele copy number via real-time quantitative PCR method

Killer cell immunoglobulin-like receptor (KIR) and human leukocyte antigen (HLA) play crucial roles in regulating NK cell activity. Here, we report a real-time quantitative PCR (qPCR) to genotype all KIR genes and their copy numbers simultaneously. With 18 pairs of locus-specific primers, we identified KIR genes by Ct values and determined KIR copy number using the 2-∆Ct method. Haplotypes were assigned based on KIR gene copy numbers. The real-time qPCR results were consistent with the NGS method, except for one sample with KIR2DL5 discrepancy. qPCR is a multiplex method that can identify KIR copy number, which helps obtain a relatively accurate haplotype structure, facilitating increased KIR research in laboratories where NGS or other high-resolution methods are not available.

Creation of Knock-In Alleles of Insulin Receptor Tagged by Fluorescent Proteins mCherry or EYFP in Fruit Fly Drosophila melanogaster

The insulin/insulin-like growth factor-like signaling (IIS) pathway is highly conserved across metazoans and regulates numerous physiological functions, including development, metabolism, fecundity, and lifespan. The insulin receptor (InR), a crucial membrane receptor in the IIS pathway, is known to be ubiquitously expressed in various tissues, albeit at generally low levels, and its subcellular localization remains incompletely characterized. In this study, we employed CRISPR-mediated mutagenesis in the fruit fly Drosophila to create knock-in alleles of InR tagged with fluorescent proteins (InR::mCherry or InR::EYFP). By inserting the coding sequence of the fluorescent proteins mCherry or EYFP near the end of the coding sequence of the endogenous InR gene, we could trace the natural InR protein through their fluorescence. As an example, we investigated epithelial cells of the male accessory gland (AG), an internal reproductive organ, and identified two distinct patterns of InR::mCherry localization. In young AG, InR::mCherry accumulated on the basal plasma membrane between cells, whereas in mature AG, it exhibited intracellular localization as multiple puncta, indicating endocytic recycling of InR during cell growth. In the AG senescence accelerated by the mutation of Diuretic hormone 31 (Dh31), the presence of InR::mCherry puncta was more pronounced compared to the wild type. These findings raise expectations for the utility of the newly created InR::mCherry/EYFP alleles for studying the precise expression levels and subcellular localization of InR. Furthermore, this fluorescently tagged allele approach can be extended to investigate other membrane receptors with low abundance, facilitating the direct examination of their true expression and localization.

Dominance and multi-locus interaction

Dominance is usually considered a constant value that describes the relative difference in fitness or phenotype between heterozygotes and the average of homozygotes at a focal polymorphic locus. However, the observed dominance can vary with the genetic background of the focal locus. Here, alleles at other loci modify the observed phenotype through position effects or dominance modifiers that are sometimes associated with pathogen resistance, lineage, sex, or mating type. Theoretical models have illustrated how variable dominance appears in the context of multi-locus interaction (epistasis). Here, we review empirical evidence for variable dominance and how the observed patterns may be captured by proposed epistatic models. We highlight how integrating epistasis and dominance is crucial for comprehensively understanding adaptation and speciation.

The novel HLA-B*44:48:02 allele characterised by two different sequencing-based typing techniques

The novel allele HLA-B*44:48:02 differs from HLA-B*44:48:01 by one synonymous nucleotide substitution in exon 3.

Family-based association of HLA-DRB1 and DQB1 alleles and haplotypes in a group of Iranian Type 1 diabetes children

This family-based study was conducted in a group of Iranians with Type 1 diabetes (T1D) to investigate the transmission from parents of risk and non-risk HLA alleles and haplotypes, and to estimate the genetic risk score for this disease within this population. A total of 240 T1D subjects including 111 parent-child trios (111 children with T1D, 133 siblings, and 222 parents) and 330 ethnically matched healthy individuals were recruited. High-resolution HLA typing for DRB1/DQB1 loci was performed for all study subjects (n = 925) using polymerase chain reaction-sequence-specific oligonucleotide probe method. The highest predisposing effect on developing T1D was conferred by the following haplotypes both in all subjects and in probands compared to controls: DRB1*04:05-DQB1*03:02 (Pc = 2.97e-06 and Pc = 6.04e-10, respectively), DRB1*04:02-DQB1*03:02 (Pc = 5.94e-17 and Pc = 3.86e-09, respectively), and DRB1*03:01-DQB1*02:01 (Pc = 8.26e-29 and Pc = 6.56e-16, respectively). Conversely, the major protective haplotypes included DRB1*13:01-DQB1*06:03 (Pc = 6.99e-08), DRB1*15:01-DQB1*06:02 (Pc = 2.97e-06) in the cases versus controls. Also, DRB1*03:01-DQB1*02:01/DRB1*04:02|05-DQB1*03:02 and DRB1*03:01-DQB1*02:01/DRB1*03:01-DQB1*02:01 diplotypes conferred the highest predisposing effect in the cases (Pc = 8.65e-17 and Pc = 6.26e-08, respectively) and in probands (Pc = 5.4e-15 and Pc = 0.001, respectively) compared to controls. Transmission disequilibrium test showed that the highest risk was conferred by DRB1*04:02-DQB1*03:02 (Pc = 3.26e-05) and DRB1*03:01-DQB1*02:01 (Pc = 1.78e-12) haplotypes and the highest protection by DRB1*14:01-DQB1*05:03 (Pc = 8.66e-05), DRB1*15:01-DQB1*06:02 (Pc = 0.002), and DRB1*11:01-DQB1*03:01 (Pc = 0.0003) haplotypes. Based on logistic regression analysis, carriage of risk haplotypes increased the risk of T1D development 24.5 times in the Iranian population (p = 5.61e-13). Also, receiver operating characteristic curve analysis revealed a high predictive power of those risk haplotypes in discrimination of susceptible from healthy individuals (area under curve: 0.88, p = 5.5e-32). Our study highlights the potential utility of genetic risk assessment based on HLA diplotypes for predicting T1D risk in individuals, particularly among family members of affected children in our population.

The nature and distribution of putative non-functional alleles suggest only two independent events at the origins of Astyanax mexicanus cavefish populations

BACKGROUND

Several studies suggested that cavefish populations of Astyanax mexicanus settled during the Late Pleistocene. This implies that the cavefish's most conspicuous phenotypic changes, blindness and depigmentation, and more cryptic characters important for cave life, evolved rapidly.

RESULTS

Using the published genomes of 47 Astyanax cavefish from la Cueva de El Pachón, El Sótano de la Tinaja, La Cueva Chica and El Sótano de Molino, we searched for putative loss-of-function mutations in previously defined sets of genes, i.e., vision, circadian clock and pigmentation genes. Putative non-functional alleles for four vision genes were identified. Then, we searched genome-wide for putative non-functional alleles in these four cave populations. Among 512 genes with segregating putative non-functional alleles in cavefish that are absent in surface fish, we found an enrichment in visual perception genes. Among cavefish populations, different levels of shared putative non-functional alleles were found. Using a subset of 12 genes for which putative loss-of-function mutations were found, we extend the analysis of shared pseudogenes to 11 cave populations. Using a subset of six genes for which putative loss-of-function mutations were found in the El Sótano del Toro population, where extensive hybridization with surface fish occurs, we found a correlation between the level of eye regression and the amount of putative non-functional alleles.

CONCLUSIONS

We confirm that very few putative non-functional alleles are present in a large set of vision genes, in accordance with the recent origin of Astyanax mexicanus cavefish. Furthermore, the genome-wide analysis indicates an enrichment of putative loss-of-function alleles in genes with vision-related GO-terms, suggesting that visual perception may be the function chiefly impacted by gene losses related to the shift from a surface to a cave environment. The geographic distribution of putative loss-of-function alleles newly suggests that cave populations from Sierra de Guatemala and Sierra de El Abra share a common origin, albeit followed by independent evolution for a long period. It also supports that populations from the Micos area have an independent origin. In El Sótano del Toro, the troglomorphic phenotype is maintained despite massive introgression of the surface genome.

GEMIN4 Variants: Risk Profiling, Bioinformatics, and Dynamic Simulations Uncover Susceptibility to Bladder Carcinoma

BACKGROUND

The relationship between GEMIN4 genetic variants and cancer, especially bladder carcinoma (BLCA), has been explored without conclusive results. This study aims to elucidate the link between GEMIN4 polymorphisms and BLCA susceptibility through genetic analyses, bioinformatics, and molecular dynamics (MD) simulations.

METHODS

A cohort of 249 participants (121 BLCA patients and 128 unrelated controls) was enrolled. PCR was employed for allelic discrimination of GEMIN4 variants, followed by subgroup stratification, haplotype analyses, structural prediction using the AlphaFold2 prediction tool, subsequent MD simulations, structural analysis, and residue interaction mapping using Desmond, UCSF ChimeraX, and Cytoscape softwares.

RESULTS

The rs.2740348*G variant demonstrated a protective role against BLCA in allelic (OR = 0.55, p = 0.002) and recessive (OR = 0.54, p = 0.017) models, whereas the rs.7813*T variant increased BLCA risk under the recessive model (OR = 1.90, p = 0.019). Haplotype analysis revealed a significant association between GEMIN4 haplotype (rs.2740348*C/rs.7813*T) with increased BLCA risk (OR = 2.01, p = 0.004). Univariate analysis revealed associations of the variants with albumin levels and absolute neutrophil count in BLCA patients. Pathogenicity evaluation categorized p.Gln450Glu as neutral and p.Arg1033Cys as deleterious. MD simulations revealed structural alterations and conformational shifts in the GEMIN4 protein induced by the Glu450 and Cys1033 mutations.

CONCLUSIONS

The study highlights the dual role of GEMIN4 variants in BLCA susceptibility, with rs.2740348 conferring protection and rs.7813 increasing risk. The Glu450 residue positively impacted protein stability, while Cys1033 had a detrimental effect on protein function. These findings underscore the significance of GEMIN4 variants in BLCA susceptibility and pave the way for future diagnostic and therapeutic initiatives.

Genotyping of Toxoplasma gondii in domestic animals from Campeche, México, reveals virulent genotypes and a recombinant ROP5 allele

Toxoplasma gondii has at least 318 genotypes distributed worldwide, and tropical regions usually have greater genetic diversity. Campeche is a state located in the southeastern region of México and has favourable climate conditions for the replication and dissemination of this protozoan, similar to those in South American countries where broad genetic diversity has been described. Thus, in this study, 4 T. gondii isolates were obtained from tissues of stray dogs and free-range chickens in Campeche, México, and were genotyped by Mn-PCR-RFLP with 10 typing markers (SAG1, altSAG2, SAG3, BTUB, GRA6, c22-8, c29-2, L358, PK1 and Apico) and 5 virulence markers (CS3, ROP16, ROP17, ROP18 and ROP5) to provide new information about the distribution and virulence prediction of T. gondii genotypes. Two isolates of T. gondii genotype #116 and 2 of genotype #38 were obtained from stray dogs and chickens, respectively. The parasite load found in these species was between <50 and more than 35 000 tachyzoites per mg of tissue. Virulence marker genotyping revealed a recombinant 1&3 ROP5 RFLP pattern in 2 ToxoDB #116 isolates with no prediction of virulence in a murine model, while in the 2 ToxoDB #38 isolates, the ROP18/ROP5 combination predicted high virulence. Considering all the typed markers, there is a predominance of type I and III alleles, as constantly reported for the isolates characterized in various regions of México. It is crucial to determine their phenotype to corroborate the genetic virulence profile of the T. gondii isolates obtained in this study.

An AluYa5 Insertion in the 3'UTR of COL4A1 and Cerebral Small Vessel Disease

Importance

Cerebral small vessel diseases (CSVDs) account for one-fifth of stroke cases. Numerous familial cases remain unresolved after routine screening of known CSVD genes.

Objective

To identify novel genes and mechanisms associated with familial CSVD.

Design, Setting, and Participants

This 2-stage study involved linkage analysis and a case-control study; linkage analysis and whole exome and genome sequencing were used to identify candidate gene variants in 2 large families with CSVD (9 patients with CSVD). Then, a case-control analysis was conducted on 246 unrelated probands, including probands from these 2 families and 244 additional probands. All probands (clinical onset

Main Outcomes and Measures

A pathogenic AluYa5 insertion was identified within the COL4A1 3'UTR in the 2 large families with CSVD. Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR), Western blot, and long-read RNA sequencing were used to investigate outcomes associated with the insertion using patient fibroblasts. Clinical and magnetic resonance imaging features of probands with variants and available relatives were assessed.

Results

Among 246 probands (141 females [57.3%]; median [IQR] age at referral, 56 [49-64] years), 7 patients of French ancestry carried the insertion. This insertion was absent in 467 healthy French individuals in a control group (odds ratio, ∞; 95% CI, 2.78 to ∞; P = 5 × 10-4) and 10 847 individuals from the gnomAD structural variant database (odds ratio, ∞; 95% CI, 64.77 to ∞; P = 2.42 × 10-12). In these 7 patients' families, 19 family members with CSVD carried the insertion. RT-qPCR and Western blot showed an upregulation of COL4A1 mRNA (10.6-fold increase; 95% CI, 1.4-fold to 17.1-fold increase) and protein levels (2.8-fold increase; 95% CI, 2.1-fold to 3.5-fold increase) in patient vs control group fibroblasts. Long-read RNA sequencing data showed that the insertion was associated with perturbation in the use of canonical COL4A1 polyadenylation signals (approximately 87% of isoforms transcribed from the wild type allele vs 5% of isoforms transcribed from the allele with the insertion used the 2 distal canonical polyadenylation signals). The main clinical feature of individuals with CSVD was the recurrence of pontine ischemic lesions starting at an early age (17 of 19 patients [89.5%]).

Conclusions and relevance

This study found a novel mechanism associated with COL4A1 upregulation and a highly penetrant adult-onset CSVD. These findings suggest that quantitative alterations of the cerebrovascular matrisome are associated with CSVD pathogenesis, with diagnostic and therapeutic implications.

Collapse

Key Words Collapse

MESH Headings

• Adult
• Female
• Humans
• Middle Aged
• 3' Untranslated Regions/genetics
• Alleles
• Case-Control Studies
• Cerebral Small Vessel Diseases/genetics
• Collagen Type IV/metabolism
• Protein Isoforms
• Mutagenesis, Insertional

Collapse

Grants Collapse

Affiliation(s)

Chaker Aloui NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France
Lisa Neumann NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France
Françoise Bergametti NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France
Eric Sartori Service de Neurologie, Centre Hospitalier Bretagne Sud, Lorient, France
Marc Herbreteau Service de Neurologie, Centre Hospitalier Bretagne Sud, Lorient, France
Arnaud Maillard Assistance Publique-Hôpitaux de Paris, Service de Génétique Moléculaire Neurovasculaire, Hôpital Saint-Louis, Paris, France
Thibault Coste NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France Assistance Publique-Hôpitaux de Paris, Service de Génétique Moléculaire Neurovasculaire, Hôpital Saint-Louis, Paris, France
Hélène Morel NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France Assistance Publique-Hôpitaux de Paris, Service de Génétique Moléculaire Neurovasculaire, Hôpital Saint-Louis, Paris, France
Dominique Hervé NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France Assistance Publique-Hôpitaux de Paris, Service de Neurologie, Hôpital Lariboisière, Paris, France
Hugues Chabriat NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France Assistance Publique-Hôpitaux de Paris, Service de Neurologie, Hôpital Lariboisière, Paris, France
Serge Timsit Service de Neurologie Vasculaire, Centre Hospitalier Régional Universitaire de Brest, Brest, France
Irina Viakhireva Service de Neurologie Vasculaire, Centre Hospitalier Régional Universitaire de Brest, Brest, France
Yves Denoyer Service de Neurologie, Centre Hospitalier Bretagne Sud, Lorient, France Université de Rennes, Laboratoire Traitement du Signal et de l'Image, Institut National de la Santé Et de la Recherche Médicale Unité Mixte de Recherche 1099, Rennes, France
Rémi Allibert Service de Neurologie, Unité Neurovasculaire, Centre Hospitalier Universitaire de Saint Etienne, Saint Etienne, France
Florence Demurger Service de Neurologie, Unité Neurovasculaire, Centre Hospitalier Bretagne Atlantique, Vannes, France
Cedric Gollion Service de Neurologie, Centre Hospitalier Universitaire de Toulouse, Toulouse, France
Patrick Vermersch Univ. Lille, Institut National de la Santé Et de la Recherche Médicale Unité Mixte de Recherche 1172 LilNCog, Centre Hospitalier Universitaire Lille, Fédérations Hospitalo-Universitaire Precise, Lille, France
Florence Marchelli Assistance Publique-Hôpitaux de Paris, Service de Génétique Moléculaire Neurovasculaire, Hôpital Saint-Louis, Paris, France
Corinne Blugeon GenomiqueENS, Institut de Biologie de l’Ecole Normale Supérieur, Département de biologie, École Normale Supérieure, Centre National de la Recherche Scientifique, Institut National de la Santé Et de la Recherche Médicale, Université Paris Sciences et Lettres, Paris, France
Sophie Lemoine GenomiqueENS, Institut de Biologie de l’Ecole Normale Supérieur, Département de biologie, École Normale Supérieure, Centre National de la Recherche Scientifique, Institut National de la Santé Et de la Recherche Médicale, Université Paris Sciences et Lettres, Paris, France
Claire Tourtier-Bellosta Service de Neurologie, Centre Hospitalier Bretagne Sud, Lorient, France
Alexis Brouazin Service de neurologie, Centre Hospitalier de Cornouaille, Quimper, France
Anne-Louise Leutenegger NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France
Eva Pipiras NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France Assistance Publique-Hôpitaux de Paris, Hôpitaux Jean Verdier et Armand Trousseau, Université Sorbonne Paris Nord, Bobigny, France
Elisabeth Tournier-Lasserve NeuroDiderot, Université Paris Cité, Institut National de la Santé Et de la Recherche Médicale, Unité Mixte de Recherche 1141, Paris, France Assistance Publique-Hôpitaux de Paris, Service de Génétique Moléculaire Neurovasculaire, Hôpital Saint-Louis, Paris, France

Collapse

Seq2Sat and SatAnalyzer toolkit: Towards comprehensive microsatellite genotyping from sequencing data

Accurate and efficient microsatellite loci genotyping is an essential process in population genetics that is also used in various demographic analyses. Protocols for next-generation sequencing of microsatellite loci enable high-throughput and cross-compatible allele scoring, common issues that are not addressed by conventional capillary-based approaches. To improve this process, we have developed an all-in-one software, called Seq2Sat (sequence to microsatellite), in C++ to support automated microsatellite genotyping. It directly takes raw reads of microsatellite amplicons and conducts read quality control before inferring genotypes based on depth-of-read, read ratio, sequence composition and length. We have also developed a module for sex identification based on sex chromosome-specific locus amplicons. To allow for greater user access and complement autoscoring, we developed SatAnalyzer (microsatellite analyzer), a user-friendly web-based platform that conducts reads-to-report analyses by calling Seq2Sat for genotype autoscoring and produces interactive genotype graphs for manual editing. SatAnalyzer also allows users to troubleshoot multiplex optimization by analysing read quality and distribution across loci and samples in support of high-quality library preparation. To evaluate its performance, we benchmarked our toolkit Seq2Sat/SatAnalyzer against a conventional capillary gel method and existing microsatellite genotyping software, MEGASAT, using two datasets. Results showed that SatAnalyzer can achieve >99.70% genotyping accuracy and Seq2Sat is ~5 times faster than MEGASAT despite many more informative tables and figures being generated. Seq2Sat and SatAnalyzer are freely available on github (https://github.com/ecogenomicscanada/Seq2Sat) and dockerhub (https://hub.docker.com/r/rocpengliu/satanalyzer).

The novel HLA-C*03:94:02 allele identified by next-generation sequencing in a Chinese individual

HLA-C*03:94:02 differs from HLA-C*03:94:01 by a single nucleotide substitution in exon 2 (codon 17 GGA->GGG).

Allele-specific antisense oligonucleotides for the treatment of BEST1-related dominantly inherited retinal diseases: An in vitro model

Retinal dystrophies are a common health problem worldwide that are currently incurable due to the inability of retinal cells to regenerate. Inherited retinal diseases (IRDs) are a diverse group of disorders characterized by progressive vision loss caused by photoreceptor cell dysfunction. The eye has always been an attractive organ for the development of novel therapies due to its independent access to the systemic pathway. Moreover, anti-sense oligonucleotides (ASOs), which facilitate manipulation of unwanted mRNAs via degradation or splicing, are undergoing rapid development and have been clinically deployed for the treatment of several diseases. The primary aim of this study was to establish a reliable in vitro model utilizing induced photoreceptor-like cells (PRCs) for assessing the efficacy and safety of ASOs targeting the BEST1 gene. Despite advances in gene therapy, effective treatments for a broad range of IRDs remain limited. An additional aim was to develop an in vitro model for evaluating RNA-based therapeutics, specifically ASOs, for the treatment in IRDs. Firstly, a cell culture model was established by induction of PRCs from dermal fibroblasts via direct programming. The induced PRCs were characterized at both the transcriptomic and protein level. Then, a common single nucleotide polymorphism (SNP) was identified in the BEST1 gene (rs1800007) for targeting with ASOs. ASOs were designed using the GapmeR strategy to target multiple alleles of this SNP, which is potentially suitable for a large proportion of the population. The efficacy and possible off-target effects of these ASOs were also analyzed in the induced PRC model. The findings show that the selected ASOs achieved allele-specific mRNA degradation with virtually no off-target effects on the global transcriptome profile, indicating their potential as safe and effective therapeutic agents. The presented in vitro model is a valuable platform for testing personalized IRD treatments and should inspire further research on RNA-based therapeutics. To the best of our knowledge this study is the first to test RNA-based therapeutics involving the use of ASOs in an induced PRC model. Based on the present findings, it will be possible to establish an ex vivo disease model using dermal fibroblast samples from affected individuals. In other words, the disease model and the ASOs that were successfully designed in this study can serve as a useful platform for the testing of personalized treatments for IRDs.

Prognostic impact of HLA supertype mismatch in single-unit cord blood transplantation

The "human leukocyte antigen (HLA) supertype" is a functional classification of HLA alleles, which was defined by structural features and peptide specificities, and has been reportedly associated with the clinical outcomes of viral infections and autoimmune diseases. Although the disparity in each HLA locus was reported to have no clinical significance in single-unit cord blood transplantation (sCBT), the clinical significance of the HLA supertype in sCBT remains unknown. Therefore, we retrospectively analyzed clinical data of 1603 patients who received sCBT in eight institutes in Japan between 2000 and 2017. Each HLA allele was categorized into 19 supertypes, and the prognostic effect of disparities was then assessed. An HLA-B supertype mismatch was identified as a poor prognostic factor (PFS: hazard ratio [HR] = 1.23, p = 0.00044) and was associated with a higher cumulative incidence (CI) of relapse (HR = 1.24, p = 0.013). However, an HLA-B supertype mismatch was not associated with the CI of acute and chronic graft-versus-host-disease. The multivariate analysis for relapse and PFS showed the significance of an HLA-B supertype mismatch independent of allelic mismatches, and other previously reported prognostic factors. HLA-B supertype-matched grafts should be selected in sCBT.

Spontaneous expression of the CIC::DUX4 fusion oncoprotein from a conditional allele potently drives sarcoma formation in genetically engineered mice

CIC::DUX4 sarcoma (CDS) is a rare but highly aggressive undifferentiated small round cell sarcoma driven by a fusion between the tumor suppressor Capicua (CIC) and DUX4. Currently, there are no effective treatments and efforts to identify and translate better therapies are limited by the scarcity of patient tumor samples and cell lines. To address this limitation, we generated three genetically engineered mouse models of CDS (Ch7CDS, Ai9CDS, and TOPCDS). Remarkably, chimeric mice from all three conditional models developed spontaneous soft tissue tumors and disseminated disease in the absence of Cre-recombinase. The penetrance of spontaneous (Cre-independent) tumor formation was complete irrespective of bi-allelic Cic function and the distance between adjacent loxP sites. Characterization of soft tissue and presumed metastatic tumors showed that they consistently expressed the CIC::DUX4 fusion protein and many downstream markers of the disease credentialing the models as CDS. In addition, tumor-derived cell lines were generated and ChIP-seq was preformed to map fusion-gene specific binding using an N-terminal HA epitope tag. These datasets, along with paired H3K27ac ChIP-sequencing maps, validate CIC::DUX4 as a neomorphic transcriptional activator. Moreover, they are consistent with a model where ETS family transcription factors are cooperative and redundant drivers of the core regulatory circuitry in CDS.

Canine coat color E locus updates: Identification of a new MC1R variant causing 'sable' coat color in English Cocker Spaniels and a proposed update to the E locus dominance hierarchy

The coat color phenotype 'sable' occurs in the English Cocker Spaniel dog breed. It closely resembles other canine color patterns known as domino/grizzle/pied (eA allele) and grizzle/domino (eG allele) determined by variants in the melanocortin 1 receptor gene (MC1R; 'extension' or E locus), a key multi-allele regulator of coat color. We examined genetic variation in MC1R, and found one new non-synonymous variant, c.250G>A (p.(Asp84Asn)), consistently associated with the English Cocker Spaniel 'sable' phenotype. We propose calling this newly identified allele eH and further show that the eA , eH and eG (previously known as EG ) alleles associate with similar phenotypes in dogs impacting genotypes regulated by beta-defensin 103 gene (CBD103; K locus) and agouti signaling protein gene (ASIP; A locus) in the absence of the EM and E alleles. This suggests that all three alleles are putative reduced-function variants of the MC1R gene. We propose the revised and updated E locus dominance hierarchy to be EM  > E > eA /eH /eG  > e1-3 .

Association of leptin receptor polymorphisms with susceptibility of non-small cell lung cancer: Evidence from 2249 subjects

Non-small cell lung cancer (NSCLC) is increasing dramatically. It is believed that energy metabolism-related genes could play an important role in etiology of NSCLC. In this study, we sought to assess the correlation between three LEPR single nucleotide polymorphisms (rs1137101, rs1137100 and rs6588147) with NSCLS susceptibility. In total, 1193 NSCLC cases and 1056 controls were included. SNPscan™ genotyping method was used to analyze the genotypes of LEPR polymorphisms. Compared to rs6588147 GG in LEPR gene, this study identified a protective role of LEPR rs6588147 GA and GA/AA for the occurrence of NSCLC (GA vs. GG [p = 0.021] and GA/AA vs. GG [p = 0.030]). As well, we found that a protective role of LEPR rs6588147 for the occurrence of non-SCC subgroup (p < 0.05). By logistic regression analysis, we found that the rs6588147 A allele related genotypes might play a protective role for the occurrence of NSCLC in drinking, BMI ≥24 kg/m2, smoking and male subgroups. We also found that the rs1137101 A allele related genotypes played a protective role for the occurrence of NSCLC in male, younger participants (under 59 years) and overweight/obesity (BMI ≥24 kg/m2) subgroups. We found that LEPR Ars1037100Ars1037101Ars6588147 haplotype might play a protective role for the occurrence of NSCLC (p = 0.013). In addition, our findings indicated that LEPR rs1137100 G>A SNP might increase the risk of lymph node metastases (p = 0.038). This study highlights that LEPR rs6588147, rs1137101 genotypes and LEPR Ars1037100Ars1037101Ars6588147 haplotype are correlated with the occurrence of NSCLC. LEPR rs1137100 G>A SNP increases the risk of lymph node metastases.

RNA Interference Effectors Selectively Silence the Pathogenic Variant GNAO1 c.607 G > A In Vitro

RNA interference (RNAi)-based therapeutics hold the potential for dominant genetic disorders, enabling sequence-specific inhibition of pathogenic gene products. We aimed to direct RNAi for the selective suppression of the heterozygous GNAO1 c.607 G > A variant causing GNAO1 encephalopathy. By screening short interfering RNA (siRNA), we showed that GNAO1 c.607G>A is a druggable target for RNAi. The si1488 candidate achieved at least twofold allelic discrimination and downregulated mutant protein to 35%. We created vectorized RNAi by incorporating the si1488 sequence into the short hairpin RNA (shRNA) in the adeno-associated virus (AAV) vector. The shRNA stem and loop were modified to improve the transcription, processing, and guide strand selection. All tested shRNA constructs demonstrated selectivity toward mutant GNAO1, while tweaking hairpin structure only marginally affected the silencing efficiency. The selectivity of shRNA-mediated silencing was confirmed in the context of AAV vector transduction. To conclude, RNAi effectors ranging from siRNA to AAV-RNAi achieve suppression of the pathogenic GNAO1 c.607G>A and discriminate alleles by the single-nucleotide substitution. For gene therapy development, it is crucial to demonstrate the benefit of these RNAi effectors in patient-specific neurons and animal models of the GNAO1 encephalopathy.

penA profile of Neisseria gonorrhoeae in Guangdong, China: Novel penA alleles are related to decreased susceptibility to ceftriaxone or cefixime

BACKGROUND

Resistance to extended-spectrum cephalosporins (ESCs) has become a public health concern with the spread of Neisseria gonorrhoeae and increasing antimicrobial resistance. Mutation of penA, encoding penicillin-binding protein 2, represents a mechanism of ESC resistance. This study sought to assess penA alleles and mutations associated with decreased susceptibility (DS) to ESCs in N. gonorrhoeae.

MATERIALS AND METHODS

In 2021, 347 gonococci were collected in Guangdong, China. Minimum inhibitory concentations (MICs) of ceftriaxone and cefixime were determined, and whole-genome sequencing and phylogenetic analysis were performed. Multi-locus sequence typing (MLST) and conventional resistance determinants such as penA, mtrR, PonA and PorB were analysed. penA was genotyped and sequence-aligned using PubMLST.

RESULTS

Genome-wide phylogenetic analysis revealed that the prevalence of DS to ESCs was highest in Clade 11.1 (100.0%), Clade 2 (66.7%) and Clade 0 (55.7%), and the leading cause was strains with penA-60.001 or new penA alleles in clades. The penA phylogenetic tree is divided into two branches: non-mosaic penA and mosaic penA. The latter contained penA-60.001, penA-10 and penA-34. penA profile analysis indicated that A311V and T483S are closely related to DS to ESCs in mosaic penA. The new alleles NEIS1753_2840 and NEIS1753_2837 are closely related to penA-60.001, with DS to ceftriaxone and cefixime of 100%. NEIS1753_2660, a derivative of penA-10 (A486V), has increased DS to ceftriaxone. NEIS1753_2846, a derivative of penA-34.007 (G546S), has increased DS to cefixime.

CONCLUSION

This study identified critical penA alleles related to elevated MICs, and trends of gonococcus-evolved mutated penA associated with DS to ESCs in Guangdong.

The novel HLA-A*30:01:23 allele characterised by two different sequencing-based typing techniques

The novel allele HLA-A*30:01:23 differs from HLA-A*30:01:01:01 by one synonymous nucleotide substitution in exon 2.

The novel HLA-DQA1*05:05:17:03 allele, identified in a potential organ donor

HLA-DQA1*05:05:17:03 differs from HLA-DQA1*05:05:01:02 by a single base substitution in exon 1 and HLA-DQA1*05:05:17:01 within introns 1 and 2.