Robust analysis of allele-specific copy number alterations from scRNA-seq data with XClone.

Somatic copy number alterations (CNAs) are major mutations that contribute to the development and progression of various cancers. Despite a few computational methods proposed to detect CNAs from single-cell transcriptomic data, the technical sparsity of such data makes it challenging to identify allele-specific CNAs, particularly in complex clonal structures. In this study, we present a statistical method, XClone, that strengthens the signals of read depth and allelic imbalance by effective smoothing on cell neighborhood and gene coordinate graphs to detect haplotype-aware CNAs from scRNA-seq data. By applying XClone to multiple datasets with challenging compositions, we demonstrated its ability to robustly detect different types of allele-specific CNAs and potentially indicate whole genome duplication, therefore enabling the discovery of corresponding subclones and the dissection of their phenotypic impacts.

The first high-altitude autotetraploid haplotype-resolved genome assembled (Rhododendron nivale subsp. boreale) provides new insights into mountaintop adaptation.

BACKGROUND: Rhododendron nivale subsp. boreale Philipson et M. N. Philipson is an alpine woody species with ornamental qualities that serve as the predominant species in mountainous scrub habitats found at an altitude of ∼4,200 m. As a high-altitude woody polyploid, this species may serve as a model to understand how plants adapt to alpine environments. Despite its ecological significance, the lack of genomic resources has hindered a comprehensive understanding of its evolutionary and adaptive characteristics in high-altitude mountainous environments. FINDINGS: We sequenced and assembled the genome of R. nivale subsp. boreale, an assembly of the first subgenus Rhododendron and the first high-altitude woody flowering tetraploid, contributing an important genomic resource for alpine woody flora. The assembly included 52 pseudochromosomes (scaffold N50 = 42.93 Mb; BUSCO = 98.8%; QV = 45.51; S-AQI = 98.69), which belonged to 4 haplotypes, harboring 127,810 predicted protein-coding genes. Conjoint k-mer analysis, collinearity assessment, and phylogenetic investigation corroborated autotetraploid identity. Comparative genomic analysis revealed that R. nivale subsp. boreale originated as a neopolyploid of R. nivale and underwent 2 rounds of ancient polyploidy events. Transcriptional expression analysis showed that differences in expression between alleles were common and randomly distributed in the genome. We identified extended gene families and signatures of positive selection that are involved not only in adaptation to the mountaintop ecosystem (response to stress and developmental regulation) but also in autotetraploid reproduction (meiotic stabilization). Additionally, the expression levels of the (group VII ethylene response factor transcription factors) ERF VIIs were significantly higher than the mean global gene expression. We suspect that these changes have enabled the success of this species at high altitudes. CONCLUSIONS: We assembled the first high-altitude autopolyploid genome and achieved chromosome-level assembly within the subgenus Rhododendron. In addition, a high-altitude adaptation strategy of R. nivale subsp. boreale was reasonably speculated. This study provides valuable data for the exploration of alpine mountaintop adaptations and the correlation between extreme environments and species polyploidization.

Allopolyploidization from two dioecious ancestors leads to recurrent evolution of sex chromosomes.

Polyploidization presents an unusual challenge for species with sex chromosomes, as it can lead to complex combinations of sex chromosomes that disrupt reproductive development. This is particularly true for allopolyploidization between species with different sex chromosome systems. Here, we assemble haplotype-resolved chromosome-level genomes of a female allotetraploid weeping willow (Salix babylonica) and a male diploid S. dunnii. We show that weeping willow arose from crosses between a female ancestor from the Salix-clade, which has XY sex chromosomes on chromosome 7, and a male ancestor from the Vetrix-clade, which has ancestral XY sex chromosomes on chromosome 15. We find that weeping willow has one pair of sex chromosomes, ZW on chromosome 15, that derived from the ancestral XY sex chromosomes in the male ancestor of the Vetrix-clade. Moreover, the ancestral 7X chromosomes from the female ancestor of the Salix-clade have reverted to autosomal inheritance. Duplicated intact ARR17-like genes on the four homologous chromosomes 19 likely have contributed to the maintenance of dioecy during polyploidization and sex chromosome turnover. Taken together, our results suggest the rapid evolution and reversion of sex chromosomes following allopolyploidization in weeping willow.

Genome-Wide Association and RNA-Seq Analyses Reveal a Potential Candidate Gene Related to Oil Content in Soybean Seeds.

