Development of multiplex allele-specific RT-qPCR assays for differentiation of SARS-CoV-2 Omicron subvariants

Quick differentiation of current circulating variants and the emerging recombinant variants of SARS-CoV-2 is essential to monitor their transmissions. However, the widely applied gene sequencing method is time-consuming and costly especially when facing recombinant variants, because a large part or whole genome sequencing is required. Allele-specific reverse transcriptase real time RT-PCR (RT-qPCR) represents a quick and cost-effective method for SNP (single nucleotide polymorphism) genotyping and has been successfully applied for SARS-CoV-2 variant screening. In the present study, we developed a panel of 5 multiplex allele-specific RT-qPCR assays targeting 20 key mutations for quick differentiation of the Omicron subvariants (BA.1 to BA.5 and their descendants) and the recombinant variants (XBB.1 and XBB.1.5). Two parallel multiplex RT-qPCR reactions were designed to separately target the prototype allele and the mutated allele of each mutation in the allele-specific RT-qPCR assay. Optimal annealing temperatures, primer and probe dosage, and time for annealing/extension for each reaction were determined by multi-factor and multi-level orthogonal test. The variation of Cp (crossing point) values (ΔCp) between the two multiplex RT-qPCR reactions was applied to determine if a mutation occurs or not. SARS-CoV-2 subvariants and related recombinant variants were differentiated by their unique mutation patterns. The developed multiplex allele-specific RT-qPCR assays exhibited excellent analytical sensitivities (with limits of detection (LoDs) of 1.47-18.52 copies per reaction), wide linear detection ranges (109-100 copies per reaction), good amplification efficiencies (88.25 to 110.68%), excellent reproducibility (coefficient of variations (CVs) < 5% in both intra-assay and inter-assay tests), and good clinical performances (99.5-100% consistencies with Sanger sequencing). The developed multiplex allele-specific RT-qPCR assays in the present study provide an alternative tool for quick differentiation of the SARS-CoV-2 Omicron subvariants and their recombinant variants. KEY POINTS: • A panel of five multiplex allele-specific RT-qPCR assays for quick differentiation of 11 SARS-CoV-2 Omicron subvariants (BA.1, BA.2, BA.4, BA.5, and their descendants) and 2 recombinant variants (XBB.1 and XBB.1.5). • The developed assays exhibited good analytical sensitivities and reproducibility, wide linear detection ranges, and good clinical performances, providing an alternative tool for quick differentiation of the SARS-CoV-2 Omicron subvariants and their recombinant variants.

Insights into the role of JAK2-I724T variant in myeloproliferative neoplasms from a unique cohort of New Zealand patients

OBJECTIVES

This study aimed to compile bioinformatic and experimental information for JAK2 missense variants previously reported in myeloproliferative neoplasms (MPN) and determine if germline JAK2-I724T, recently found to be common in New Zealand Polynesians, associates with MPN.

METHODS

For all JAK2 variants found in the literature, gnomAD_exome allele frequencies were extracted and REVEL scores were calculated using the dbNSFP database. We investigated the prevalence of JAK2-I724T in a cohort of 111 New Zealand MPN patients using a TaqMan assay, examined its allelic co-occurrence with JAK2-V617F using Oxford Nanopore sequencing, and modelled the impact of I724T on JAK2 using I-Mutant and ChimeraX software.

RESULTS

Several non-V617F JAK2 variants previously reported in MPN had REVEL scores greater than 0.5, suggesting pathogenicity. JAK2-I724T (REVEL score 0.753) was more common in New Zealand Polynesian MPN patients (n = 2/27; 7.4%) than in other New Zealand patients (n = 0/84; 0%) but less common than expected for healthy Polynesians (n = 56/377; 14.9%). Patients carrying I724T (n = 2), one with polycythaemia vera and one with essential thrombocythaemia, had high-risk MPN. Both patients with JAK2-I724T were also positive for JAK2-V617F, found on the same allele as I724T, as well as separately. In silico modelling did not identify noticeable structural changes that would give JAK2-I724T a gain-of-function.

CONCLUSION

Several non-canonical JAK2 variants with high REVEL scores have been reported in MPN, highlighting the need to further understand their relationship with disease. The JAK2-I724T variant does not drive MPN, but additional investigations are required to exclude any potential modulatory effect on the MPN phenotype.

Introgression of OsAP47 by marker-assisted selection enhanced resistance against southern rice black-streaked dwarf virus disease

Southern rice black-streaked dwarf virus disease (SRBSDVD) is the most destructive viral disease in rice. In order to breeding resistant cultivars, Insertion-Deletion (InDel) markers were developed linked to OsAP47, the first isolated major resistance gene against SRBSDVD. Marker-assisted selection (MAS) was conducted to introduce this gene into the commercial variety. A rice line carrying homozygous resistance allele of OsAP47 was selected and named Kanghei No. 201 (KH201). Evaluated by artificial inoculation, KH201 showed significantly higher resistance than the recurrent parent Suxiu No.867 (SX867). And no significant differences were detected for KH201 in the yield-related components, including spikelets per panicle (SPP), ripened grains per panicle (RGPP), 1000-grain weight (TGW) and panicles per square meter (PPSM), leading to stable theoretical yield. The results indicated that introgression of OsAP47 improved rice resistance and can avoid yield losses produced by SRBSDVD. KH201 was demonstrated as a resistance material that could be used in rice breeding.

Large-scale DNA sequencing identifies rare variants associated with Systemic Lupus Erythematosus susceptibility in known risk genes

The identification of rare genetic variants associated to Systemic Lupus Erythematosus (SLE) could also help to understand the pathogenic mechanisms at the basis of the disease. In this study we have analyzed a cohort of 200 Italian SLE patients in order to explore the rare protein-coding variants in five genes (TNFAIP3, STAT4, IL10, TRAF3IP2, and HCP5) already investigated for commons variants found associated in our previous studies. Genomic DNA of 200 SLE patients was sequenced by whole exome sequencing. The identified variants were filtered by frequency and evaluated by in silico predictions. Allelic association analysis was performed with standard Fisher's exact test. Introducing a cutoff at MAF < 0.01, a total of 19 rare variants were identified. Seven of these variants were ultra-rare (MAF < 0.001) and six were absent in the GnomAD database. For TNFAIP3 gene, the variant c.A1939C was observed in 4 SLE patients and it is located in a region enriched in phosphorylation sites and affects the predict affinity of specific kinases. In TRAF3IP2 gene, we observed 5 different rare variants, including the novel variant c.G410A, located in the region that mediates interaction with TRAF6, and therefore a possible risk factor for SLE development. In STAT4 gene, we identified 6 different rare variants. Among these, three missense variants decrease the stability of this protein. Moreover, 3 novel rare variants were detected in 3 SLE patients. In particular, c.A767T variant was predicted as damaging by six prediction tools. Concluding, we have observed that even in genes whose common variability is associated with SLE susceptibility, it is possible to identify rare variants that could have a strong effect in the disease development and could therefore allow a better understanding of the functional domain involved.

Monoallelic Loss-of-Function IFT140 Pathogenic Variants Cause Autosomal Dominant Polycystic Kidney Disease: A Confirmatory Study With Suspicion of an Additional Cardiac Phenotype

Autosomal dominant polycystic kidney disease (ADPKD) is the most common inherited kidney disease. While biallelic variants affecting IFT140 are responsible for Mainzer-Saldino syndrome (characterized by severe ciliopathy causing skeletal abnormalities, kidney disease, and cysts), monoallelic loss-of-function (LoF) variants have been recently reported as an important cause of ADPKD beyond PKD1/2 genes. Herein, we report 6 non-family-related cases of monoallelic IFT140 LoF variants, identified from 1,340 exomes sequenced for nephrological indications in our local database. Every patient presented with polycystic kidney disease. Furthermore, the mother of a boy diagnosed with Mainzer-Saldino syndrome with a biallelic variant affecting IFT140 presented with several bilateral cysts, revealed after kidney imaging, and was found to carry a pathologic frameshift IFT140 variation. As well as this particular Mainzer-Saldino case, our 6 additional patients confirm that heterozygous IFT140 frameshift variants are responsible for the cystic phenotype and kidney failure. Interestingly, of the 6 patients, 2 also exhibited dilated cardiomyopathy, which was of unknown origin, as no genetic cause was found after exome sequencing analysis, suggesting a potential connection between IFT140 and heart disease.

Replication Study of the Association of GAS6 and PROS1 Polymorphisms with Behçet's Disease in a Japanese Population

PURPOSE

To investigate whether polymorphisms of GAS6 and PROS1, which each encode protein ligands for a family of tyrosine kinase receptors, are associated with Behçet's disease (BD) in a Japanese population.

METHODS

We recruited 734 Japanese patients with BD and 1789 Japanese healthy controls. In all participants, we genotyped two single-nucleotide polymorphisms (SNPs) reportedly associated with BD: rs9577873 in GAS6 and rs4857037 in PROS1.

RESULTS

We found that GAS6 rs9577873 was not significantly associated with BD. In contrast, PROS1 rs4857037, specifically the A allele, was associated with increased risk for BD. The A allele was also significantly associated with BD under additive and recessive genetic models. Expression analysis revealed that this allele was significantly associated with increased PROS1 expression.

CONCLUSIONS

