Genomic Regions Research Articles

Characterizing Genetic Susceptibility to Colorectal Cancer in Taiwan Through Genome-Wide Association Study.

We conducted the first genome-wide association study (GWAS) of colorectal cancer (CRC) in Taiwan with 5342 cases and 61,015 controls. Ninety-two SNPs in three genomic regions reached genome-wide significance (p < 5 × 10-8). The lead SNPs in these three regions were: rs12778523 (OR = 1.18, 95% CI, 1.15-1.23, p = 4.51 × 10-13), an intergenic SNP between RNA5SP299 and LINC02676 at chromosome 10p14; rs647161 (OR = 1.14, 95% CI, 1.09-1.19, p = 2.21 × 10-9), an intronic SNP in PITX1 at 5q31.1, and rs10427139 (OR = 1.20, 95% CI, 1.14-1.28, p = 3.62 × 10-9), an intronic SNP in GPATCH1 at 19q13.1. We further validated CRC susceptibility SNPs previously identified through GWAS in other populations. A total of 61 CRC susceptibility SNPs were confirmed in Taiwanese. The top validated putative CRC susceptibility genes included: POU2AF2, HAO1, LAMC1, EIF3H, BMP2, ZMIZ1, BMP4, POLD3, CDKN1A, PREX1, CDKN2B, CDH1, and LRIG1. The top enriched pathways included TGF-β signaling, BMP signaling, extracellular matrix organization, DNA repair, and cell cycle control. We could not validate SNPs in HLA-G at 6p22.1 and in NOTCH4 at 6p21.32. We generated a weighted genetic risk score (GRS) using the 61 SNPs and constructed receiver operating characteristic (ROC) curves using the GRS to predict CRC. The area under the ROC curve (AUC) was 0.589 for GRS alone and 0.645 for GRS, sex, and age. These susceptibility SNPs and genes provide important insights into the molecular mechanisms of CRC development and help identify high-risk individuals for CRC in Taiwan.

Regulatory genome annotation of 33 insect species

Annotation of newly sequenced genomes frequently includes genes, but rarely covers important non-coding genomic features such as the cis-regulatory modules—e.g., enhancers and silencers—that regulate gene expression. Here, we begin to remedy this situation by developing a workflow for rapid initial annotation of insect regulatory sequences, and provide a searchable database resource with enhancer predictions for 33 genomes. Using our previously developed SCRMshaw computational enhancer prediction method, we predict over 2.8 million regulatory sequences along with the tissues where they are expected to be active, in a set of insect species ranging over 360 million years of evolution. Extensive analysis and validation of the data provides several lines of evidence suggesting that we achieve a high true-positive rate for enhancer prediction. One, we show that our predictions target specific loci, rather than random genomic locations. Two, we predict enhancers in orthologous loci across a diverged set of species to a significantly higher degree than random expectation would allow. Three, we demonstrate that our predictions are highly enriched for regions of accessible chromatin. Four, we achieve a validation rate in excess of 70% using in vivo reporter gene assays. As we continue to annotate both new tissues and new species, our regulatory annotation resource will provide a rich source of data for the research community and will have utility for both small-scale (single gene, single species) and large-scale (many genes, many species) studies of gene regulation. In particular, the ability to search for functionally related regulatory elements in orthologous loci should greatly facilitate studies of enhancer evolution even among distantly related species.

Regulatory genome annotation of 33 insect species.

Annotation of newly sequenced genomes frequently includes genes, but rarely covers important non-coding genomic features such as the cis-regulatory modules-e.g., enhancers and silencers-that regulate gene expression. Here, we begin to remedy this situation by developing a workflow for rapid initial annotation of insect regulatory sequences, and provide a searchable database resource with enhancer predictions for 33 genomes. Using our previously developed SCRMshaw computational enhancer prediction method, we predict over 2.8 million regulatory sequences along with the tissues where they are expected to be active, in a set of insect species ranging over 360 million years of evolution. Extensive analysis and validation of the data provides several lines of evidence suggesting that we achieve a high true-positive rate for enhancer prediction. One, we show that our predictions target specific loci, rather than random genomic locations. Two, we predict enhancers in orthologous loci across a diverged set of species to a significantly higher degree than random expectation would allow. Three, we demonstrate that our predictions are highly enriched for regions of accessible chromatin. Four, we achieve a validation rate in excess of 70% using in vivo reporter gene assays. As we continue to annotate both new tissues and new species, our regulatory annotation resource will provide a rich source of data for the research community and will have utility for both small-scale (single gene, single species) and large-scale (many genes, many species) studies of gene regulation. In particular, the ability to search for functionally related regulatory elements in orthologous loci should greatly facilitate studies of enhancer evolution even among distantly related species.