Soybean is a crucial crop globally, serving as a significant source of unsaturated fatty acids and protein in the human diet. However, further enhancements are required for the related genes that regulate soybean oil synthesis. In this study, 155 soybean germplasms were cultivated under three different environmental conditions, followed by phenotypic identification and genome-wide association analysis using simplified sequencing data. Genome-wide association analysis was performed using SLAF-seq data. A total of 36 QTLs were significantly associated with oil content (-log10(p) > 3). Out of the 36 QTLs associated with oil content, 27 exhibited genetic overlap with previously reported QTLs related to oil traits. Further transcriptome sequencing was performed on extreme high-low oil soybean varieties. Combined with transcriptome expression data, 22 candidate genes were identified (|log2FC| ≥ 3). Further haplotype analysis of the potential candidate genes showed that three potential candidate genes had excellent haplotypes, including Glyma.03G186200, Glyma.09G099500, and Glyma.18G248900. The identified loci harboring beneficial alleles and candidate genes likely contribute significantly to the molecular network's underlying marker-assisted selection (MAS) and oil content.

VolcanoSV enables accurate and robust structural variant calling in diploid genomes from single-molecule long read sequencing.

Structural variants (SVs) significantly contribute to human genome diversity and play a crucial role in precision medicine. Although advancements in single-molecule long-read sequencing offer a groundbreaking resource for SV detection, identifying SV breakpoints and sequences accurately and robustly remains challenging. We introduce VolcanoSV, an innovative hybrid SV detection pipeline that utilizes both a reference genome and local de novo assembly to generate a phased diploid assembly. VolcanoSV uses phased SNPs and unique k-mer similarity analysis, enabling precise haplotype-resolved SV discovery. VolcanoSV is adept at constructing comprehensive genetic maps encompassing SNPs, small indels, and all types of SVs, making it well-suited for human genomics studies. Our extensive experiments demonstrate that VolcanoSV surpasses state-of-the-art assembly-based tools in the detection of insertion and deletion SVs, exhibiting superior recall, precision, F1 scores, and genotype accuracy across a diverse range of datasets, including low-coverage (10x) datasets. VolcanoSV outperforms assembly-based tools in the identification of complex SVs, including translocations, duplications, and inversions, in both simulated and real cancer data. Moreover, VolcanoSV is robust to various evaluation parameters and accurately identifies breakpoints and SV sequences.

HINT1 Gene Polymorphisms, Smoking Behaviour, and Personality Traits: A Haplotype Case-Control Study.