Our findings suggest that increased PROS1 expression related to the A risk allele of rs4857037 affects tyrosine kinase receptor signaling, contributing to the development of BD.

Genetic etiology of truncus arteriosus excluding 22q11.2 deletion syndrome and identification of c.1617del, a prevalent variant in TMEM260, in the Japanese population

Truncus Arteriosus (TA) is a congenital heart disease characterized by a single common blood vessel emerging from the right and left ventricles instead of the main pulmonary artery and aorta. TA accounts for 4% of all critical congenital heart diseases. The most common cause of TA is 22q11.2 deletion syndrome, accounting for 12-35% of all TA cases. However, no major causes of TA other than 22q11.2 deletion have been reported. We performed whole-genome sequencing of 11 Japanese patients having TA without 22q11.2 deletion. Among five patients, we identified pathogenic variants in TMEM260; the biallelic loss-of-function variants of which have recently been associated with structural heart defects and renal anomalies syndrome (SHDRA). In one patient, we identified a de novo pathogenic variant in GATA6, and in another patient, we identified a de novo probably pathogenic variant in NOTCH1. Notably, we identified a prevalent variant in TMEM260 (ENST00000261556.6), c.1617del (p.Trp539Cysfs*9), in 8/22 alleles among the 11 patients. The c.1617del variant was estimated to occur approximately 23 kiloyears ago. Based on the allele frequency of the c.1617del variant in the Japanese population (0.36%), approximately 26% of Japanese patients afflicted with TA could harbor homozygous c.1617del variants. This study highlights TMEM260, especially c.1617del, as a major genetic cause of TA in the Japanese population.

DNA recovery after sequential processing of latent fingerprints on black polyethylene plastic

Latent fingerprints on plastic substrates can be visualized by using sequential treatments to enhance the contrast between the fingerprint residues and underlying substrate; however, the extent to which these processes affect subsequent DNA analysis is mostly unknown. Latent fingerprints deposited on black plastic by one donor were visualized with single-process fingerprint powders (i.e., white powder, bichromatic powder, or bichromatic magnetic powder) or sequential treatments (i.e., laser → reflected ultraviolet imaging system (RUVIS) → CA fuming → RUVIS → Rhodamine 6G, Ardrox, and MBD (RAM) or CA fuming → RAM/laser → bichromatic magnetic powder). Samples were examined after the addition of each treatment. DNA was collected using cotton swabs, extracted, quantified, and amplified. DNA yields, peak heights, number of alleles obtained, and percentage of DNA profiles eligible for CODIS upload were examined. Latent fingerprints processed with the laser and up to three sequential treatments generated DNA profiles with significantly higher peaks heights than those of the untreated samples. Fingerprints processed with the laser and up to two sequential treatments generated DNA profiles with significantly more alleles. All methods beginning with laser enhancement generated more CODIS-eligible profiles. Additional research is needed to determine the extent to which initial laser enhancement impacts the success of downstream DNA profiling results. Although DNA profile development is not guaranteed due to the variable quantities of DNA contained within latent fingerprints, the selection of an appropriate latent fingerprint visualization method could maximize both fingerprint detection and the generation of CODIS-eligible DNA profiles.

Identifying the genetic associations among the psoriasis patients in eastern India

Psoriasis is a multifactorial genetic disorder manifested by hyperproliferation and abnormal differentiation of epidermal keratinocytes, along with the infiltration of inflammatory cells into the skin. Although ~80 genetic susceptibility variants were reported in psoriasis, many loci showed population-specific associations, warranting the need for more population-specific association studies in psoriasis. We determined the association of forty single nucleotide polymorphisms (SNPs) among 2136 psoriasis patients and normal individuals from eastern India. We investigated the expression of corresponding genes and evaluated the protein structure stability for the genes with susceptible coding variants. We found fifteen SNPs significantly associated with psoriasis, while additional three SNPs showed significant association when we classified the patients based on the presence of HLA-Cw6 allele. Epistatic interaction between HLA-Cw6 and other associated loci showed significant association with the SNPs at PSORS1 region, along with other five SNPs outside PSORS1. Three genes showed significant differential expression in psoriatic tissues compared to the adjacent normal skin tissues but were not differential when classified the patients based on their genotypes. SNP rs495337 at SPATA2 (Spermatogenesis Associated 2) showed a 1.2-fold increased risk among the HLA-Cw6 patients compared to combined samples. We found significant downregulation of SPATA2 among the patients with risk genotypes and HLA-Cw6 allele compared to the non-risk genotypes. Protein structure stability analysis showed reduced structural stability for all the mutant residues caused by the associated coding variants. Our study evaluated the genetic associations of psoriasis-susceptible variants in India and evaluated the possible functional significance of these associated variants in psoriasis.

HLA Class II (DR, DQ, DP) Genes Were Separately Associated With the Progression From Seroconversion to Onset of Type 1 Diabetes Among Participants in Two Diabetes Prevention Trials (DPT-1 and TN07)

OBJECTIVE

To explore associations of HLA class II genes (HLAII) with the progression of islet autoimmunity from asymptomatic to symptomatic type 1 diabetes (T1D).

RESEARCH DESIGN AND METHODS

Next-generation targeted sequencing was used to genotype eight HLAII genes (DQA1, DQB1, DRB1, DRB3, DRB4, DRB5, DPA1, DPB1) in 1,216 participants from the Diabetes Prevention Trial-1 and Randomized Diabetes Prevention Trial with Oral Insulin sponsored by TrialNet. By the linkage disequilibrium, DQA1 and DQB1 are haplotyped to form DQ haplotypes; DP and DR haplotypes are similarly constructed. Together with available clinical covariables, we applied the Cox regression model to assess HLAII immunogenic associations with the disease progression.

RESULTS

First, the current investigation updated the previously reported genetic associations of DQA1*03:01-DQB1*03:02 (hazard ratio [HR] = 1.25, P = 3.50*10-3) and DQA1*03:03-DQB1*03:01 (HR = 0.56, P = 1.16*10-3), and also uncovered a risk association with DQA1*05:01-DQB1*02:01 (HR = 1.19, P = 0.041). Second, after adjusting for DQ, DPA1*02:01-DPB1*11:01 and DPA1*01:03-DPB1*03:01 were found to have opposite associations with progression (HR = 1.98 and 0.70, P = 0.021 and 6.16*10-3, respectively). Third, DRB1*03:01-DRB3*01:01 and DRB1*03:01-DRB3*02:02, sharing the DRB1*03:01, had opposite associations (HR = 0.73 and 1.44, P = 0.04 and 0.019, respectively), indicating a role of DRB3. Meanwhile, DRB1*12:01-DRB3*02:02 and DRB1*01:03 alone were found to associate with progression (HR = 2.6 and 2.32, P = 0.018 and 0.039, respectively). Fourth, through enumerating all heterodimers, it was found that both DQ and DP could exhibit associations with disease progression.

CONCLUSIONS

These results suggest that HLAII polymorphisms influence progression from islet autoimmunity to T1D among at-risk subjects with islet autoantibodies.

Template-specific optimization of NGS genotyping pipelines reveals allele-specific variation in MHC gene expression

Using high-throughput sequencing for precise genotyping of multi-locus gene families, such as the major histocompatibility complex (MHC), remains challenging, due to the complexity of the data and difficulties in distinguishing genuine from erroneous variants. Several dedicated genotyping pipelines for data from high-throughput sequencing, such as next-generation sequencing (NGS), have been developed to tackle the ensuing risk of artificially inflated diversity. Here, we thoroughly assess three such multi-locus genotyping pipelines for NGS data, the DOC method, AmpliSAS and ACACIA, using MHC class IIβ data sets of three-spined stickleback gDNA, cDNA and "artificial" plasmid samples with known allelic diversity. We show that genotyping of gDNA and plasmid samples at optimal pipeline parameters was highly accurate and reproducible across methods. However, for cDNA data, the gDNA-optimal parameter configuration yielded decreased overall genotyping precision and consistency between pipelines. Further adjustments of key clustering parameters were required tο account for higher error rates and larger variation in sequencing depth per allele, highlighting the importance of template-specific pipeline optimization for reliable genotyping of multi-locus gene families. Through accurate paired gDNA-cDNA typing and MHC-II haplotype inference, we show that MHC-II allele-specific expression levels correlate negatively with allele number across haplotypes. Lastly, sibship-assisted cDNA-typing of MHC-I revealed novel variants linked in haplotype blocks, and a higher-than-previously-reported individual MHC-I allelic diversity. In conclusion, we provide novel genotyping protocols for the three-spined stickleback MHC-I and -II genes, and evaluate the performance of popular NGS-genotyping pipelines. We also show that fine-tuned genotyping of paired gDNA-cDNA samples facilitates amplification bias-corrected MHC allele expression analysis.

Uncoupling of Cytochrome P450 2B6 and stimulation of reactive oxygen species production in pharmacogenomic alleles affected by interethnic variability