Region-based analyses of existing genome-wide association studies identifies novel potential genetic susceptibility regions for glioma.

Glioma is a rare and debilitating brain cancer with one of the lowest cancer survival rates. Genome-wide association studies have identified 34 genetic susceptibility regions. We sought to discover novel susceptibility regions using approaches which test groups of contiguous genetic markers simultaneously. We analyzed data from three independent glioma studies of European ancestry, GliomaScan (1,316 cases/1,293 controls), AGOG (560 cases/2,237 controls), and GICC (4,000 cases/2,411 controls), using the machine-learning algorithm DEPTH and a region-based regression method based on the generalized Berk-Jones (GBJ) statistic, to assess the association of glioma with genomic regions by glioma type and sex. Summary statistics from the UCSF/Mayo Clinic study were used for independent validation. We conducted a meta-analysis using GliomaScan, AGOG, GICC and UCSF/Mayo. We identified 11 novel candidate genomic regions for glioma risk common to multiple studies. Two of the 11 regions, 16p13.3 containing RBFOX1 and 1p36.21 containing PRDM2, were significantly associated with female and male glioma risk respectively, based on results of the meta-analysis. Both regions have been previously linked to glioma tumor progression. Three of the 11 regions contain neurotransmitter receptor genes (7q31.33 GRM8, 5q35.2 DRD1, 15q13.3 CHRNA7). Our region-based approach identified 11 genomic regions that suggest association with glioma risk of which two regions, 16p13.3 and 1p36.21, warrant further investigation as genetic susceptibility regions for female and male risk respectively. Our analyses suggest that genetic susceptibility to glioma may differ by sex and highlights the possibility that synapse-related genes play a role in glioma susceptibility.

PAM-relaxed and temperature-tolerant CRISPR-Mb3Cas12a single transcript unit systems for efficient singular and multiplexed genome editing in rice, maize, and tomato.

Class 2 Type V-A CRISPR-Cas (Cas12a) nucleases are powerful genome editing tools, particularly effective in A/T-rich genomic regions, complementing the widely used CRISPR-Cas9 in plants. To enhance the utility of Cas12a, we investigate three Cas12a orthologs-Mb3Cas12a, PrCas12a, and HkCas12a-in plants. Protospacer adjacent motif (PAM) requirements, editing efficiencies, and editing profiles are compared in rice. Among these orthologs, Mb3Cas12a exhibits high editing efficiency at target sites with a simpler, relaxed TTV PAM which is less restrictive than the canonical TTTV PAM of LbCas12a and AsCas12a. To optimize Mb3Cas12a, we develop an efficient single transcription unit (STU) system by refining the linker between Mb3Cas12a and CRISPR RNA (crRNA), nuclear localization signal (NLS), and direct repeat (DR). This optimized system enables precise genome editing in rice, particularly for fine-tuning target gene expression by editing promoter regions. Further, we introduced Arginine (R) substitutions at Aspartic acid (D) 172, Asparagine (N) 573, and Lysine (K) 579 of Mb3Cas12a, creating two temperature-tolerant variants: Mb3Cas12a-R (D172R) and Mb3Cas12a-RRR (D172R/N573R/K579R). These variants demonstrate significantly improved editing efficiency at lower temperatures (22 °C and 28 °C) in rice cells, with Mb3Cas12a-RRR showing the best performance. We extend this approach by developing efficient Mb3Cas12a-RRR STU systems in maize and tomato, achieving biallelic mutants targeting single or multiple genes in T0 lines cultivated at 28 °C and 25 °C, respectively. This study significantly expands Cas12a's targeting capabilities in plant genome editing, providing valuable tools for future research and practical applications.