The factors influencing the development and maintenance of nicotine dependence are numerous and complex. Recent studies indicate that smokers exhibit distinct genetic predispositions to nicotine dependence. We aimed to analyse (1) the association between rs2551038 and cigarette smoking, (2) the association of between the rs3864236-rs2526303-rs2551038 haplotype and cigarette smoking, and (3) the personality traits measured by the NEO Five-Factor Inventory in cigarette users and never-smokers. No significant differences were present in the frequency of rs2551038 genotypes and alleles in the studied cigarette users compared to the control group. Cigarette users, compared to the control group, had higher scores on the NEO-FFI Extraversion scale (p = 0.0011), and lower scores were obtained by the cigarette users for the NEO-FFI Openness (p = 0.0060), Agreeability (p ≤ 0.000), and Conscientiousness (p ≤ 0.000) scales. There was a significant positive Pearson's linear correlation between the age and the Fagestrom test (r = 0.346; p < 0.0001) and the NEO-FFI Openness scale (r = 0.180; p < 0.0001) in the group of cigarette users. We observed significant linkage disequilibrium between rs2526303 and rs3864236 (D' = 0.3581; p < 2.2204 × 10-16) and between rs2526303 and rs2551038 (D' = 0.9993; p < 2.2204 × 10-16) in the tested sample. The sex-stratified haplotype analysis revealed that in the group of male never-smokers, the GTC haplotype was significantly more frequent than in the group of cigarette users (38% vs. 22%; p = 0.0039). The presented study reveals significant differences in personality trait scores between cases and controls. Moreover, the sex-stratified analysis showed significant differences in haplotype distribution. These results underscore the interplay between genetic predisposition, sex, and personality in nicotine-using individuals.

Association between CYP1A2 gene variants -163 C/A (rs762551) and -3860 G/A (rs2069514) and bladder cancer susceptibility.

BACKGROUND: Bladder cancer (BLCA) poses a significant global health challenge due to its high incidence, poor prognosis, and limited treatment options. AIMS AND OBJECTIVES: This study aims to investigate the association between two specific polymorphisms, CYP1A2-163 C/A and CYP1A2-3860G/A, within the Cytochrome P450 1A2 (CYP1A2) gene and susceptibility to BLCA. METHODS: The study employed a case-control design, genotyping 340 individuals using Polymerase Chain Reaction-High-Resolution Melting Curve (PCR-HRM). Various genetic models were applied to evaluate allele and genotype frequencies. Genetic linkage analysis was facilitated using R packages. RESULTS: The study reveals a significant association with the - 163 C/A allele, particularly in the additive model. Odds ratio (OR) analysis links CYP1A2-163 C/A (rs762551) and CYP1A2-3860G/A(rs2069514) polymorphisms to BLCA susceptibility. The rs762551 C/A genotype is prevalent in 55% of BLCA cases and exhibits an OR of 2.21. The A/A genotype has an OR of 1.54. Regarding CYP1A2-3860G/A, the G/A genotype has an OR of 1.54, and the A/A genotype has an OR of 2.08. Haplotype analysis shows a predominant C-C haplotype at 38.2%, followed by a C-A haplotype at 54.7%, and a less frequent A-A haplotype at 7.1%. This study underscores associations between CYP1A2 gene variants, particularly rs762551 (CYP1A2-163 C/A), and an increased susceptibility to BLCA. Haplotype analysis of 340 individuals reveals a predominant C-C haplotype at 38.2%, followed by a C-A haplotype at 54.7%, and a less frequent A-A haplotype at 7.1%. CONCLUSION: In conclusion, the - 163 C/A allele, C/A genotype of rs762551, and G/A genotype of rs2069514 emerge as potential genetic markers associated with elevated BLCA risk.

Concurrent of compound heterozygous variant of a novel in-frame deletion and the common hypomorphic haplotype in TBX6 and inherited 17q12 microdeletion in a fetus.

BACKGROUND: TBX6, a member of the T-box gene family, encodes the transcription factor box 6 that is critical for somite segmentation in vertebrates. It is known that the compound heterozygosity of disruptive variants in trans with a common hypomorphic risk haplotype (T-C-A) in the TBX6 gene contribute to 10% of congenital scoliosis (CS) cases. The deletion of chromosome 17q12 is a rare cytogenetic abnormality, which often leads to renal cysts and diabetes mellitus. However, the affected individuals often exhibit clinical heterogeneity and incomplete penetrance. METHODS: We here present a Chinese fetus who was shown to have CS by ultrasound examination at 17 weeks of gestation. Trio whole-exome sequencing (WES) was performed to investigate the underlying genetic defects of the fetus. In vitro functional experiments, including western-blotting and luciferase transactivation assay, were performed to determine the pathogenicity of the novel variant of TBX6. RESULTS: WES revealed the fetus harbored a compound heterozygous variant of c.338_340del (p.Ile113del) and the common hypomorphic risk haplotype of the TBX6 gene. In vitro functional study showed the p.Ile113del variant had no impact on TBX6 expression, but almost led to complete loss of its transcriptional activity. In addition, we identified a 1.85 Mb deletion on 17q12 region in the fetus and the mother. Though there is currently no clinical phenotype associated with this copy number variation in the fetus, it can explain multiple renal cysts in the pregnant woman. CONCLUSIONS: This study is the first to report a Chinese fetus with a single amino acid deletion variant and a T-C-A haplotype of TBX6. The clinical heterogeneity of 17q12 microdeletion poses significant challenges for prenatal genetic counseling. Our results once again suggest the complexity of prenatal genetic diagnosis.

Detecting haplotype-specific transcript variation in long reads with FLAIR2.

BACKGROUND: RNA-seq has brought forth significant discoveries regarding aberrations in RNA processing, implicating these RNA variants in a variety of diseases. Aberrant splicing and single nucleotide variants (SNVs) in RNA have been demonstrated to alter transcript stability, localization, and function. In particular, the upregulation of ADAR, an enzyme that mediates adenosine-to-inosine editing, has been previously linked to an increase in the invasiveness of lung adenocarcinoma cells and associated with splicing regulation. Despite the functional importance of studying splicing and SNVs, the use of short-read RNA-seq has limited the community's ability to interrogate both forms of RNA variation simultaneously. RESULTS: We employ long-read sequencing technology to obtain full-length transcript sequences, elucidating cis-effects of variants on splicing changes at a single molecule level. We develop a computational workflow that augments FLAIR, a tool that calls isoform models expressed in long-read data, to integrate RNA variant calls with the associated isoforms that bear them. We generate nanopore data with high sequence accuracy from H1975 lung adenocarcinoma cells with and without knockdown of ADAR. We apply our workflow to identify key inosine isoform associations to help clarify the prominence of ADAR in tumorigenesis. CONCLUSIONS: Ultimately, we find that a long-read approach provides valuable insight toward characterizing the relationship between RNA variants and splicing patterns.

Candidate Genes and Favorable Haplotypes Associated with Iron Toxicity Tolerance in Rice.

Iron (Fe) toxicity is a major issue adversely affecting rice production worldwide. Unfortunately, the physiological and genetic mechanisms underlying Fe toxicity tolerance in rice remain relatively unknown. In this study, we conducted a genome-wide association study using a diverse panel consisting of 551 rice accessions to identify genetic mechanisms and candidate genes associated with Fe toxicity tolerance. Of the 29 quantitative trait loci (QTL) for Fe toxicity tolerance detected on chromosomes 1, 2, 5, and 12, five (qSH_Fe5, qSFW_Fe2.3, qRRL5.1, qRSFW1.1, and qRSFW12) were selected to identify candidate genes according to haplotype and bioinformatics analyses. The following five genes were revealed as promising candidates: LOC_Os05g40160, LOC_Os05g40180, LOC_Os12g36890, LOC_Os12g36900, and LOC_Os12g36940. The physiological characteristics of rice accessions with contrasting Fe toxicity tolerance reflected the importance of reactive oxygen species-scavenging antioxidant enzymes and Fe homeostasis for mitigating the negative effects of Fe toxicity on rice. Our findings have clarified the genetic and physiological mechanisms underlying Fe toxicity tolerance in rice. Furthermore, we identified valuable genetic resources for future functional analyses and the development of Fe toxicity-tolerant rice varieties via marker-assisted selection.

Association between OX40L polymorphism and type 2 diabetes mellitus in Iranians.

INTRODUCTION: Diabetes mellitus (DM) is one of the leading causes of morbidity and mortality worldwide. It is a multifactorial disease that genetic and environmental factors contribute to its development. The aim of the study was to investigate the association of OX40L promoter gene polymorphisms with type 2 diabetes mellitus (T2DM) in Iranians. MATERIALS AND METHODS: Three hundred and sixty-eight subjects including 184 healthy subjects and 184 T2DM patients were enrolled in our study. Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) was applied to detect genotype and allele frequencies of rs3850641, rs1234313 and rs10912580. In addition, SNPStats web tool was applied to estimate haplotype frequency and linkage disequilibrium (LD). RESULTS: The distribution of tested polymorphisms was statistically different between the T2DM patients and healthy subjects (P < 0.01). rs1234313 AG (OR = 0.375, 95% CI = 0.193-0.727, P = 0.004) and rs10912580 AG (OR = 0.351, 95% CI = 0.162-0.758, P = 0.008) genotypes were associated with the decreased risk of T2DM in Iranians. Moreover, our prediction revealed that AAG (OR = 0.46, 95% CI= (0.28-0.76), P = 0.0028) and GAG (OR = 0.24, 95% CI= (0.13-0.45), P < 0.0001) haplotypes were related to the reduced risk of the disease. However, the tested polymorphisms had no effect on biochemical parameters and body mass index (BMI) in the patient group (P > 0.05). CONCLUSION: Our findings revealed that OX40L promoter gene polymorphisms are associated with T2DM. Moreover, genotype and allelic variations were related to the decreased risk of T2DM in Iranians. Further studies are recommended to show whether these polymorphic variations could affect OX40/OX40L interaction or OX40L phenotype.