Cytochrome P450 mediated substrate metabolism is generally characterized by the formation of reactive intermediates. In vitro and in vivo reaction uncoupling, results in the accumulation and dissociation of reactive intermediates, leading to increased ROS formation. The susceptibility towards uncoupling and altered metabolic activity is partly modulated by pharmacogenomic alleles resulting in amino acid substitutions. A large variability in the prevalence of these alleles has been demonstrated in CYP2B6, with some being predominantly unique to African populations. The aim of this study is to characterize the uncoupling potential of recombinant CYP2B6*1, CYP2B6*6 and CYP2B6*34 metabolism of specific substrates. Therefore, functional effects of these alterations on enzyme activity were determined by quantification of bupropion, efavirenz and ketamine biotransformation using HPLC-MS/MS. Determination of H2O2 levels was performed by the AmplexRed/horseradish peroxidase assay. Our studies of the amino acid substitutions Q172H, K262R and R487S revealed an exclusive use of the peroxide shunt for the metabolism of bupropion and ketamine by CYP2B6*K262R. Ketamine was also identified as a trigger for the peroxide shunt in CYP2B6*1 and all variants. Concurrently, ketamine acted as an uncoupler for all enzymes. We further showed that the expressed CYP2B6*34 allele results in the highest H2O2 formation. We therefore conclude that the reaction uncoupling and peroxide shunt are directly linked and can be substrate specifically induced with K262R carriers being most likely to use the peroxide shunt and R487S carrier being most prone to reaction uncoupling. This elucidates the functional diversity of pharmacogenomics in drug metabolism and safety.

Expression analysis of NF1-mutated alleles in a rare compound heterozygous spinal NF1 patient by digital PCR

BACKGROUD

Neurofibromatosis type 1 (NF1) is a heterogeneous neurocutaneous disorder. Spinal neurofibromatosis (SNF) is a distinct clinical entity of NF1, characterized by bilateral neurofibromas involving all spinal nerve roots. Although both forms are caused by intragenic heterozygous variants of NF1, missense variants have been associated with SNF, according to a dominant inheritance model causing haploinsufficiency. Most patients carry pathogenic variants in one of the NF1 alleles; nevertheless, patients with both NF1-mutated copies have been described. Interestingly, all NF1 variants carried by the known SNF compound heterozygotes were missense/splicing variants or in-frame insertion-deletions.

AIMS

To investigate whether there is a differential expression of NF1 variant alleles in an NF1 compound heterozygous SNF patient possibly contributing to clinical phenotype.

MATERIALS & METHODS

We performed an allele-specific expression study, by chip-based digital PCR, in an SNF family carrying two NF1 missense variants. We evaluated the expression levels of the two NF1-mutated alleles both carried by the compound heterozygous SNF patient and his relatives.

RESULTS

Both alleles were expressed at comparable levels in the patient and hyper-expressed compared to the wild-type alleles of healthy controls.

DISCUSSION

Here we provide new insights into expression studies of NF1-mutated transcripts suggesting that a novel pathogenetic mechanism, caused by gain-of-function variants, could be associated with SNF.

CONCLUSIONS

Further studies should be performed in larger cohorts, opening new perspectives in the NF1 pathogenesis comprehension.

An exploratory view into allelic drop-out of sequenced autosomal STRs

As massively parallel sequencing is implemented in forensic genetics, an understanding of sequence data must accompany these advancements, that is, accurate modeling of data for proper statistical analysis. Allelic drop-out, a common stochastic effect seen in genetic data, is often modeled in statistical analysis of STR results. This proof-of-concept study sequenced several serial dilutions of a standard sample ranging from 4 ng to 7.82 pg to evaluate allelic drop-out trends on a select panel of autosomal STRs using the ForenSeq™ DNA Signature Prep Kit, Primer Set A on the Illumina MiSeq FGx. Parameters assessed included locus, profile, and run specific information. A majority of the allelic drop-out occurred in DNA concentrations less than 31.25 pg. Statistical results indicated a need for locus-specific modeling based on STR descriptors, like simple versus compound repeat patterns. No correlation was seen between average read count of scored alleles and allelic drop-out at a locus. A statistical correlation was observed between the amount of allelic drop-out and the starting amount of DNA in a sample, average read count of a sample, and total read count generated on a flow cell. This study supports using common allelic drop-out factors used in fragment length analysis on sequenced STRs while including additional locus, sample, and run specific information. Results demonstrate multiple factors that can be considered when developing probability of allelic drop-out models for sequenced autosomal STRs including locus-specific analysis, total read count of a profile, and total read count sequenced on a flow cell.

Identification of Antibodies to DQβ:DRα Interisotypic Heterodimers in Human Sera

BACKGROUND

HLA class II antigens, DR, DQ, and DP, comprised an α and β chains, which typically combine, within the same isotype, to form the major histocompatibility complex:peptide complex. Interisotypic pairing is not commonly observed. Although reports of DQβ:DRα heterodimers exist, the pairing was reported to be unstable and, therefore, not studied to any extent.

METHODS

DQβ:DRα single antigens were produced through transfectant cell lines and used to identify and characterize positive reactive human sera by a multiplex bead-based assay.

RESULTS

Stable DQβ:DRα transfectants were constructed. Cell surface staining with class II-specific monoclonal antibodies revealed that some DQB1 alleles appear to be more efficient in expressing DQβ:DRα heterodimers. Interestingly, alleles within the same serological group varied in their efficiency of forming dimers on the cell surface. For example, DQβ0601:DRα had the highest transfection and cell membrane expression efficiency among 16 common DQB1 alleles tested. In contrast, DQβ0603:DRα-positive transfectants demonstrated minimal surface expression. Assembly of DQβ0601:DRα was not affected by the presence of a DQα chain. DQβ0601:DRα and DQβ0603:DRα single-antigen beads were used to screen human sera. Positive sera were identified that reacted to the unique epitopes of DQβ0601:DRα protein on the cell surface of the transfectants.

CONCLUSIONS

Our studies have demonstrated that unique DQβ:DRα heterodimers can be formed and are stably expressed on the cell surface. Such antigenic combinations, presented on single-antigen beads, demonstrated that patient sera can react with such heterodimers. Investigations on the potential clinical roles of antibodies against such interisotypic heterodimers are now possible.

Harmonizing the forensic nomenclature for STR loci D6S474 and DYS612

The autosomal STR D6S474 and the Y-chromosomal STR DYS612 have been reported in multiple ways in the forensic literature, with differences in both the bracketed repeat structures and counting of numerical length-based capillary electrophoresis (CE) alleles. These issues often come to light when STR loci are introduced in commercial assays and results compared with historical publications of allele frequency data, or multiple assays are characterized with reference materials. We review the forensic literature and other relevant information, and provide suggestions for the future treatment of each STR.

Allelic heterogeneity in a patient with postzygotic MTOR-related hypomelanosis of Ito with neurodevelopmental abnormalities

A case of mosaic MTOR-associated hemimegalencephaly and hypomelanosis of Ito, died at 33 probably because of sudden unexpected death in epilepsy. Assessment of the variant allele fraction (VAF) in different tissues postmortem showed high variability not correlated with clinical features, representing the most detailed assessment of VAFs in different tissues to date.

High number of candidate gene variants are identified as disease-causing in a period of 4 years

Advances in bioinformatic tools paired with the ongoing accumulation of genetic knowledge and periodic reanalysis of genomic sequencing data have led to an improvement in genetic diagnostic rates. Candidate gene variants (CGVs) identified during sequencing or on reanalysis but not yet implicated in human disease or associated with a phenotypically distinct condition are often not revisited, leading to missed diagnostic opportunities. Here, we revisited 33 such CGVs from our previously published study and determined that 16 of them are indeed disease-causing (novel or phenotype expansion) since their identification. These results emphasize the need to focus on previously identified CGVs during sequencing or reanalysis and the importance of sharing that information with researchers around the world, including relevant functional analysis to establish disease causality.

Investigation of the Differences in the Pharmacokinetics of CYP2D6 Substrates, Desipramine, and Dextromethorphan in Healthy African Subjects Carrying the Allelic Variants CYP2D6*17 and CYP2D6*29, When Compared with Normal Metabolizers