Purifying Selection Influences the Comparison of Heterozygosities between Populations.

Heterozygosity is a fundamental measure routinely used to compare between populations to infer the level of genetic variation and their relative effective population sizes. However, such comparison is highly influenced by the magnitude of selection pressure on the genomic regions used. Using over 2 million Single Nucleotide Variants (SNVs) from chimpanzee and mouse populations, this study shows that the heterozygosities estimated using neutrally evolving sites of large populations were two times higher than those of small populations. However, this difference was only ~1.6 times for the heterozygosities estimated using nonsynonymous sites. This suggests an excess in the nonsynonymous heterozygosities due to the segregation of deleterious variants in small populations. This excess in the nonsynonymous heterozygosities of the small populations was estimated to be 23-31%. Further analysis revealed that the magnitude of the excess is modulated by effective population size (Ne) and selection intensity (s). Using chimpanzee populations, this investigation found that the excess in nonsynonymous diversity in the small population was little (6%) when the difference between the Ne values of large and small populations was small (2.4 times). Conversely, this was high (23%) when the difference in Ne was large (5.9 times). Analysis using mouse populations showed that the excess in the nonsynonymous diversity of highly constrained genes of the small population was much higher (38%) than that observed for the genes under relaxed selective constraints (21%). Similar results were observed when the expression levels of genes were used as a proxy for selection intensity. These results emphasize the use of neutral regions, less constrained genes, or lowly expressed genes when comparing the heterozygosities between populations.

Genetic Dissection of Genomic Regions Associated with Resistance to Urdbean Leaf Crinkle Disease and Yield Related Traits in Black Gram (Vigna mungo (L.) Hepper)

Genetic Dissection of Genomic Regions Associated with Resistance to Urdbean Leaf Crinkle Disease and Yield Related Traits in Black Gram (Vigna mungo (L.) Hepper)

Benchmarking reveals superiority of deep learning variant callers on bacterial nanopore sequence data.

Variant calling is fundamental in bacterial genomics, underpinning the identification of disease transmission clusters, the construction of phylogenetic trees, and antimicrobial resistance detection. This study presents a comprehensive benchmarking of variant calling accuracy in bacterial genomes using Oxford Nanopore Technologies (ONT) sequencing data. We evaluated three ONT basecalling models and both simplex (single-strand) and duplex (dual-strand) read types across 14 diverse bacterial species. Our findings reveal that deep learning-based variant callers, particularly Clair3 and DeepVariant, significantly outperform traditional methods and even exceed the accuracy of Illumina sequencing, especially when applied to ONT's super-high accuracy model. ONT's superior performance is attributed to its ability to overcome Illumina's errors, which often arise from difficulties in aligning reads in repetitive and variant-dense genomic regions. Moreover, the use of high-performing variant callers with ONT's super-high accuracy data mitigates ONT's traditional errors in homopolymers. We also investigated the impact of read depth on variant calling, demonstrating that 10× depth of ONT super-accuracy data can achieve precision and recall comparable to, or better than, full-depth Illumina sequencing. These results underscore the potential of ONT sequencing, combined with advanced variant calling algorithms, to replace traditional short-read sequencing methods in bacterial genomics, particularly in resource-limited settings.

The pathogenic germline ETV4 P433L mutation identified in multiple primary lung cancer affect tumor stem-like property by Wnt/β-catenin pathway