Two decades of molecular surveillance in Senegal reveal rapid changes in known drug resistance mutations over time.

BACKGROUND: Drug resistance in Plasmodium falciparum is a major threat to malaria control efforts. Pathogen genomic surveillance could be invaluable for monitoring current and emerging parasite drug resistance. METHODS: Data from two decades (2000-2020) of continuous molecular surveillance of P. falciparum parasites from Senegal were retrospectively examined to assess historical changes in malaria drug resistance mutations. Several known drug resistance markers and their surrounding haplotypes were profiled using a combination of single nucleotide polymorphism (SNP) molecular surveillance and whole genome sequence based population genomics. RESULTS: This dataset was used to track temporal changes in drug resistance markers whose timing correspond to historically significant events such as the withdrawal of chloroquine (CQ) and the introduction of sulfadoxine-pyrimethamine (SP) in 2003. Changes in the mutation frequency at Pfcrt K76T and Pfdhps A437G coinciding with the 2014 introduction of seasonal malaria chemoprevention (SMC) in Senegal were observed. In 2014, the frequency of Pfcrt K76T increased while the frequency of Pfdhps A437G declined. Haplotype-based analyses of Pfcrt K76T showed that this rapid increase was due to a recent selective sweep that started after 2014. DISCUSSION (CONCLUSION): The rapid increase in Pfcrt K76T is troubling and could be a sign of emerging amodiaquine (AQ) resistance in Senegal. Emerging AQ resistance may threaten the future clinical efficacy of artesunate-amodiaquine (ASAQ) and AQ-dependent SMC chemoprevention. These results highlight the potential of molecular surveillance for detecting rapid changes in parasite populations and stress the need to monitor the effectiveness of AQ as a partner drug for artemisinin-based combination therapy (ACT) and for chemoprevention.