This study investigated the differences in the pharmacokinetics (PK) of dextromethorphan and desipramine in healthy African volunteers to understand the effect of allelic variants of the human cytochrome P450 2D6 (CYP2D6) enzyme, namely the diplotypes of CYP2D6*1/*2 (*1*1, *1*2, *2*2) and the genotypes of CYP2D6*17*17 and CYP2D6*29*29. Overall, 28 adults were included and split into 3 cohorts after genotype screening: CYP2D6*1/*2 (n = 12), CYP2D6*17*17 (n = 12), and CYP2D6*29*29 (n = 4). Each subject received a single oral dose of dextromethorphan 30 mg syrup on day 1 and desipramine 50 mg tablet on day 8. The PK parameters of area under the plasma concentration-time curve from time of dosing to time of last quantifiable concentration (AUClast), and extrapolated to infinity (AUCinf), and the maximum plasma concentration (Cmax) were determined. For both dextromethorphan and desipramine, AUCinf and Cmax were higher in subjects of the CYP2D6*29*29 and CYP2D6*17*17 cohorts, as compared with subjects in the CYP2D6*1/*2 diplotype cohort and with normal metabolizers from the literature. All PK parameters, including AUCinf, Cmax, and the elimination half-life, followed a similar trend: CYP2D6*17*17 > CYP2D6*29*29 > CYP2D6*1/*2. The plasma and urinary drug/metabolite exposure ratios of both drugs were higher in subjects of the CYP2D6*17*17 and CYP2D6*29*29 cohorts, when compared with subjects in the CYP2D6*1/*2 diplotype cohort. All adverse events were mild, except in 1 subject with CYP2D6*17*17 who had moderately severe headache with desipramine. These results indicate that subjects with CYP2D6*17*17 and CYP2D6*29*29 genotypes were 5-10 times slower metabolizers than those with CYP2D6*1/*2 diplotypes. These findings suggest that dose optimization may be required when administering CYP2D6 substrate drugs in African patients. Larger studies can further validate these findings.

Self- and non-self-DNA on hands and sleeve cuffs

Studying DNA transfer and persistence has become increasingly important over the last decade, due to the impressive sensitivity of modern DNA detection methods in forensic genetics. To improve our understanding of background DNA that could also potentially be transferred, we analyzed the DNA composition on the outside of sleeve cuffs and sampled DNA directly from the hands of four different collaborators upon their arrival at work during 25 working days. Sampling of their hands was repeated after several hours working in our department. The shedder status of the participants, as assumed from previous internal studies, was well re-produced in the study. However, we noticed that the DNA shedding capacity could also change drastically during the day, with one participant showing a more than sixfold increase between hands sampled in the morning and hands sampled in the afternoon. As expected, poor DNA shedders carry more relative amounts of non-self-DNA on their hands than good shedders. Non-self-alleles were detected in 95% of the samples. We also observed potential effects of hand washing and the mode of transport to get to work on the DNA amount. People living with family members occasionally carried their DNA on their hands and more frequently on their sleeve cuffs. Sleeve cuffs, as being close to our hands, have a large potential to transfer DNA from one place to another, yet they have sparsely been studied as DNA transfer intermediates so far. In general, we collected consistently more DNA from the sleeve cuffs than from the hands of the participants, demonstrating their importance as potential transfer vectors. More DNA was recovered from sleeve cuffs made of synthetic fabric than from cuffs made of cotton or leather. In the afternoon, DNA from co-habitant family members could not be detected on the hands anymore and the detection of profiles from colleagues became more frequent. From two out of 100 analyzed sleeve cuffs and two out of 200 sampled hands, we established unknown major DNA profiles that would have been suitable for an entry in the national DNA database. This finding demonstrates the possibility to transfer DNA that has most likely been picked up somewhere in the public space.

CALR-mutated patients with low allele burden represent a specific subtype of essential thrombocythemia: A study on behalf of FIM and GBMHM

A low allele burden (i.e., <20%) of the CALR driver mutation is found in 10.8% of CALR-mutated MPNs, mostly in essential thrombocythemia, and correlates with a milder phenotype and a more indolent evolution compared to patients with an allele burden ≥20%.

Estimates of heterozygosity from single nucleotide polymorphism markers are context-dependent and often wrong

Genetic diversity is frequently described using heterozygosity, particularly in a conservation context. Often, it is estimated using single nucleotide polymorphisms (SNPs); however, it has been shown that heterozygosity values calculated from SNPs can be biased by both study design and filtering parameters. Though solutions have been proposed to address these issues, our own work has found them to be inadequate in some circumstances. Here, we aimed to improve the reliability and comparability of heterozygosity estimates, specifically by investigating how sample size and missing data thresholds influenced the calculation of autosomal heterozygosity (heterozygosity calculated from across the genome, i.e. fixed and variable sites). We also explored how the standard practice of tri- and tetra-allelic site exclusion could bias heterozygosity estimates and influence eventual conclusions relating to genetic diversity. Across three distinct taxa (a frog, Litoria rubella; a tree, Eucalyptus microcarpa; and a grasshopper, Keyacris scurra), we found heterozygosity estimates to be meaningfully affected by sample size and missing data thresholds, partly due to the exclusion of tri- and tetra-allelic sites. These biases were inconsistent both between species and populations, with more diverse populations tending to have their estimates more severely affected, thus having potential to dramatically alter interpretations of genetic diversity. We propose a modified framework for calculating heterozygosity that reduces bias and improves the utility of heterozygosity as a measure of genetic diversity, whilst also highlighting the need for existing population genetic pipelines to be adjusted such that tri- and tetra-allelic sites be included in calculations.

Diverse INOSITOL PHOSPHORYLCERAMIDE SYNTHASE mutant alleles of Physcomitrium patens offer new insight into complex sphingolipid metabolism

Sphingolipids are widespread, abundant, and essential lipids in plants and in other eukaryotes. Glycosyl inositol phosphorylceramides (GIPCs) are the most abundant class of plant sphingolipids, and are enriched in the plasma membrane of plant cells. They have been difficult to study due to lethal or pleiotropic mutant phenotypes. To overcome this, we developed a CRISPR/Cas9-based method for generating multiple and varied knockdown and knockout populations of mutants in a given gene of interest in the model moss Physcomitrium patens. This system is uniquely convenient due to the predominantly haploid state of the Physcomitrium life cycle, and totipotency of Physcomitrium protoplasts used for transformation. We used this approach to target the INOSITOL PHOSPHORYLCERAMIDE SYNTHASE (IPCS) gene family, which catalyzes the first, committed step in the synthesis of GIPCs. We isolated knockout single mutants and knockdown higher-order mutants showing a spectrum of deficiencies in GIPC content. Remarkably, we also identified two mutant alleles accumulating inositol phosphorylceramides, the direct products of IPCS activity, and provide our best explanation for this unexpected phenotype. Our approach is broadly applicable for studying essential genes and gene families, and for obtaining unusual lesions within a gene of interest.

Copy number and selection of MHC genes in ruminants are related to habitat, average life span and diet

The ruminants, as the main group of livestock, have been extensively studied in terms of their physiology, endocrinology, biochemistry, genetics, and nutrition. Despite the wide geographic distribution and habitat diversity of animals in this group, their ecology and evolution remain poorly understood. In this study, we analyzed the gene copy number, selection, and ecological and evolutionary processes that have affected the evolution of major histocompatibility complex (MHC) genes across ruminant lineages based on available genomic data. The 51 species analyzed represented all six families of ruminants. Our finding indicated that the architecture of the MHC region is conserved in ruminants, but with variable copy numbers of MHC-I, MHC-IIA, and MHC-IIB genes. No lineage-specific gene duplication was observed in the MHC genes. The phylogenetic generalized least squares regression (PGLS) model revealed association between ecological and biological factors (habitat and lifespan) and gene duplication in DQA and DQB, but not in DRB. The selection pressure of DQA and DQB were related with lifespan, diet, and the ratio of genetic repeat elements. These results suggest that the MHC evolution in ruminants, including copy number and selection, has been influenced by genetic repeat elements, pathogen exposure risk, and intrinsic cost of possessing multiple MHC genes.

Exploring hematological alterations and genetics linked to SNV rs10974944 in myeloproliferative neoplasms among Amazon patients

BCR::ABL1-negative myeloproliferative neoplasms are hematopoietic disorders characterized by panmyelosis. JAK2 V617F is a frequent variant in these diseases and often occurs in the 46/1 haplotype. The G allele of rs10974944 has been shown to be associated with this variant, specifically its acquisition, correlations with familial cases, and laboratory alterations. This study evaluated the association between the 46/1 haplotype and JAK2 V617F in patients with myeloproliferative neoplasms in a population from the Brazilian Amazon. Clinical, laboratory and molecular sequencing analyses were considered. Carriers of the G allele of rs10974944 with polycythemia vera showed an increase in mean corpuscular volume and mean corpuscular hemoglobin, while in those with essential thrombocythemia, there was an elevation in red blood cells, hematocrit, and hemoglobin. Associations were observed between rs10974944 and the JAK2 V617F, in which the G allele (OR 3.4; p < 0.0001) and GG genotype (OR 4.9; p = 0.0016) were associated with JAK2 V617F + and an increase in variant allele frequency (GG: OR 15.8; p =  < 0.0001; G: OR 6.0; p = 0.0002). These results suggest an association between rs10974944 (G) and a status for JAK2 V617F, JAK2 V617F + _VAF ≥ 50%, and laboratory alterations in the erythroid lineage.

Reversal of C9orf72 mutation-induced transcriptional dysregulation and pathology in cultured human neurons by allele-specific excision

Efforts to genetically reverse C9orf72 pathology have been hampered by our incomplete understanding of the regulation of this complex locus. We generated five different genomic excisions at the C9orf72 locus in a patient-derived induced pluripotent stem cell (iPSC) line and a non-diseased wild-type (WT) line (11 total isogenic lines), and examined gene expression and pathological hallmarks of C9 frontotemporal dementia/amyotrophic lateral sclerosis in motor neurons differentiated from these lines. Comparing the excisions in these isogenic series removed the confounding effects of different genomic backgrounds and allowed us to probe the effects of specific genomic changes. A coding single nucleotide polymorphism in the patient cell line allowed us to distinguish transcripts from the normal vs. mutant allele. Using digital droplet PCR (ddPCR), we determined that transcription from the mutant allele is upregulated at least 10-fold, and that sense transcription is independently regulated from each allele. Surprisingly, excision of the WT allele increased pathologic dipeptide repeat poly-GP expression from the mutant allele. Importantly, a single allele was sufficient to supply a normal amount of protein, suggesting that the C9orf72 gene is haplo-sufficient in induced motor neurons. Excision of the mutant repeat expansion reverted all pathology (RNA abnormalities, dipeptide repeat production, and TDP-43 pathology) and improved electrophysiological function, whereas silencing sense expression did not eliminate all dipeptide repeat proteins, presumably because of the antisense expression. These data increase our understanding of C9orf72 gene regulation and inform gene therapy approaches, including antisense oligonucleotides (ASOs) and CRISPR gene editing.

Correlation analysis of serum TLR4 protein levels and TLR4 gene polymorphisms in gouty arthritis patients

OBJECTIVE

The Toll-like receptor (TLR) 4-mediated nuclear factor kappa B (NF-κB) signaling pathway regulates the production of inflammatory factors and plays a key role in the pathogenesis of gouty arthritis. The aim of the present study was to investigate the link among TLR4 gene polymorphisms at various loci, protein expression, and gouty arthritis susceptibility.

METHODS

Between 2016 and 2021, a case-control study was used to collect a total of 1207 study subjects, including 317 male patients with gouty arthritis (gout group) and 890 healthy males (control group). The association between gout susceptibility and different genetic models was analyzed by typing three loci of the TLR4 gene (rs2149356, rs2737191, and rs10759932) using a multiplex point mutation rapid assay, and the association between protein expression and gout was confirmed by measuring TLR4 protein concentrations using enzyme-linked immunosorbent assays (ELISAs).

RESULTS

In a codominant models AA and AG, the rs2737191 polymorphism in the gout group increased the risk of gout compared to the AA genotype (OR = 2.249, 95%CI 1.010~5.008), and the risk of gout was higher for those carrying the G allele compared to the A allele (OR = 2.227, 95%CI 1.006~4.932). TLR4 protein expression was different between the two groups with different locus genotypes. The differences in TLR4 protein expression between the gout group and control group were statistically significant between the following genotypes: the GG and GT genotypes of the rs2149356 polymorphism; the AA and AG genotypes of the rs2737191 polymorphism; and the TT and TC genotypes of the rs10759932 polymorphism(P<0.05). The TLR4 protein level in the gout group (19.19±3.09 ng/ml) was significantly higher than that in the control group (15.85±4.75 ng/ml).

CONCLUSION

The AG genotype of the TLR4 gene rs2737191 polymorphism may be correlated with the development of gouty arthritis. The level of TLR4 protein expression is significantly higher in patients with gouty arthritis than in controls, and there is a correlation between high TLR4 protein expression and the development of gouty arthritis.

Transcriptional priming and chromatin regulation during stochastic cell fate specification

Stochastic cell fate specification, in which a cell chooses between two or more fates with a set probability, diversifies cell subtypes in development. Although this is a vital process across species, a common mechanism for these cell fate decisions remains elusive. This review examines two well-characterized stochastic cell fate decisions to identify commonalities between their developmental programmes. In the fly eye, two subtypes of R7 photoreceptors are specified by the stochastic ON/OFF expression of a transcription factor, spineless. In the mouse olfactory system, olfactory sensory neurons (OSNs) randomly select to express one copy of an olfactory receptor (OR) gene out of a pool of 2800 alleles. Despite the differences in these sensory systems, both stochastic fate choices rely on the dynamic interplay between transcriptional priming, chromatin regulation and terminal gene expression. The coupling of transcription and chromatin modifications primes gene loci in undifferentiated neurons, enabling later expression during terminal differentiation. Here, we compare these mechanisms, examine broader implications for gene regulation during development and posit key challenges moving forward. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Impact of APOE-ε alleles on brain structure and cognitive function in healthy older adults: A VBM and DTI replication study

BACKGROUND

The Apolipoprotein E (APOE) gene has been established in the Alzheimer's disease (AD) literature to impact brain structure and function and may also show congruent effects in healthy older adults, although findings in this population are much less consistent. The current study aimed to replicate and expand the multimodal approach employed by Honea et al. Structural magnetic resonance imaging (MRI), diffusion tensor imaging (DTI), and neuropsychological measures were used to investigate the impact of APOE-ε status on grey matter structure, white matter integrity, and cognitive functioning.

METHODS

Data were obtained from the Alzheimer's Disease Initiative Phase 3 (ADNI3) database. Baseline MRI, DTI and cognitive composite scores for memory (ADNI-Mem) and executive function (ADNI-EF) were acquired from 116 healthy controls. Participants were grouped according to APOE allele presence (APOE-ε2+ N = 17, APOE-ε3ε3 N = 64, APOE-ε4+ N = 35). Voxel-based morphometry (VBM) and tract-based spatial statistics (TBSS) were used to compare grey matter volume (GMV) and white matter integrity, respectively, between APOE-ε2+ and APOE-ε3ε3 controls, and again between APOE-ε4+ and APOE-ε3ε3 controls. Multivariate analysis of covariance (MANCOVA) was used to examine the effects of APOE polymorphism on memory and EF across all APOE groups with age, sex and education as regressors of no interest. Cognitive scores were correlated (Pearson r) with imaging metrics within groups.

RESULTS

No significant differences were seen across groups, within groups in MRI metrics, or cognitive performance (p>0.05, corrected for multiple comparisons).

CONCLUSIONS

The current study partially replicated and extended previous findings from an earlier multimodal study (Honea 2009). Future studies should clarify APOE mechanisms in healthy ageing by adding other imaging, cognitive, and lifestyle metrics and longitudinal design in larger sample sizes.

[Molecular biology analysis of 2 rare RhD variant individuals with RHD*DEL37]

The Rh blood grouping system is a critical standardized test in transfusion medicine, especially for the cases related to haemolytic transfusion reactions and neonatal haemolytic disease caused by clinical RhD blood group incompatibility. In the present case report, we presented two cases with the uncommon RHD gene variation RHD*DEL37. The blood samples of the two subjects were mistakenly identified as RhD-negative through conventional serological testing. Firstly, both blood samples were tested negative for the RhD antigen using traditional tube test and gel microcolumn methods. The phenotyping of RhCE were identified as ccEe and ccee for each sample, respectively. Secondly, genetic analysis was performed using polymerase chain reaction-sequence specific prime (PCR-SSP) which revealed that neither sample belonging to the several common RHD gene variants which was found in Asia. Moreover, they turned out to be positive for the RHD haplotype, which indicated that exons 1-10 on one of the RHD alleles were entirely absent. In addition, a T>C mutation was observed at bases 1154-31 in intron 8 of the other allele, which was located at the intron 8 breakpoint. This result was obtained after further Sanger sequencing of exons 1-10 of the RHD gene. The mutant allele was designated as RHD*DEL37 by the International Society of Blood Transfusion (ISBT) and was identified as D-elute(Del) by phenotype ana-lysis. Both samples were genotyped as RHD*DEL37 and showed positive results. In summary, the true genotype of the two blood samples, of which the screening results only using serological testing method was negative, were RHD*DEL37 /RHD-(RHD*01N.01). Notably, this kind of genotype was reported for the first time in Chinese population. Moreover, the two individuals did not have ties of consanguinity, indicating that some of the Chinese individuals could be carriers of the genetic mutation. Therefore, it might be necessary to further confirm the frequency of this mutation in the Chinese population and the possibility of homozygosity for this mutation. This report identifies infrequent RHD gene mutation samples by coupling molecular biology and serological methods to prevent misclassification of blood groups. Combining serological and molecular biology test results to determine blood group is critical in protecting patients during clinical transfusion procedures.

KASP: a high-throughput genotyping system and its applications in major crop plants for biotic and abiotic stress tolerance

Advances in plant molecular breeding have resulted in the development of new varieties with superior traits, thus improving the crop germplasm. Breeders can screen a large number of accessions without rigorous and time-consuming phenotyping by marker-assisted selection (MAS). Molecular markers are one of the most imperative tools in plant breeding programmes for MAS to develop new cultivars possessing multiple superior traits. Single nucleotide polymorphisms (SNPs) are ideal for MAS due to their low cost, low genotyping error rates, and reproducibility. Kompetitive Allele Specific PCR (KASP) is a globally recognized technology for SNP genotyping. KASP is an allele-specific oligo extension-based PCR assay that uses fluorescence resonance energy transfer (FRET) to detect genetic variations such as SNPs and insertions/deletions (InDels) at a specific locus. Additionally, KASP allows greater flexibility in assay design, which leads to a higher success rate and the capability to genotype a large population. Its versatility and ease of use make it a valuable tool in various fields, including genetics, agriculture, and medical research. KASP has been extensively used in various plant-breeding applications, such as the identification of germplasm resources, quality control (QC) analysis, allele mining, linkage mapping, quantitative trait locus (QTL) mapping, genetic map construction, trait-specific marker development, and MAS. This review provides an overview of the KASP assay and emphasizes its validation in crop improvement related to various biotic and abiotic stress tolerance and quality traits.

17q21 Variants Disturb Mucosal Host Defense in Childhood Asthma

Rationale: The strongest genetic risk factor for childhood-onset asthma, the 17q21 locus, is associated with increased viral susceptibility and disease-promoting processes.Objectives: To identify biological targets underlying the escalated viral susceptibility associated with the clinical phenotype mediated by the 17q21 locus.Methods: Genome-wide transcriptome analysis of nasal brush samples from 261 children (78 healthy, 79 with wheezing at preschool age, 104 asthmatic) within the ALLIANCE (All-Age-Asthma) cohort, with a median age of 10.0 (range, 1.0-20.0) years, was conducted to explore the impact of their 17q21 genotype (SNP rs72163891). Concurrently, nasal secretions from the same patients and visits were collected, and high-sensitivity mesoscale technology was employed to measure IFN protein levels.Measurements and Main Results: This study revealed that the 17q21 risk allele induces a genotype- and asthma/wheeze phenotype-dependent enhancement of mucosal GSDMB expression as the only relevant 17q21-encoded gene in children with preschool wheeze. Increased GSDMB expression correlated with the activation of a type-1 proinflammatory, cell-lytic immune, and natural killer signature, encompassing key genes linked to an IFN type-2-signature (IFNG, CXCL9, CXCL10, KLRC1, CD8A, GZMA). Conversely, there was a reduction in IFN type 1 and type 3 expression signatures at the mRNA and protein levels.Conclusions: This study demonstrates a novel disease-driving mechanism induced by the 17q21 risk allele. Increased mucosal GSDMB expression is associated with a cell-lytic immune response coupled with compromised airway immunocompetence. These findings suggest that GSDMB-related airway cell death and perturbations in the mucosal IFN signature account for the increased vulnerability of 17q21 risk allele carriers to respiratory viral infections during early life, opening new options for future biological interventions.The All-Age-Asthma (ALLIANCE) cohort is registered at www.clinicaltrials.gov (pediatric arm, NCT02496468).

Genetic mapping reveals new loci and alleles for flowering time and plant height using the double round-robin population of barley

Flowering time and plant height are two critical determinants of yield potential in barley (Hordeum vulgare). Despite their role in plant physiological regulation, a complete overview of the genetic complexity of flowering time and plant height regulation in barley is still lacking. Using a double round-robin population originated from the crossings of 23 diverse parental inbred lines, we aimed to determine the variance components in the regulation of flowering time and plant height in barley as well as to identify new genetic variants by single and multi-population QTL analyses and allele mining. Despite similar genotypic variance, we observed higher environmental variance components for plant height than flowering time. Furthermore, we detected new QTLs for flowering time and plant height. Finally, we identified a new functional allelic variant of the main regulatory gene Ppd-H1. Our results show that the genetic architecture of flowering time and plant height might be more complex than reported earlier and that a number of undetected, small effect, or low-frequency genetic variants underlie the control of these two traits.

Exploring the origin of a unique mutant allele in twin-tail goldfish using CRISPR/Cas9 mutants

Artificial selection has been widely applied to genetically fix rare phenotypic features in ornamental domesticated animals. For many of these animals, the mutated loci and alleles underlying rare phenotypes are known. However, few studies have explored whether these rare genetic mutations might have been fixed due to competition among related mutated alleles or if the fixation occurred due to contingent stochastic events. Here, we performed genetic crossing with twin-tail ornamental goldfish and CRISPR/Cas9-mutated goldfish to investigate why only a single mutated allele-chdS with a E127X stop codon (also called chdAE127X)-gives rise to the twin-tail phenotype in the modern domesticated goldfish population. Two closely related chdS mutants were generated with CRISPR/Cas9 and compared with the E127X allele in F2 and F3 generations. Both of the CRISPR/Cas9-generated alleles were equivalent to the E127X allele in terms of penetrance/expressivity of the twin-tail phenotype and viability of carriers. These findings indicate that multiple truncating mutations could have produced viable twin-tail goldfish. Therefore, the absence of polymorphic alleles for the twin-tail phenotype in modern goldfish likely stems from stochastic elimination or a lack of competing alleles in the common ancestor. Our study is the first experimental comparison of a singular domestication-derived allele with CRISPR/Cas9-generated alleles to understand how genetic fixation of a unique genotype and phenotype may have occurred. Thus, our work may provide a conceptual framework for future investigations of rare evolutionary events in domesticated animals.

Genotype-phenotype associations in CRB1 bi-allelic patients: a novel mutation, a systematic review and meta-analysis

PURPOSE

The goal of the study was to search for novel bi-allelic CRB1 mutations, and then to analyze the CRB1 literature at the genotypic and phenotypic levels.

APPROACH

We screened various variables such as the CRB1 mutation types, domains, exons, and genotypes and their relation with specific ocular phenotypes. An emphasis was given to the bi-allelic missense and nonsense mutations because of their high prevalence compared to other mutation types. Finally, we quantified the effect of various non-modifiable factors over the best-corrected visual acuity oculus uterque (BCVA OU) using multivariate linear regression models and identified genetic interactions.

RESULTS

A novel bi-allelic missense in the exon 9 of CRB1; c.2936G > A; p.(Gly979Asp) was found to be associated with rod-cone dystrophy (RCD). CRB1 mutation type, exons, domains, and genotype distribution varied significantly according to fundus characteristics, such as peripheral pigmentation and condition, optic disc, vessels, macular condition, and pigmentation (P < 0.05). Of the 154 articles retrieved from PubMed, 96 studies with 439 bi-allelic CRB1 patients were included. Missense mutations were significantly associated with an absence of macular pigments, pale optic disc, and periphery pigmentation, resulting in a higher risk of RCD (P < 0.05). In contrast, homozygous nonsense mutations were associated with macular pigments, periphery pigments, and a high risk of LCA (P < 0.05) and increased BCVA OU levels. We found that age, mutation types, and inherited retinal diseases were critical determinants of BCVA OU as they significantly increased it by 33% 26%, and 38%, respectively (P < 0.05). Loss of function alleles additively increased the risk of LCA, with nonsense having a more profound effect than indels. Finally, our analysis showed that p.(Cys948Tyr) and p.(Lys801Ter) and p.(Lys801Ter); p.(Cys896Ter) might interact to modify BCVA OU levels.

CONCLUSION

This meta-analysis updated the literature and identified genotype-phenotype associations in bi-allelic CRB1 patients.

Whole genome sequencing identifies associations for nonsyndromic sagittal craniosynostosis with the intergenic region of BMP2 and noncoding RNA gene LINC01428

Craniosynostosis (CS) is a major birth defect resulting from premature fusion of cranial sutures. Nonsyndromic CS occurs more frequently than syndromic CS, with sagittal nonsyndromic craniosynostosis (sNCS) presenting as the most common CS phenotype. Previous genome-wide association and targeted sequencing analyses of sNCS have identified multiple associated loci, with the strongest association on chromosome 20. Herein, we report the first whole-genome sequencing study of sNCS using 63 proband-parent trios. Sequencing data for these trios were analyzed using the transmission disequilibrium test (TDT) and rare variant TDT (rvTDT) to identify high-risk rare gene variants. Sequencing data were also examined for copy number variants (CNVs) and de novo variants. TDT analysis identified a highly significant locus at 20p12.3, localized to the intergenic region between BMP2 and the noncoding RNA gene LINC01428. Three variants (rs6054763, rs6054764, rs932517) were identified as potential causal variants due to their probability of being transcription factor binding sites, deleterious combined annotation dependent depletion scores, and high minor allele enrichment in probands. Morphometric analysis of cranial vault shape in an unaffected cohort validated the effect of these three single nucleotide variants (SNVs) on dolichocephaly. No genome-wide significant rare variants, de novo loci, or CNVs were identified. Future efforts to identify risk variants for sNCS should include sequencing of larger and more diverse population samples and increased omics analyses, such as RNA-seq and ATAC-seq.

Genetic analysis of 37 cases with primary periodic paralysis in Chinese patients

BACKGROUND

Primary periodic paralysis (PPP) is an inherited disorders of ion channel dysfunction characterized by recurrent episodes of flaccid muscle weakness, which can classified as hypokalemic (HypoPP), normokalemic (NormoPP), or hyperkalemic (HyperPP) according to the potassium level during the paralytic attacks. However, PPP is charactered by remarkable clinical and genetic heterogeneity, and the diagnosis of suspected patients is based on the characteristic clinical presentation then confirmed by genetic testing. At present, there are only limited cohort studies on PPP in the Chinese population.

RESULTS

We included 37 patients with a clinical diagnosis of PPP. Eleven (29.7%) patients were tested using a specific gene panel and 26 (70.3%) by the whole-exome sequencing (WES). Twenty-two cases had a genetic variant identified, representing a diagnostic rate of 59.5% (22/37). All the identified mutations were either in the SCN4A or the CACNA1S gene. The overall detection rate was comparable between the panel (54.5%: 6/11) and WES (61.5%: 16/26). The remaining patients unresolved through panel sequencing were further analyzed by WES, without the detection of any mutation. The novel atypical splicing variant c.2020-5G > A affects the normal splicing of the SCN4A mRNA, which was confirmed by minigene splicing assay. Among 21 patients with HypoPP, 15 patients were classified as HypoPP-2 with SCN4A variants, and 6 HypoPP-1 patients had CACNA1S variants.

CONCLUSIONS

Our results suggest that SCN4A alleles are the main cause in our cohort, with the remainder caused by CACNA1S alleles, which are the predominant cause in Europe and the United States. Additionally, this study identified 3 novel SCN4A and 2 novel CACNA1S variants, broadening the mutation spectrum of genes associated with PPP.

Genomic basis determining root system architecture in maize

KEY MESSAGE

A total of 389 and 344 QTLs were identified by GWAS and QTL mapping explaining accumulatively 32.2-65.0% and 23.7-63.4% of phenotypic variation for 14 shoot-borne root traits using more than 1300 individuals across multiple field trails. Efficient nutrient and water acquisition from soils depends on the root system architecture (RSA). However, the genetic determinants underlying RSA in maize remain largely unexplored. In this study, we conducted a comprehensive genetic analysis for 14 shoot-borne root traits using 513 inbred lines and 800 individuals from four recombinant inbred line (RIL) populations at the mature stage across multiple field trails. Our analysis revealed substantial phenotypic variation for these 14 root traits, with a total of 389 and 344 QTLs identified through genome-wide association analysis (GWAS) and linkage analysis, respectively. These QTLs collectively explained 32.2-65.0% and 23.7-63.4% of the trait variation within each population. Several a priori candidate genes involved in auxin and cytokinin signaling pathways, such as IAA26, ARF2, LBD37 and CKX3, were found to co-localize with these loci. In addition, a total of 69 transcription factors (TFs) from 27 TF families (MYB, NAC, bZIP, bHLH and WRKY) were found for shoot-borne root traits. A total of 19 genes including PIN3, LBD15, IAA32, IAA38 and ARR12 and 19 GWAS signals were overlapped with selective sweeps. Further, significant additive effects were found for root traits, and pyramiding the favorable alleles could enhance maize root development. These findings could contribute to understand the genetic basis of root development and evolution, and provided an important genetic resource for the genetic improvement of root traits in maize.

Replication study identified EFEMP1 association with varicose vein predisposition among Indians

BACKGROUND

Varicose vein is a chronic condition that affects the lower extremities of the human body. Several factors have been implicated in the development of this disease, viz age, gender, weight, height and prolonged standing. Recently, genome-wide studies have identified genetic biomarkers that are associated with varicose veins in different ethnic groups. Such genetic studies are lacking in South Asians specifically in Indians where the prevalence of varicose veins is high, and it is important to replicate these variants in the stated population. The study aimed to replicate the association of genetic variants associated with varicose veins in this target population, which were found to be associated with the other ethnic groups.

METHODOLOGY

The studied cohort is of the Indian population comprising unrelated 104 varicose veins cases and 448 non-varicose vein controls. The samples were genotyped using the Illumina Global Screening Array. Using the genomic data from UK BioBank and 23andMe studied cohorts; eight genetic variants were selected to replicate in our dataset. The allelic association was performed to identify the effective allele and risk was estimated using odds ratio and p-value as level of significance. Multifactor Dimensionality Reduction was used to estimate the cumulative effect of variants in Indians.

RESULT

Variant rs3791679 of EFEMP1 was found to be associated with varicose veins in Indians. After observing the association of the EFEMP1 with varicose veins, we further ensued to identify all genetic variants within EFEMP1 to uncover the additional variants associated with this trait. Interestingly, we identified six new variants of EFEMP1 gene that have shown association. Moreover, the cumulative effect of all associated variations was estimated and the risk was 2.7 times higher in cases than controls whereas independently their effect ranges from 0.37-1.58.

CONCLUSION

This study identifies EFEMP1 as a potential gene related to the risk of varicose veins in Indians. It also highlights that evaluating the maximum number of variants of a gene rather than focusing solely on replicating single variations offers a more comprehensive and nuanced understanding of the genetic factors contributing to a complex trait like varicose veins.

Simultaneous analysis of pMHC binding and reactivity unveils virus-specific CD8 T cell immunity to a concise epitope set

CD8 T cells provide immunity to virus infection through recognition of epitopes presented by peptide major histocompatibility complexes (pMHCs). To establish a concise panel of widely recognized T cell epitopes from common viruses, we combined analysis of TCR down-regulation upon stimulation with epitope-specific enumeration based on barcode-labeled pMHC multimers. We assess CD8 T cell binding and reactivity for 929 previously reported epitopes in the context of 1 of 25 HLA alleles representing 29 viruses. The prevalence and magnitude of CD8 T cell responses were evaluated in 48 donors and reported along with 137 frequently recognized virus epitopes, many of which were underrepresented in the public domain. Eighty-four percent of epitope-specific CD8 T cell populations demonstrated reactivity to peptide stimulation, which was associated with effector and long-term memory phenotypes. Conversely, nonreactive T cell populations were associated primarily with naive phenotypes. Our analysis provides a reference map of epitopes for characterizing CD8 T cell responses toward common human virus infections.

Genetic characterization of a locus responsible for low pungency using EMS-induced mutants in Capsicum annuum L

KEY MESSAGE

The pepper mutants ('221-2-1a' and '1559-1-2h') with very low pungency were genetically characterized. The Pun4 locus, responsible for the reduced pungency of the mutant fruits, was localized to a 208 Mb region on chromosome 6. DEMF06G16460, encoding 3-ketoacyl-CoA synthase, was proposed as a strong candidate gene based on the genetic analyses of bulked segregants, DEG, and expression analyses. Capsaicinoids are unique alkaloids present in pepper (Capsicum spp.), synthesized through the condensation of by-products from the phenylpropanoid and branched-chain fatty acid pathways, and accumulating in the placenta. In this study, we characterized two allelic ethyl methanesulfonate-induced mutant lines with extremely low pungency ('221-2-1a' and '1559-1-2h'). These mutants, derived from the pungent Korean landrace 'Yuwolcho,' exhibited lower capsaicinoid content than Yuwolcho but still contained a small amount of capsaicinoid with functional capsaicinoid biosynthetic genes. Genetic crosses between the mutants and Yuwolcho or pungent lines indicated that a single recessive mutation was responsible for the low-pungency phenotype of mutant 221-2-1a; we named the causal locus Pungency 4 (Pun4). To identify Pun4, we combined genome-wide polymorphism analysis and transcriptome analysis with bulked-segregant analysis. We narrowed down the location of Pun4 to a 208-Mb region on chromosome 6 containing five candidate genes, of which DEMF06G16460, encoding a 3-ketoacyl-CoA synthase associated with branched-chain fatty acid biosynthesis, is the most likely candidate for Pun4. The expression of capsaicinoid biosynthetic genes in placental tissues in Yuwolcho and the mutant was consistent with the branched-chain fatty acid pathway playing a pivotal role in the lower pungency observed in the mutant. We also obtained a list of differentially expressed genes in placental tissues between the mutant and Yuwolcho, from which we selected candidate genes using gene co-expression analysis. In summary, we characterized the capsaicinoid biosynthesis-related locus Pun4 through integrated of genetic, genomic, and transcriptome analyses. These findings will contribute to our understanding of capsaicinoid biosynthesis in pepper.

Protein restriction during pregnancy alters Cdkn1c silencing, dopamine circuitry and offspring behaviour without changing expression of key neuronal marker genes

We tracked the consequences of in utero protein restriction in mice throughout their development and life course using a luciferase-based allelic reporter of imprinted Cdkn1c. Exposure to gestational low-protein diet (LPD) results in the inappropriate expression of paternally inherited Cdkn1c in the brains of embryonic and juvenile mice. These animals were characterised by a developmental delay in motor skills, and by behavioural alterations indicative of reduced anxiety. Exposure to LPD in utero resulted in significantly more tyrosine hydroxylase positive (dopaminergic) neurons in the midbrain of adult offspring as compared to age-matched, control-diet equivalents. Positron emission tomography (PET) imaging revealed an increase in striatal dopamine synthesis capacity in LPD-exposed offspring, where elevated levels of dopamine correlated with an enhanced sensitivity to cocaine. These data highlight a profound sensitivity of the developing epigenome to gestational protein restriction. Our data also suggest that loss of Cdkn1c imprinting and p57KIP2 upregulation alters the cellular composition of the developing midbrain, compromises dopamine circuitry, and thereby provokes behavioural abnormalities in early postnatal life. Molecular analyses showed that despite this phenotype, exposure to LPD solely during pregnancy did not significantly change the expression of key neuronal- or dopamine-associated marker genes in adult offspring.

Biallelic non-productive enhancer-promoter interactions precede imprinted expression of Kcnk9 during mouse neural commitment

It is only partially understood how constitutive allelic methylation at imprinting control regions (ICRs) interacts with other regulation levels to drive timely parental allele-specific expression along large imprinted domains. The Peg13-Kcnk9 domain is an imprinted domain with important brain functions. To gain insights into its regulation during neural commitment, we performed an integrative analysis of its allele-specific epigenetic, transcriptomic, and cis-spatial organization using a mouse stem cell-based corticogenesis model that recapitulates the control of imprinted gene expression during neurodevelopment. We found that, despite an allelic higher-order chromatin structure associated with the paternally CTCF-bound Peg13 ICR, enhancer-Kcnk9 promoter contacts occurred on both alleles, although they were productive only on the maternal allele. This observation challenges the canonical model in which CTCF binding isolates the enhancer and its target gene on either side and suggests a more nuanced role for allelic CTCF binding at some ICRs.

[Genetic analysis of two patients with a rare Ael subtype]

OBJECTIVE

To analyze the genetic sequences of two patients with a rare Ael blood subgroup.

METHODS

Two female patients undergoing treatment respectively for adenomyoma of the uterus and gastritis at the Second Affiliated Hospital, Yuying Children's Hospital of Wenzhou Medical University in June 2019 and September 2020 were selected as the study subjects. Their Ael subtypes were identified with a saline tube agglutination assay and absorption-emission assay. Sequence of the ABO gene Ael subtypes was determined by the Sanger method. The impact of genetic variants on the structural stability of N-acetylgalactosaminyl transferase (GTA) was analyzed with PyMOL software by constructing a structure predicted model.

RESULTS

Both patients were determined as Ael blood subgroup. Sequencing result of patient 1 was ABO*O.01.02/ABO*O.01.02, which has resulted in a p.Thr88Profs*31 amino acid substitution. The sequencing result of patient 2 was ABO*Ael.06/ABO*O.01.02, in which c.425C>T and c.467C>T variants in exon 7 have led to p.Met142Thr and p.Pro156Leu substitutions. Prediction of the protein model speculated that the p.Met142Thr not only can change the binding of GTA protein with water molecules, but also the local hydrogen bond network of GTA, which may lead to decreased enzymatic activity. By contrast, the p.Pro156Leu variant has trivial effect on the structural stability of GTA.

CONCLUSION

The molecular structure of Ael subtypes can be diverse. The genotypes of the two patients have been respectively determined as ABO*O.01.02/ABO*O.01.02 with a G261 deletion and ABO*Ael.06/ABO*O.01.02.

An averaging model for analysis and interpretation of high-order genetic interactions

While combinatorial genetic data collection from biological systems in which quantitative phenotypes are controlled by active and inactive alleles of multiple genes (multi-gene systems) is becoming common, a standard analysis method for such data has not been established. The currently common approaches have three major drawbacks. First, although it is a long tradition in genetics, modeling the effect of an inactive allele (a null mutant allele) contrasted against that of the active allele (the wild-type allele) is not suitable for mechanistic understanding of multi-gene systems. Second, a commonly-used additive model (ANOVA with interaction) mathematically fails in estimation of interactions among more than two genes when the phenotypic response is not linear. Third, interpretation of higher-order interactions defined by an additive model is not intuitive. I derived an averaging model based on algebraic principles to solve all these problems within the framework of a general linear model. In the averaging model: the effect of the active allele is contrasted against the effect of the inactive allele for easier mechanistic interpretations; there is mathematical stability in estimation of higher-order interactions even when the phenotypic response is not linear; and interpretations of higher-order interactions are intuitive and consistent-interactions are defined as the mean effects of the last active genes added to the system. Thus, the key outcomes of this study are development of the averaging model, which is suitable for analysis of multi-gene systems, and a new, intuitive, and mathematically and interpretationally consistent definition of a genetic interaction, which is central to the averaging model.

[Molecular study of an individual with Bel subtype due to a novel c.620T>C variant]

OBJECTIVE

To explore the molecular basis for an individual with Bel subtype of the ABO blood type due to a novel c.620T>C variant gene, and assess its impact on the structure of GTB transferase.

METHODS

An individual who had visited the First Affiliated Hospital of Zhengzhou University on February 11, 2023 was selected as the study subject. ABO phenotyping was initially conducted with serological methods, which was followed by direct sequencing of 7 exons of the ABO gene. Subsequently, single-strand sequencing was carried out by using allele-specific primers, and the variant in the B transferase was homology-modeled using the Modeller software. The impact of the variant on the transferase's spatial structure was analyzed with the PyMOL software.

RESULTS

The serological phenotype of the patient was identified as the Bel subtype. Direct sequencing revealed that she has harbored a novel c.620T>C variant, resulting in a p.Leu207Pro substitution in the polypeptide chain. Combined with single-strand sequencing, her genotype was ultimately determined as ABO*BELnew/ABO*O.01.02. Three-dimensional protein structure modeling showed that, compared with the wild type, the distance of one hydrogen bond between Proline and Glycine at position 272 has increased, along with disappearance of another hydrogen bond.

CONCLUSION

The novel c.620T>C (p.Leu207Pro) variant of B allele may affect the structural stability of the glycosyltransferase. The weakened enzyme activity in turn may lead to reduced B antigen expression, manifesting as the Bel subtype by serological analysis.

BRAF and RET polymorphism association with thyroid cancer risk, a preliminary study from Khyber Pakhtunkhwa population

BACKGROUND

Thyroid cancer, originating in the neck's thyroid gland, encompasses various types. Genetic mutations, particularly in BRAF and RET genes are crucial in its development. This study investigates the association between BRAF (rs113488022) and RET (rs77709286) polymorphisms and thyroid cancer risk in the Khyber Pakhtunkhwa (KP) population.

METHODS

Blood samples from 100 thyroid cancer patients and 100 healthy controls were genotyped using ARMS-PCR followed by gel electrophoresis and statistical analysis.

RESULTS

Analysis revealed a significant association between the minor allele T of BRAF (rs113488022) and thyroid cancer risk (P = 0.0001). Both genotypes of BRAF (rs113488022) showed significant associations with thyroid cancer risk (AT; P = 0.0012 and TT; P = 0.045). Conversely, the minor allele G of RET (rs77709286) exhibited a non-significant association with thyroid cancer risk (P = 0.2614), and neither genotype showed significant associations (CG; P = 0.317, GG; P = 0.651). Demographic and clinical parameters analysis using SPSS showed a non-significant association between BRAF and RET variants and age group (P = 0.878 and P = 0.536), gender (P = 0.587 and P = 0.21), tumor size (P = 0.796 and P = 0.765), or tumor localization (P = 0.689 and P = 0.727).

CONCLUSION

In conclusion, this study emphasizes the significant association between BRAF polymorphism and thyroid cancer risk, while RET polymorphism showed a less pronounced impact. Further validation using larger and specific datasets is essential to establish conclusive results.

Investigation of HLA susceptibility alleles and genotypes with hematological disease among Chinese Han population

Several genes involved in the pathogenesis have been identified, with the human leukocyte antigen (HLA) system playing an essential role. However, the relationship between HLA and a cluster of hematological diseases has received little attention in China. Blood samples (n = 123913) from 43568 patients and 80345 individuals without known pathology were genotyped for HLA class I and II using sequencing-based typing. We discovered that HLA-A*11:01, B*40:01, C*01:02, DQB1*03:01, and DRB1*09:01 were prevalent in China. Furthermore, three high-frequency alleles (DQB1*03:01, DQB1*06:02, and DRB1*15:01) were found to be hazardous in malignant hematologic diseases when compared to controls. In addition, for benign hematologic disorders, 7 high-frequency risk alleles (A*01:01, B*46:01, C*01:02, DQB1*03:03, DQB1*05:02, DRB1*09:01, and DRB1*14:54) and 8 high-frequency susceptible genotypes (A*11:01-A*11:01, B*46:01-B*58:01, B*46:01-B*46:01, C*01:02-C*03:04, DQB1*03:01-DQB1*05:02, DQB1*03:03-DQB1*06:01, DRB1*09:01-DRB1*15:01, and DRB1*14:54-DRB1*15:01) were observed. To summarize, our findings indicate the association between HLA alleles/genotypes and a variety of hematological disorders, which is critical for disease surveillance.

Plasmodium falciparum genetic diversity and multiplicity of infection based on msp-1, msp-2, glurp and microsatellite genetic markers in sub-Saharan Africa: a systematic review and meta-analysis

BACKGROUND

In sub-Saharan Africa (SSA), Plasmodium falciparum causes most of the malaria cases. Despite its crucial roles in disease severity and drug resistance, comprehensive data on Plasmodium falciparum genetic diversity and multiplicity of infection (MOI) are sparse in SSA. This study summarizes available information on genetic diversity and MOI, focusing on key markers (msp-1, msp-2, glurp, and microsatellites). The systematic review aimed to evaluate their influence on malaria transmission dynamics and offer insights for enhancing malaria control measures in SSA.

METHODS

The review was conducted following the Preferred Reporting Items for Systematic Review and Meta-Analysis (PRISMA) guidelines. Two reviewers conducted article screening, assessed the risk of bias (RoB), and performed data abstraction. Meta-analysis was performed using the random-effects model in STATA version 17.

RESULTS

The review included 52 articles: 39 cross-sectional studies and 13 Randomized Controlled Trial (RCT)/cohort studies, involving 11,640 genotyped parasite isolates from 23 SSA countries. The overall pooled mean expected heterozygosity was 0.65 (95% CI: 0.51-0.78). Regionally, values varied: East (0.58), Central (0.84), Southern (0.74), and West Africa (0.69). Overall pooled allele frequencies of msp-1 alleles K1, MAD20, and RO33 were 61%, 44%, and 40%, respectively, while msp-2 I/C 3D7 and FC27 alleles were 61% and 55%. Central Africa reported higher frequencies (K1: 74%, MAD20: 51%, RO33: 48%) than East Africa (K1: 46%, MAD20: 42%, RO33: 31%). For msp-2, East Africa had 60% and 55% for I/C 3D7 and FC27 alleles, while West Africa had 62% and 50%, respectively. The pooled allele frequency for glurp was 66%. The overall pooled mean MOI was 2.09 (95% CI: 1.88-2.30), with regional variations: East (2.05), Central (2.37), Southern (2.16), and West Africa (1.96). The overall prevalence of polyclonal Plasmodium falciparum infections was 63% (95% CI: 56-70), with regional prevalences as follows: East (62%), West (61%), Central (65%), and South Africa (71%).

CONCLUSION

The study shows substantial regional variation in Plasmodium falciparum parasite genetic diversity and MOI in SSA. These findings suggest a need for malaria control strategies and surveillance efforts considering regional-specific factors underlying Plasmodium falciparum infection.