The occurrence of multiple primary lung cancer (MPLC) has witnessed a significant surge in recent years within the Chinese population. MPLC is distinguished by its potential genetic susceptibility and notable genetic heterogeneity. Investigating the etiology of MPLC holds substantial clinical importance.The whole genome sequencing (WGS) and genome-wide linkage analysis were performed in a family affected by a dominant form of lung abnormalities. Specifically, five family members were diagnosed with MPLC, while nine members had pulmonary nodules and one normal member. To confirm the potential pathogenic germline mutations sites, Sanger sequencing was performed in an additional 162 MPLC family patients. Furthermore, molecular biology experiments were conducted to investigate the function and the mechanism of the identified pathogenic mutation site in lung cancer A549 and H322, both in vitro and in vivo. Linkage analysis revealed the presence of shared genomic regions among affected family members. Subsequent exome sequencing identified a deleterious variant within these linkage intervals, specifically a heterozygous mutation in ETS-oncogene transcription factors 4 (ETV4). This particular variant was found in affected family members at a rate of 13 out of 15 individuals. Furthermore, ETV4 P433L mutation could be detected in an additional MPLC family patients and mutation frequency was 3.7% (6 out of 162). The ETV4 P433L mutations site was introduced into lung cancer cell lines, resulting in altered migration and stem-like properties of the cancer cells. Further investigation revealed that the activation of the Wnt/β-catenin signaling pathway, which is associated with stemness, could be attributed to the presence of the ETV4 P433L mutation, suggesting its involvement in tumor promotion. A novel pathogenic germline mutation, ETV4 P433L, was identified in a dominant MPLC family, with a mutation rate of 3.7% among MPLC family patients. The ETV4 P433L mutation was found to impact the stem-like properties and migration of tumors through Wnt/β-catenin signaling pathway.

Drought-responsive dynamics of H3K9ac-marked 3D chromatin interactions are integrated by OsbZIP23-associated super-enhancer-like promoter regions in rice

BackgroundIn response to drought stress (DS), plants undergo complex processes that entail significant transcriptome reprogramming. However, the intricate relationship between the dynamic alterations in the three-dimensional (3D) genome and the modulation of gene co-expression in drought responses remains a relatively unexplored area.ResultsIn this study, we reconstruct high-resolution 3D genome maps based on genomic regions marked by H3K9ac, an active histone modification that dynamically responds to soil water variations in rice. We discover a genome-wide disconnection of 3D genome contact upon DS with over 10,000 chromatin loops lost, which are partially recovered in the subsequent re-watering. Loops integrating promoter–promoter interactions (PPI) contribute to gene expression in addition to basal H3K9ac modifications. Moreover, H3K9ac-marked promoter regions with high affinities in mediating PPIs, termed as super-promoter regions (SPRs), integrate spatially clustered PPIs in a super-enhancer-like manner. Interestingly, the knockout mutation of OsbZIP23, a well-defined DS-responsive transcription factor, leads to the disassociation of over 80% DS-specific PPIs and decreased expression of the corresponding genes under DS. As a case study, we show how OsbZIP23 integrates the PPI cluster formation and the co-expression of four dehydrin genes, RAB16A–D, through targeting the RAB16C SPR in a stress signaling-dependent manner.ConclusionsOur high-resolution 3D genome maps unveil the principles and details of dynamic genome folding in response to water supply variations and illustrate OsbZIP23 as an indispensable integrator of the yet unique 3D genome organization that is essential for gene co-expression under DS in rice.

Genome-Wide Association Studies (GWAS) and Transcriptome Analysis Reveal Male Heterogametic Sex-Determining Regions and Candidate Genes in Northern Snakeheads (Channa argus).

The Northern snakehead (Channa argus) is a significant economic aquaculture species in China. Exhibiting sexual dimorphism in the growth rate between females and males, mono-sex breeding holds substantial value for aquaculture. This study employed GWAS and transcriptome analysis were applied to identify sex determination genomic regions and develop sex-specific markers. A total of 270 single-nucleotide polymorphisms (SNPs) and 31 insertion-deletions (InDels) were identified as being sexually dimorphic through GWAS and fixation index (Fst) scanning. Based on GWAS results, two sex-specific InDel markers were developed, effectively distinguishing genetic sex for XX females, XY males, and YY super-males via (polymerase chain reaction) PCR amplification. A major genomic segment of approximately 115 kb on chromosome 3 (Chr 03) was identified as the sex-determination region. A comparative transcriptome analysis of gonads for three sexes identified 158 overlapping differentially expressed genes (DEGs). Additionally, three sex-related candidate genes were identified near the sex determination region, including id2, sox11, and rnf144a. Further studies are required to elucidate the functions of these genes. Overall, two sex-specific InDel markers support a male heterogametic XX/XY sex-determination system in Northern snakeheads and three candidate genes offer new insights into sex determination and the evolution of sex chromosomes in teleost fish.

MethylSeqLogo: DNA methylation smart sequence logos