Optimisation of global stem cell donor recruitment based on analysis of unsuccessful donor searches.

Despite over 41 million registered potential volunteer stem cell donors worldwide, many patients in need of a transplant do not find an HLA-matched unrelated donor or cord blood units, with the respective odds differing significantly between various populations. In this study, we analysed data of 2205 unsuccessful real-life donor searches sent to the DKMS Registry to identify populations in which further donor recruitment would be associated with particularly large patient benefits. For that purpose, we estimated haplotype frequencies of 67 donor populations at various sample sizes and entered them into two different mathematical models. These models assessed patient benefits from population-specific donor recruitment, operationalised by the number of originally unsuccessful searches that may become successful due to new donors. Consistently, across the different mathematical models and sample sizes, we obtained several countries from East and Southeast Asia (Thailand, Vietnam, China, and the Philippines) and the population of Asians in the USA as countries/populations where donor recruitment activities would be particularly beneficial for patients. We also identified various countries in Southeast and Central Europe as possible target regions for donor recruitment with above-average patient benefits. The results presented are registry-specific in the sense that they were obtained by optimising unsuccessful searches that had been sent to the DKMS Registry. Therefore, it would be desirable to apply the presented methods to a global data set that includes all unsuccessful stem cell donor searches worldwide and uses population-specific haplotype frequencies based on all donors available in the WMDA Search & Match Service.

Structural and genetic diversity in the secreted mucins MUC5AC and MUC5B.

The secreted mucins MUC5AC and MUC5B are large glycoproteins that play critical defensive roles in pathogen entrapment and mucociliary clearance. Their respective genes contain polymorphic and degenerate protein-coding variable number tandem repeats (VNTRs) that make the loci difficult to investigate with short reads. We characterize the structural diversity of MUC5AC and MUC5B by long-read sequencing and assembly of 206 human and 20 nonhuman primate (NHP) haplotypes. We find that human MUC5B is largely invariant (5,761-5,762 amino acids [aa]); however, seven haplotypes have expanded VNTRs (6,291-7,019 aa). In contrast, 30 allelic variants of MUC5AC encode 16 distinct proteins (5,249-6,325 aa) with cysteine-rich domain and VNTR copy-number variation. We group MUC5AC alleles into three phylogenetic clades: H1 (46%, ∼5,654 aa), H2 (33%, ∼5,742 aa), and H3 (7%, ∼6,325 aa). The two most common human MUC5AC variants are smaller than NHP gene models, suggesting a reduction in protein length during recent human evolution. Linkage disequilibrium and Tajima's D analyses reveal that East Asians carry exceptionally large blocks with an excess of rare variation (p < 0.05) at MUC5AC. To validate this result, we use Locityper for genotyping MUC5AC haplogroups in 2,600 unrelated samples from the 1000 Genomes Project. We observe a signature of positive selection in H1 among East Asians and a depletion of the likely ancestral haplogroup (H3). In Europeans, H3 alleles show an excess of common variation and deviate from Hardy-Weinberg equilibrium (p < 0.05), consistent with heterozygote advantage and balancing selection. This study provides a generalizable strategy to characterize complex protein-coding VNTRs for improved disease associations.