BackgroundSome transcription factors, MYC for example, bind sites of potentially methylated DNA. This may increase binding specificity as such sites are (1) highly under-represented in the genome, and (2) offer additional, tissue specific information in the form of hypo- or hyper-methylation. Fortunately, bisulfite sequencing data can be used to investigate this phenomenon.MethodWe developed MethylSeqLogo, an extension of sequence logos which includes new elements to indicate DNA methylation and under-represented dimers in each position of a set binding sites. Our method displays information from both DNA strands, and takes into account the sequence context (CpG or other) and genome region (promoter versus whole genome) appropriate to properly assess the expected background dimer frequency and level of methylation. MethylSeqLogo preserves sequence logo semantics—the relative height of nucleotides within a column represents their proportion in the binding sites, while the absolute height of each column represents information (relative entropy) and the height of all columns added together represents total informationResultsWe present figures illustrating the utility of using MethylSeqLogo to summarize data from several CpG binding transcription factors. The logos show that unmethylated CpG binding sites are a feature of transcription factors such as MYC and ZBTB33, while some other CpG binding transcription factors, such as CEBPB, appear methylation neutral.ConclusionsOur software enables users to explore bisulfite and ChIP sequencing data sets—and in the process obtain publication quality figures.

Histone modifications of circulating nucleosomes are associated with changes in cell-free DNA fragmentation patterns

The analysis of tissues of origin of cell-free DNA (cfDNA) is of research and diagnostic interest. Many studies focused on bisulfite treatment or immunoprecipitation protocols to assess the tissues of origin of cfDNA. DNA loss often occurs during such processes. Fragmentomics of cfDNA molecules has uncovered a wealth of information related to tissues of origin of cfDNA. There is still much room for the development of tools for assessing contributions from various tissues into plasma using fragmentomic features. Hence, we developed an approach to analyze the relative contributions of DNA from different tissues into plasma, by identifying characteristic fragmentation patterns associated with selected histone modifications. We named this technique as FRAGmentomics-based Histone modification Analysis (FRAGHA). Deduced placenta-specific histone H3 lysine 27 acetylation (H3K27ac)-associated signal correlated well with the fetal DNA fraction in maternal plasma (Pearson's r = 0.96). The deduced liver-specific H3K27ac-associated signal correlated with the donor-derived DNA fraction in liver transplantation recipients (Pearson's r = 0.92) and was significantly increased in patients with hepatocellular carcinoma (HCC) (P < 0.01, Wilcoxon rank-sum test). Significant elevations of erythroblasts-specific and colon-specific H3K27ac-associated signals were observed in patients with β-thalassemia major and colorectal cancer, respectively. Furthermore, using the fragmentation patterns from tissue-specific H3K27ac regions, a machine learning algorithm was developed to enhance HCC detection, with an area under the curve (AUC) of up to 0.97. Finally, genomic regions with H3K27ac or histone H3 lysine 4 trimethylation (H3K4me3) were found to exhibit different fragmentomic patterns of cfDNA. This study has shed light on the relationship between cfDNA fragmentomics and histone modifications, thus expanding the armamentarium of liquid biopsy.

Genetic Heterogeneity Across Dimensions of Alcohol Use Behaviors.

Increasingly large samples in genome-wide association studies (GWASs) for alcohol use behaviors (AUBs) have led to an influx of implicated genes, yet the clinical and functional understanding of these associations remains low, in part because most GWASs do not account for the complex and varied manifestations of AUBs. This study applied a multidimensional framework to investigate the latent genetic structure underlying heterogeneous dimensions of AUBs. Multimodal assessments (self-report, interview, electronic health records) were obtained from approximately 400,000 UK Biobank participants. GWAS was conducted for 18 distinct AUBs, including consumption, drinking patterns, alcohol problems, and clinical sequelae. Latent genetic factors were identified and carried forward to GWAS using genomic structural equation modeling, followed by functional annotation, genetic correlation, and enrichment analyses to interpret the genetic associations. Four latent factors were identified: Problems, Consumption, BeerPref (declining alcohol consumption with a preference for drinking beer), and AtypicalPref (drinking fortified wine and spirits). The latent factors were moderately correlated (rg values, 0.12-0.57) and had distinct patterns of associations, with BeerPref in particular implicating many novel genomic regions. Patterns of regional and cell type-specific gene expression in the brain also differed between the latent factors. Deep phenotyping is an important next step to improve understanding of the genetic etiology of AUBs, in addition to increasing sample size. Further effort is required to uncover the genetic heterogeneity underlying AUBs using methods that account for their complex, multidimensional nature.