The novel CFTR haplotype E583G/F508del in CFTR-related disorder.

BACKGROUND: CFTR-related disorder (CFTR-RD) is a clinical entity associated to complex diagnostic paths and newly upgraded standard of care. In CFTR-RD, CFTR genotyping represents a diagnostic surrogate marker. In case of novel haplotype, the diagnosis could represents an area of concern. We described the molecular evaluation of the rare CFTR variant E583G identified in trans with the F508del in a novel haplotype. METHODS AND RESULTS: An adult woman was referred to our pulmonary unit for persistent respiratory symptoms. CFTR Next Generation Sequencing was performed to evaluate full-gene mutational status. The variant identified was evaluated for its pathogenicity integrating clinical evidences with dedicated bioinformatics analyses. Clinical evaluation of patient matched with a mono-organ CFTR-RD diagnosis. Genotyping revealed the novel CFTR haplotype F508del/E583G. Multiple evidences of a deleterious effect of the CFTR E583G rare variant emerged from the bioinformatics analyses performed. CONCLUSIONS: Guidelines for CFTR-RD are available with the purpose of harmonizing clinical and molecular investigations. In such context, the identification of novel CFTR haplotype need to a deeper evaluation with a combination of skills. The novel E583G variant could be considered of clinical interest and overall a CFTR-RD Variants of Varying Clinical Consequences.

Orthanq: transparent and uncertainty-aware haplotype quantification with application in HLA-typing.

BACKGROUND: Identification of human leukocyte antigen (HLA) types from DNA-sequenced human samples is important in organ transplantation and cancer immunotherapy and remains a challenging task considering sequence homology and extreme polymorphism of HLA genes. RESULTS: We present Orthanq, a novel statistical model and corresponding application for transparent and uncertainty-aware quantification of haplotypes. We utilize our approach to perform HLA typing while, for the first time, reporting uncertainty of predictions and transparently observing mutations beyond reported HLA types. Using 99 gold standard samples from 1000 Genomes, Illumina Platinum Genomes and Genome In a Bottle projects, we show that Orthanq can provide overall superior accuracy and shorter runtimes than state-of-the-art HLA typers. CONCLUSIONS: Orthanq is the first approach that allows to directly utilize existing pangenome alignments and type all HLA loci. Moreover, it can be generalized for usages beyond HLA typing, e.g. for virus lineage quantification. Orthanq is available under https://orthanq.github.io .

Mutation accumulation in H. sapiens F508del CFTR countermands dN/dS type genomic analysis.

Understanding the mechanisms that underlie de novo mutations (DNMs) can be essential for interpreting human evolution, including aspects such as rapidly diverging genes, conservation of non-coding regulatory elements, and somatic DNA adaptation, among others. DNM accumulation in Homo sapiens is often limited to evaluation of human trios or quads across a single generation. Moreover, human SNPs in exons, pseudogenes, or other non-coding elements can be ancient and difficult to date, including polymorphisms attributable to founder effects and identity by descent. In this report, we describe multigenerational evolution of a human coding locus devoid of natural selection, and delineate patterns and principles by which DNMs have accumulated over the past few thousand years. We apply a data set comprising cystic fibrosis transmembrane conductance regulator (CFTR) alleles from 2,393 individuals homozygous for the F508del defect. Additional polymorphism on the F508del background diversified subsequent to a single mutational event during recent human history. Because F508del CFTR is without function, SNPs observed on this haplotype are effectively attributable to factors that govern accumulating de novo mutations. We show profound enhancement of transition, synonymous, and positionally repetitive polymorphisms, indicating appearance of DNMs in a manner evolutionarily designed to protect protein coding DNA against mutational attrition while promoting diversity.

Common and novel haplotype structures between different types of cancer.