Changes in the genome of tick-borne encephalitis virus during cultivation

The tick-borne encephalitis virus (TBEV) of the Siberian genotype was previously isolated from the brain of a deceased person. TBEV variants obtained at passages 3 and 8 on SPEV cells were inoculated into the brains of white mice for subsequent passages. Full-genome sequences of all virus variants were analyzed by high-throughput sequencing. The analysis showed the presence of 41 point nucleotide substitutions, which were mainly localized in the genes of non-structural proteins NS3 and NS5 of TBEV. In the deduced virus protein sequences, 12 amino acid substitutions were identified. After three passages through mouse brains, reverse nucleotide and amino acid substitutions were detected. Most of them were mapped in the NS5 protein gene, where 5 new nucleotide substitutions also appeared. At the same time, there was an increase in the length of the 3′-untranslated region (3′-UTR) of the viral genome by 306 nucleotides. The Y3 and Y2 3’-UTR elements were found to contain imperfect L and R repeats, which probably associated with inhibition of the activity of cellular XRN1 RNase and thus involved in the formation of sfRNA1 and sfRNA2. For all TBEV variants, high levels of infectious virus were detected both in cell culture and in the brain of white mice. The revealed changes in the genome that occur during successive passages of TBEV are most likely due to the significant genetic variability of the virus, which ensures its efficient reproduction in different hosts and active circulation in nature.

Hidden evolutionary constraints dictate the retention of coronavirus accessory genes.

Coronaviruses exhibit many mechanisms of genetic innovation, including the acquisition of accessory genes that originate by capture of cellular genes or through duplication of existing viral genes. Accessory genes influence viral host range and cellular tropism, but little is known about how selection acts on these variable regions of virus genomes. We used experimental evolution of mouse hepatitis virus (MHV) encoding a cellular AKAP7 phosphodiesterase and an inactive native phosphodiesterase, NS2 to model the evolutionary fate of accessory genes. After courses of serial infection, the gene encoding inactive NS2, ORF2, unexpectedly remained intact, suggesting it is under cryptic constraint uncoupled from the function of NS2. In contrast, AKAP7 was retained under strong selection but rapidly lost under relaxed selection. Experimental evolution also led to altered viral replication in a cell type-specific manner and changed the relative proportions of subgenomic viral RNA in plaque-purified viral isolates, revealing additional mechanisms of adaptation. Guided by the retention of ORF2 and similar patterns in related betacoronaviruses, we analyzed ORF8 of SARS-CoV-2, which arose via gene duplication and contains premature stop codons in several globally successful lineages. As with MHV ORF2, the coding-defective SARS-CoV-2 ORF8 gene remains largely intact, mirroring patterns observed during MHV experimental evolution, challenging assumptions on the dynamics of gene loss in virus genomes and extending these findings to viruses currently adapting to humans.

The multifaceted roles of phosphoethanolamine-modified lipopolysaccharides: from stress response and virulence to cationic antimicrobial resistance.

SUMMARYLipopolysaccharides (LPS) are an integral part of the outer membrane of Gram-negative bacteria and play essential structural and functional roles in maintaining membrane integrity as well as in stress response and virulence. LPS comprises a membrane-anchored lipid A group, a sugar-based core region, and an O-antigen formed by repeating oligosaccharide units. 3-Deoxy-D-manno-octulosonic acid-lipid A (Kdo2-lipid A) is the minimum LPS component required for bacterial survival. While LPS modifications are not essential, they play multifaceted roles in stress response and host-pathogen interactions. Gram-negative bacteria encode several distinct LPS-modifying phosphoethanolamine transferases (PET) that add phosphoethanolamine (pEtN) to lipid A or the core region of LPS. The pet genes differ in their genomic locations, regulation mechanisms, and modification targets of the encoded enzyme, consistent with their various roles in different growth niches and under varied stress conditions. The discovery of mobile colistin resistance genes, which represent lipid A-modifying pet genes that are encoded on mobile elements and associated with resistance to the last-resort antibiotic colistin, has led to substantial interest in PETs and pEtN-modified LPS over the last decade. Here, we will review the current knowledge of the functional diversity of pEtN-based LPS modifications, including possible roles in niche-specific fitness advantages and resistance to host-produced antimicrobial peptides, and discuss how the genetic and structural diversities of PETs may impact their function. An improved understanding of the PET group will further enhance our comprehension of the stress response and virulence of Gram-negative bacteria and help contextualize host-pathogen interactions.