BACKGROUND: Background: Genome-wide association studies (GWAS) have identified hundreds of genetic variants associated with cancer risk. GWAS data are important for cancer prevention and understanding the underlying mechanisms of cancer. AIMS: This study aimed to investigate the genetic association between different types of cancer using GWAS data and a bioinformatics approach. METHODS AND RESULTS: The significant GWAS variants associated with more than one cancer type were identified. Common linkage disequilibrium (LD) variants between different types of cancer were identified by 1000 genomes phase 3 LD data. Haplotype blocks were identified by analyzing 1000 Genomes phase 3 genotyping data in the GWAS populations. Subsequent analyses included functional SNP analyses and TCGA gene expression. The results associated with significant GWAS variants (P<5E-8) showed the following haplotype associations in European population: GT rs4808075-rs8170 haplotype on BABAM1 with breast and ovarian cancers, GC rs16857609-rs11693806 haplotype on DIRC3 with breast and thyroid cancers, GCG rs380286-rs401681-rs31487 haplotype on CLPTM1L with skin and lung cancers, GGG rs4430796-rs11651052-rs11263763 haplotype on HNF1B with prostate and endometrial cancers, and GT rs10505477-rs6983267 haplotype on CASC8 associated with colorectal and prostate cancers. All these genes had significantly different expressions in tumor tissues (P<1E-3). In addition, the rs11693806 variant is located in the hsa-miR-873-5p binding site and has an enhancing effect on the hsa-miR-873-5p:DIRC3 interaction. CONCLUSION: These novel haplotype structures and miRNA:lncRNA interactions are important for understanding the common genetic link between cancers. These results can potentially be used in genetic panels.

[Genetic study of a rare Chinese pedigree with a recombination occurring between the HLA-A/C loci in both parents].

OBJECTIVE: To analyze a Chinese pedigree with a recombination occurring between the HLA-A/C loci in both parents. METHODS: A patient who was planning to undergo hematopoietic stem cell transplantation due to "aplastic anemia" in February 2022 was selected as the study subject. Peripheral blood samples were collected from the patient, his parents and brother. HLA-A/C/B/DRB1/DQB1 high-resolution typing was carried out by using sequence-based typing and sequence-specific oligonucleotides. The recombination was identified by pedigree analysis. The HLA haplotype of each individual was identified by genealogical analysis. The parentage possibility was determined by short tandem repeat analysis. HLA-A/C/B/DRB1/DRB345/DQA1/DQB1/DPA1/DPB1 were determined with next-generation high-throughput sequence-based typing. The recombination sites were analyzed by family study. RESULTS: The high parentage possibilities of the family was confirmed by short tandem repeat analysis. Recombination was found between the HLA-A*24:02 A*33:03/C*14:03 in the paternally transmitted haplotype, whilst HLA-A*01:01 A*03:01/C*08:02 was found in the maternally transmitted haplotype, which had resulted in two novel HLA haplotypes in the proband. CONCLUSION: A rare case with simultaneous recombination of the paternal and maternal HLA-A/C loci has been discovered, which may facilitate further study of the mechanisms of the HLA recombination.

The haplotype-resolved genome of diploid Chrysanthemum indicum unveils new acacetin synthases genes and their evolutionary history.

Acacetin, a flavonoid compound, possesses a wide range of pharmacological effects, including antimicrobial, immune regulation, and anticancer effects. Some key steps in its biosynthetic pathway were largely unknown in flowering plants. Here, we present the first haplotype-resolved genome of Chrysanthemum indicum, whose dried flowers contain abundant flavonoids and have been utilized as traditional Chinese medicine. Various phylogenetic analyses revealed almost equal proportion of three tree topologies among three Chrysanthemum species (C. indicum, C. nankingense, and C. lavandulifolium), indicating that frequent gene flow among Chrysanthemum species or incomplete lineage sorting due to rapid speciation might contribute to conflict topologies. The expanded gene families in C. indicum were associated with oxidative functions. Through comprehensive candidate gene screening, we identified five flavonoid O-methyltransferase (FOMT) candidates, which were highly expressed in flowers and whose expressional levels were significantly correlated with the content of acacetin. Further experiments validated two FOMTs (CI02A009970 and CI03A006662) were capable of catalyzing the conversion of apigenin into acacetin, and these two genes are possibly responsible acacetin accumulation in disc florets and young leaves, respectively. Furthermore, combined analyses of ancestral chromosome reconstruction and phylogenetic trees revealed the distinct evolutionary fates of the two validated FOMT genes. Our study provides new insights into the biosynthetic pathway of flavonoid compounds in the Asteraceae family and offers a model for tracing the origin and evolutionary routes of single genes. These findings will facilitate in vitro biosynthetic production of flavonoid compounds through cellular and metabolic engineering and expedite molecular breeding of C. indicum cultivars.