Genome-wide association studies of COVID-19 vaccine seroconversion and breakthrough outcomes in UK Biobank

Understanding the genetic basis of COVID-19 vaccine seroconversion is crucial to study the role of genetics on vaccine effectiveness. In our study, we used UK Biobank data to find the genetic determinants of COVID-19 vaccine-induced seropositivity and breakthrough infections. We conducted four genome-wide association studies among vaccinated participants for COVID-19 vaccine seroconversion and breakthrough susceptibility and severity. Our findings confirmed a link between the HLA region and seroconversion after the first and second doses. Additionally, we identified 10 genomic regions associated with breakthrough infection (SLC6A20, ST6GAL1, MUC16, FUT6, MXI1, MUC4, HMGN2P18-KRTCAP2, NFKBIZ and APOC1), and one with breakthrough severity (APOE). No significant evidence of genetic colocalisation was found between those traits. Our study highlights the roles of individual genetic make-up in the varied antibody responses to COVID-19 vaccines and provides insights into the potential mechanisms behind breakthrough infections occurred even after the vaccination.

Phenotypic, genetic, variation, and molecular function of CaMYB113 in pepper (Capsicum annuum L.).

Pepper (Capsicum annuum L.) is widely consumed vegetables worldwide, and F1 hybrids are highly sought after in the pepper seed industry. However, studies on gene mutations affecting the color of cotyledon are rare, and the same is true for peppers. In this study, a segregating population was developed by crossing the pepper accession 21C1344 with purple cotyledon and accession 21C912 with green cotyledon. Initially, a target genomic region was identified by screening polymorphic SSR markers distributed across 12 chromosomes. Subsequently, polymorphic markers were developed based on resequencing data from the two parental lines, and genetic linkage analysis was performed. This approach ultimately identified Capana10g001433 (CaMYB113) as the candidate gene responsible for the purple cotyledons. The gene mutation type in 21C912 represents a new mutation type distinct from the reported missense mutation types, and this mutation affects the biosynthesis of anthocyanins. Virus-induced gene silencing (VIGS) of CaMYB113 substantially decreased anthocyanin accumulation in the cotyledons. Subsequent overexpression of CaMYB113 resulted in purple callus and leaves of pepper, and changed the expression levels of downstream genes involved in anthocyanin synthesis. Yeast one-hybrid and dual-luciferase transient expression assays demonstrated the binding of CaMYB113 to anthocyanin biosynthesis-related genes, thereby regulating anthocyanin accumulation in pepper cotyledons.

A T cell receptor specific for an HLA-A*03:01-restricted epitope in the endogenous retrovirus ERV-K-Env exhibits limited recognition of its cognate epitope

Transposable elements (TEs) are often expressed at higher levels in tumor cells than normal cells, implicating these genomic regions as an untapped pool of tumor-associated antigens. In ovarian cancer (OC), protein from the TE ERV-K is frequently expressed by tumor cells. Here we determined whether the targeting of previously identified epitope in the envelope gene (env) of ERV-K resulted in target antigen specificity against cancer cells. We found that transducing healthy donor T cells with an ERV-K-Env-specific T cell receptor construct resulted in antigen specificity only when co-cultured with HLA-A*03:01 B lymphoblastoid cells. Furthermore, in vitro priming of several healthy donors with this epitope of ERV-K-Env did not result in target antigen specificity. These data suggest that the T cell receptor is a poor candidate for targeting this specific ERV-K-Env epitope and has limited potential as a T cell therapy for OC.