Genomic Regions Research Articles

Integrated meta-QTL and in silico transcriptome assessment pinpoint major genomic regions responsible for spike length in wheat (Triticum aestivum L.).

Spike length (SL) is one of the major contributors to wheat yield. Uncovering major genetic regions affecting SL is an integral part of elucidating the genetic basis of wheat yield traits and goes further pivotal for marker-assisted selection breeding. A genome-wide meta-quantitative trait locus (MQTL) analysis of wheat SL resulted in the refinement of 48 MQTLs using 227 initial QTLs retrieved from previous studies published over the past decades. The average confidence interval (CI) of these MQTLs amounted to a 5.16-fold reduction compared to the mean CI of the initial QTLs. As many as 2240 putative candidate genes (CGs) were identified from the MQTL intervals using transcriptomics data in silico of wheat, of which 58 CGs were identified based on wheat-rice homology analysis. For the key CG affecting SL, a functional kompetitive allele-specific PCR (KASP) marker, TaPP2C-3B-KASP, was developed to distinguish TaPP2C-3B-Hap I and TaPP2C-3B-Hap II based on the single nucleotide polymorphism at the 272 bp (A/G). The frequency of the elite allelic variation TaPP2C-3B-Hap II with high SL remained relatively stable at about 49.62% from the 1960s to 1990s. Integration of MQTL analysis and in silico transcriptome data led to a significant increase in the reliability of CGs for the genetic regulation of wheat SL, and the haplotype analysis for key CGs TaPP2C-3B of SL provided insights into the biological function of the TaPP2C-3B gene.

Quantitative trait loci concentrate in specific regions of the Mexican cavefish genome and reveal key candidate genes for cave-associated evolution.

A major goal of modern biology is connecting phenotype with its underlying genetic basis. The Mexican cavefish (Astyanax mexicanus), a characin fish species comprised of a surface ecotype and a cave-derived ecotype, is well suited as a model to study the genetic mechanisms underlying adaptation to extreme environments. Here we map 206 previously published quantitative trait loci (QTL) for cave-derived traits in A. mexicanus to the newest version of the surface fish genome assembly, AstMex3. These analyses revealed that QTL cluster in the genome more than expected by chance, and this clustering is not explained by the distribution of genes in the genome. To investigate whether certain characteristics of the genome facilitate phenotypic evolution, we tested whether genomic characteristics associated with increased opportunities for mutation, such as highly mutagenic CpG sites, are reliable predictors of the sites of trait evolution but did not find any significant trends. Finally, we combined the QTL map with previously collected expression and selection data to identify 36 candidate genes that may underlie the repeated evolution of cave phenotypes, including rgrb, which is predicted to be involved in phototransduction. We found this gene has disrupted exons in all non-hybrid cave populations but intact reading frames in surface fish. Overall, our results suggest specific regions of the genome may play significant roles in driving adaptation to the cave environment in Astyanax mexicanus and demonstrate how this compiled dataset can facilitate our understanding of the genetic basis of repeated evolution in the Mexican cavefish.

Genome-wide association study for polledness, horn shape, and wool traits in Original Valachian sheep

Abstract. The Original Valachian sheep is an endangered Slovak national breed that is well adapted to high-altitude pastures. The sheep can be horned with various shapes and can have multi-coloured or completely white or black wool. Breeders are interested in learning about the genetic basis of these traits. We conducted a genome-wide association study based on the genomic information of 96 sheep genotyped by the GeneSeek GGP Ovine 50K SNP (single-nucleotide polymorphism) chip and on the following traits: polledness (presence or absence of horns), horn shape, and wool colour (completely white and completely black). The univariate linear mixed model was used to discover genetic variants significantly associated with tested traits. The Bonferroni correction and the false-discovery rate were used as significance thresholds. The RXFP2 gene (chromosome 10, 29.5 Mb) was identified as a strong candidate for polledness. In addition, when compared to animals with sideways-turned horns vs. polled, the region around the ADAMTS3 gene (chromosome 6, 88.47 Mb) was significant. A total of nine significant genomic regions were found when comparing the sideways-turned spiral horns with the backwards-curled horns, the two most frequent horn types in Original Valachian sheep. The RXFP2 may also contribute to the genetic control of horn shape. Genes identified in other regions were involved to osteogenic differentiation and osteoblast proliferation (PCP4, chromosome 1, 260.7 Mb), bone mineral density and mineral content (NKX1-2, chromosome 22, 43.75 Mb). The significant genetic variants close to the region of MC1R (chromosome 14, at 14.2 Mb) were associated with the wool colour of sheep that were fully white or fully black animals. The results of this study will contribute to a better understanding of the phenotypic variability of the Original Valachian sheep, especially regarding traits that are very important for breeders of this endangered breed.

Genetic diversity of United States Rambouillet, Katahdin and Dorper sheep

BackgroundManaging genetic diversity is critically important for maintaining species fitness. Excessive homozygosity caused by the loss of genetic diversity can have detrimental effects on the reproduction and production performance of a breed. Analysis of genetic diversity can facilitate the identification of signatures of selection which may contribute to the specific characteristics regarding the health, production and physical appearance of a breed or population. In this study, breeds with well-characterized traits such as fine wool production (Rambouillet, N = 745), parasite resistance (Katahdin, N = 581) and environmental hardiness (Dorper, N = 265) were evaluated for inbreeding, effective population size (Ne), runs of homozygosity (ROH) and Wright’s fixation index (FST) outlier approach to identify differential signatures of selection at 36,113 autosomal single nucleotide polymorphisms (SNPs).ResultsKatahdin sheep had the largest current Ne at the most recent generation estimated with both the GONe and NeEstimator software. The most highly conserved ROH Island was identified in Rambouillet with a signature of selection on chromosome 6 containing 202 SNPs called in an ROH in 50 to 94% of the individuals. This region contained the DCAF16, LCORL and NCAPG genes that have been previously reported to be under selection and have biological roles related to milk production and growth traits. The outlier regions identified through the FST comparisons of Katahdin with Rambouillet and Dorper contained genes with known roles in milk production and mastitis resistance or susceptibility, and the FST comparisons of Rambouillet with Katahdin and Dorper identified genes related to wool growth, suggesting these traits have been under natural or artificial selection pressure in these populations. Genes involved in the cytokine-cytokine receptor interaction pathways were identified in all FST breed comparisons, which indicates the presence of allelic diversity between these breeds in genomic regions controlling cytokine signaling mechanisms.ConclusionsIn this paper, we describe signatures of selection within diverse and economically important U.S. sheep breeds. The genes contained within these signatures are proposed for further study to understand their relevance to biological traits and improve understanding of breed diversity.

Exploring the roles of RNAs in chromatin architecture using deep learning

Recent studies have highlighted the impact of both transcription and transcripts on 3D genome organization, particularly its dynamics. Here, we propose a deep learning framework, called AkitaR, that leverages both genome sequences and genome-wide RNA-DNA interactions to investigate the roles of chromatin-associated RNAs (caRNAs) on genome folding in HFFc6 cells. In order to disentangle the cis- and trans-regulatory roles of caRNAs, we have compared models with nascent transcripts, trans-located caRNAs, open chromatin data, or DNA sequence alone. Both nascent transcripts and trans-located caRNAs improve the models’ predictions, especially at cell-type-specific genomic regions. Analyses of feature importance scores reveal the contribution of caRNAs at TAD boundaries, chromatin loops and nuclear sub-structures such as nuclear speckles and nucleoli to the models’ predictions. Furthermore, we identify non-coding RNAs (ncRNAs) known to regulate chromatin structures, such as MALAT1 and NEAT1, as well as several new RNAs, RNY5, RPPH1, POLG-DT and THBS1-IT1, that might modulate chromatin architecture through trans-interactions in HFFc6. Our modeling also suggests that transcripts from Alus and other repetitive elements may facilitate chromatin interactions through trans R-loop formation. Our findings provide insights and generate testable hypotheses about the roles of caRNAs in shaping chromatin organization.

Cell-type-specific effects of autism-associated 15q duplication syndrome in the human brain

Recurrent copy-number variation represents one of the most well-established genetic drivers in neurodevelopmental disorders, including autism spectrum disorder. Duplication of 15q11-q13 (dup15q) is a well-described neurodevelopmental syndrome that increases the risk of autism more than 40-fold. However, the effects of this duplication on gene expression and chromatin accessibility in specific cell types in the human brain remain unknown. To identify the cell-type-specific transcriptional and epigenetic effects of dup15q in the human frontal cortex, we conducted single-nucleus RNA sequencing and multi-omic sequencing on dup15q-affected individuals (n= 6) as well as individuals with non-dup15q autism (n= 7) and neurotypical control individuals (n= 7). Cell-type-specific differential expression analysis identified significantly regulated genes, critical biological pathways, and differentially accessible genomic regions. Although there was overall increased gene expression across the duplicated genomic region, cellular identity represented an important factor mediating gene-expression changes. As compared to other cell types, neuronal subtypes showed greater upregulation of gene expression across a critical region within the duplication. Genes that fell within the duplicated region and had high baseline expression in control individuals showed only modest changes in dup15q, regardless of cell type. Of note, dup15q and autism had largely distinct signatures of chromatin accessibility but shared the majority of transcriptional regulatory motifs, suggesting convergent biological pathways. However, the transcriptional binding-factor motifs implicated in each condition implicated distinct biological mechanisms: neuronal JUN and FOS networks in autism vs. an inflammatory transcriptional network in dup15q microglia. This work provides a cell-type-specific analysis of how dup15q changes gene expression and chromatin accessibility in the human brain, and it finds evidence of marked cell-type-specific effects of this genetic driver. These findings have implications for guiding therapeutic development in dup15q syndrome, as well as understanding the functional effects of copy-number variants more broadly in neurodevelopmental disorders.

Features of the CD1 gene family in rodents and the uniqueness of the immune system of naked mole-rat

Cluster of Differentiation 1 (CD1) proteins are widely expressed throughout jawed vertebrates and present lipid antigens to specific CD1-restricted T lymphocytes. CD1 molecules play an important role in immune defense with the presence or absence of particular CD1 proteins frequently associated with the functional characteristics of the immune system. Here, we show the evolution of CD1 proteins in the Rodentia family and the diversity among its members. Based on the analysis of CD1 protein-coding regions in rodent genomes and the reconstruction of protein structures, we found that Heterocephalus glaber represents a unique member of the suborder Hystricomorpha with significant changes in protein sequences and structures of the CD1 family. Multiple lines of evidence point to the absence of CD1d and CD1e and probably a dysfunctional CD1b protein in Heterocephalus glaber. In addition, the impact of CD1d loss on the CD1d/Natural killer T (NKT) cell axis in the naked mole-rat and its potential implications for immune system function are discussed in detail.

GWAS breakthroughs: mapping the journey from one locus to 393 significant CAD associations.

Coronary artery disease (CAD) poses a substantial threat to global health, leading to significant morbidity and mortality worldwide. It has a significant genetic component that has been studied through genome-wide association studies (GWAS) over the past 17 years. These studies have made progress with larger sample sizes, diverse ancestral backgrounds, and the discovery of multiple genomic regions related to CAD risk. In this review, we provide a comprehensive overview of CAD GWAS, including information about the genetic makeup of the disease and the importance of ethnic diversity in these studies. We also discuss challenges of identifying causal genes and variants within GWAS loci with a focus on non-coding regions. Additionally, we highlight tissues and cell types relevant to CAD, and discuss clinical implications of GWAS findings including polygenic risk scores, sex-specific differences in CAD genetics, ethnical aspects of personalized interventions, and GWAS guided drug development.

Incorporating prior information in gene expression network-based cancer heterogeneity analysis.

Cancer is molecularly heterogeneous, with seemingly similar patients having different molecular landscapes and accordingly different clinical behaviors. In recent studies, gene expression networks have been shown as more effective/informative for cancer heterogeneity analysis than some simpler measures. Gene interconnections can be classified as "direct" and "indirect," where the latter can be caused by shared genomic regulators (such as transcription factors, microRNAs, and other regulatory molecules) and other mechanisms. It has been suggested that incorporating the regulators of gene expressions in network analysis and focusing on the direct interconnections can lead to a deeper understanding of the more essential gene interconnections. Such analysis can be seriously challenged by the large number of parameters (jointly caused by network analysis, incorporation of regulators, and heterogeneity) and often weak signals. To effectively tackle this problem, we propose incorporating prior information contained in the published literature. A key challenge is that such prior information can be partial or even wrong. We develop a two-step procedure that can flexibly accommodate different levels of prior information quality. Simulation demonstrates the effectiveness of the proposed approach and its superiority over relevant competitors. In the analysis of a breast cancer dataset, findings different from the alternatives are made, and the identified sample subgroups have important clinical differences.

Beyond cyanogenesis: Temperature gradients drive environmental adaptation in North American white clover (Trifolium repens L.).

Species that repeatedly evolve phenotypic clines across environmental gradients have been highlighted as ideal systems for characterizing the genomic basis of local environmental adaptation. However, few studies have assessed the importance of observed phenotypic clines for local adaptation: conspicuous traits that vary clinally may not necessarily be the most critical in determining local fitness. The present study was designed to fill this gap, using a plant species characterized by repeatedly evolved adaptive phenotypic clines. White clover is naturally polymorphic for its chemical defence cyanogenesis (HCN release with tissue damage); climate-associated cyanogenesis clines have evolved throughout its native and introduced range worldwide. We performed landscape genomic analyses on 415 wild genotypes from 43 locations spanning much of the North American species range to assess the relative importance of cyanogenesis loci vs. other genomic factors in local climatic adaptation. We find clear evidence of local adaptation, with temperature-related climatic variables best describing genome-wide differentiation between sampling locations. The same climatic variables are also strongly correlated with cyanogenesis frequencies and gene copy number variations (CNVs) at cyanogenesis loci. However, landscape genomic analyses indicate no significant contribution of cyanogenesis loci to local adaptation. Instead, several genomic regions containing promising candidate genes for plant response to seasonal cues are identified - some of which are shared with previously identified QTLs for locally adaptive fitness traits in North American white clover. Our findings suggest that local adaptation in white clover is likely determined primarily by genes controlling the timing of growth and flowering in response to local seasonal cues. More generally, this work suggests that caution is warranted when considering the importance of conspicuous phenotypic clines as primary determinants of local adaptation.

Detection and evaluation of parameters influencing the identification of heterozygous-enriched regions in Holstein cattle based on SNP chip or whole-genome sequence data

BackgroundA heterozygous-enriched region (HER) is a genomic region with high variability generated by factors such as balancing selection, introgression, and admixture processes. In this study, we evaluated the genomic background of HERs and the impact of different parameters (i.e., minimum number of SNPs in a HER, maximum distance between two consecutive SNPs, minimum length of a HER, maximum number of homozygous allowed in a HER) and scenarios [i.e., different SNP panel densities and whole-genome sequence (WGS)] on the detection of HERs. We also compared HERs characterized in Holstein cattle with those identified in Angus, Jersey, and Norwegian Red cattle using WGS data.ResultsThe parameters used for the identification of HERs significantly impact their detection. The maximum distance between two consecutive SNPs did not impact HERs detection as the same average of HERs (269.31 ± 787.00) was observed across scenarios. However, the minimum number of markers, maximum homozygous markers allowed inside a HER, and the minimum length size impacted HERs detection. For the minimum length size, the 10 Kb scenario showed the highest average number of HERs (1,364.69 ± 1,483.64). The number of HERs decreased as the minimum number of markers increased (621.31 ± 1,271.83 to 6.08 ± 21.94), and an opposite pattern was observed for the maximum homozygous markers allowed inside a HER (54.47 ± 195.51 to 494.89 ± 1,169.35). Forty-five HER islands located in 23 chromosomes with high Tajima’s D values and differential among the observed and estimated heterozygosity were detected in all evaluated scenarios, indicating their ability to potentially detect regions under balancing selection. In total, 3,440 markers and 28 genes previously related to fertility (e.g., TP63, ZSCAN23, NEK5, ARHGAP44), immunity (e.g., TP63, IGC, ARHGAP44), residual feed intake (e.g., MAYO9A), stress sensitivity (e.g., SERPINA6), and milk fat percentage (e.g., NOL4) were identified. When comparing HER islands among breeds, there were substantial overlaps between Holstein with Angus (95.3%), Jersey (94.3%), and Norwegian Red cattle (97.1%), indicating conserved HER across taurine breeds.ConclusionsThe detection of HERs varied according to the parameters used, but some HERs were consistently identified across all scenarios. Heterozygous genotypes observed across generations and breeds appear to be conserved in HERs. The results presented could serve as a guide for defining HERs detection parameters and further investigating their biological roles in future studies.

Distinct saliva DNA methylation profiles in relation to treatment outcome in youth with posttraumatic stress disorder

In youth with posttraumatic stress disorder (PTSD) non-response rates after treatment are often high. Epigenetic mechanisms such as DNA methylation (DNAm) have previously been linked to PTSD pathogenesis, additionally DNAm may affect response to (psychological) therapies. Besides investigating the direct link between DNAm and treatment response, it might be helpful to investigate the link between DNAm and previously associated biological mechanisms with treatment outcome. Thereby gaining a deeper molecular understanding of how psychotherapy (reflecting a change in the environment) relates to epigenetic changes and the adaptability of individuals. To date, limited research is done in clinical samples and no studies have been conducted in youth. Therefore we conducted a study in a Dutch cohort of youth with and without PTSD (n = 87, age 8–18 years). We examined the cross-sectional and longitudinal changes of saliva-based genome-wide DNA methylation (DNAm) levels, and salivary cortisol secretion. The last might reflect possible abbreviations on the hypothalamic–pituitary– adrenal (HPA) axis. The HPA-axis is previously linked to DNAm and the development and recovery of PTSD. Youth were treated with 8 sessions of either Eye Movement Reprocessing Therapy (EMDR) or Trauma Focused Cognitive behavioral Therapy (TF-CBT). Our epigenome wide approach showed distinct methylation between treatment responders and non-responders on C18orf63 gene post-treatment. This genomic region is related to the PAX5 gene, involved in neurodevelopment and inflammation response. Additionally, our targeted approach indicated that there were longitudinal DNAm changes in successfully treated youth at the CRHR2 gene. Methylation at this gene was further correlated with cortisol secretion pre- and post-treatment. Awaiting replication, findings of this first study in youth point to molecular pathways involved in stress response and neuroplasticity to be associated with treatment response.

Research Progress of Long Non-coding RNA-ZFAS1 in Malignant Tumors.

Long non-coding RNAs (lncRNAs), although incapable of encoding proteins, play crucial roles in multiple layers of gene expression regulation, epigenetic modifications, and post-transcriptional regulation. Zinc finger antisense 1 (ZFAS1), a lncRNA located in the 20q13 region of the human genome, exhibits dual functions as an oncogene or tumor suppressor in various human malignancies. ZFAS1 plays a crucial role in cancer progression, metastasis, invasion, apoptosis, cell cycle regulation, and drug resistance through complex molecular mechanisms. Additionally, ZFAS1 has a long half-life of over 16 h, demonstrating exceptional stability, and making it a potential biomarker. This review integrates recent studies on the role and molecular mechanisms of ZFAS1 in malignancies and summarizes its clinical significance. By summarizing the role of ZFAS1 in cancer, we aim to highlight its potential as an anti-cancer biomarker and therapeutic target.

Genome of the early spider-orchid Ophrys sphegodes provides insights into sexual deception and pollinator adaptation

Pollinator-driven evolution of floral traits is thought to be a major driver of angiosperm speciation and diversification. Ophrys orchids mimic female insects to lure male pollinators into pseudocopulation. This strategy, called sexual deception, is species-specific, thereby providing strong premating reproductive isolation. Identifying the genomic architecture underlying pollinator adaptation and speciation may shed light on the mechanisms of angiosperm diversification. Here, we report the 5.2 Gb chromosome-scale genome sequence of Ophrys sphegodes. We find evidence for transposable element expansion that preceded the radiation of the O. sphegodes group, and for gene duplication having contributed to the evolution of chemical mimicry. We report a highly differentiated genomic candidate region for pollinator-mediated evolution on chromosome 2. The Ophrys genome will prove useful for investigations into the repeated evolution of sexual deception, pollinator adaptation and the genomic architectures that facilitate evolutionary radiations.

Understanding natural genetic variation for nutritional quality in grain and identification of superior haplotypes in deepwater rice genotypes of Assam, India

Rice grown under deepwater ecosystem is considered to be natural farming and hence they are considered to be input efficient. Thus, to identify gene responsible for nutritional content under natural conditions, a genome-wide association study (GWAS)was performed. GWAS identified single nucleotide polymorphisms (SNPs) significantly associated with various nutritional quality traits such as Zn (mg/kg), Fe (mg/kg), Protein (%), Oil (%), Amylose (%), Starch (%), Phytic acid (%), Phenol (%) and TDF (%) in 184 deepwater rice accessions evaluated over 2 consecutive years. A total of 278 SNPs distributed across 12 chromosomes were found to be significantly associated with Zn, Oil and Phenol content. Among them, eight high confidence SNPs were significant and identified on chr1 (AX-95933712), chr7 (AX-95957036), and chr8 (AX-95965181) for Zn content. Similarly, on chr2 (AX-95945186), chr8 (AX-95964718), and chr11 (AX-95961099) have been found to be associated with Oil content and on chr3 (AX-95922121) and chr4 (AX-95963889) for Phenol content. Genomic regions of ± 220 kb flanking the three consistent lowest p value containing SNPs for each trait were considered for finding superior haplotypes. These SNPs showed significant phenotypic variations with different identified haplotype blocks. The allelic variations with phenotypes were considered to be superior haplotypes i.e., Block 1: Hap 1 (ACCC) for high Zn content, Block 2: Hap 1 (CT) for high Oil content, and Block 2: Hap 1(CGGG) for low Phenol content. The discovered superior haplotype with high nutritional content could be important for understanding the mechanisms involving nutrient use efficiency. Thus, the present study demonstrated that developing rice varieties with appropriate nutritional quality traits will be possible through the incorporation of such superior haplotypes in breeding programs.

Chromatin plasticity predetermines neuronal eligibility for memory trace formation.

Memories are encoded by sparse populations of neurons but how such sparsity arises remains largely unknown. We found that a neuron's eligibility to be recruited into the memory trace depends on its epigenetic state prior to encoding. Principal neurons in the mouse lateral amygdala display intrinsic chromatin plasticity, which when experimentally elevated favors neuronal allocation into the encoding ensemble. Such chromatin plasticity occurred at genomic regions underlying synaptic plasticity and was accompanied by increased neuronal excitability in single neurons in real time. Lastly, optogenetic silencing of the epigenetically altered neurons prevented memory expression, revealing a cell-autonomous relationship between chromatin plasticity and memory trace formation. These results identify the epigenetic state of a neuron as a key factor enabling information encoding.

Exploring cross-tissue DNA methylation patterns: blood–brain CpGs as potential neurodegenerative disease biomarkers

The difficulty of obtaining samples from certain human tissues has led to efforts to find accessible sources to analyze molecular markers derived from DNA. In this study, we look for DNA methylation patterns in blood samples and its association with the brain methylation pattern in neurodegenerative disorders. Using data from methylation databases, we selected 18,293 CpGs presenting correlated methylation levels between blood and brain (bb-CpGs) and compare their methylation level between blood samples from patients with neurodegenerative diseases (Alzheimer’s disease, Parkinson’s disease, Multiple Sclerosis, and X Fragile Syndrome) and healthy controls. Sixty-four bb-CpGs presented significant distinct methylation levels in patients, being: nine for Alzheimer’s disease, nine for Parkinson’s disease, 28 for Multiple Sclerosis, and 18 for Fragile X Syndrome. Similar differences in methylation pattern for the nine Alzheimer’s bb-CpGs was also observed when comparing brain tissue from patients vs. controls. The genomic regions of some of these 64 bb-CpGs are placed close to or inside genes previously associated with the respective condition. Our findings support the rationale of using blood DNA as a surrogate of brain tissue to analyze changes in CpG methylation level in patients with neurodegenerative diseases, opening the possibility for characterizing new biomarkers.

KPNA2 promotes the progression of gastric cancer by regulating the alternative splicing of related genes

RNA-binding proteins (RBPs) play critical roles in genome regulation. In this study, we explored the latent function of KPNA2, which is an essential member of the RBP family, in the regulation of alternative splicing (AS) in gastric cancer (GC). We analyzed the role of KPNA2 in regulating differential expression and AS via RNA sequencing (RNA-seq) and improved RNA immunoprecipitation sequencing (iRIP-seq). Clinical specimens were used to analyze the associations between KPNA2 expression and clinicopathological characteristics. CCK8 assays, transwell assays and wound healing assays were performed to explore the effect of KPNA2/WDR62 on GC cell progression. KPNA2 was shown to be highly expressed in GC cells and tissues and associated with lymph node metastases. KPNA2 promoted the proliferation, migration and invasion of GC cells and primarily regulated exon skipping, alternative 3's splice sites (A3SSs), alternative 5' splice sites (A5SSs), and cassette exons. We further revealed that KPNA2 participated in biological processes related to cell proliferation, and the immune response in GC via the regulation of transcription. In addition, KPNA2 preferentially bound to intron regions. Notably, KPNA2 regulated the A3SS AS mode of WDR62, and upregulation of WDR62 reversed the KPNA2 downregulation-induced inhibition of GC cell proliferation, migration and invasion. Finally, we discovered that the AS of immune-related molecules could be regulated by KPNA2. Overall, our results demonstrated for the first time that KPNA2 functions as an oncogenic splicing factor in GC that regulated the AS and differential expression of GC-related genes, and KPNA2 may be a potential target for GC treatment.

The evolutionary fate of Neanderthal DNA in 30,780 admixed genomes with recent African-like ancestry.

Following introgression, Neanderthal DNA was initially purged from non-African genomes, but the evolutionary fate of remaining introgressed DNA has not been explored yet. To fill this gap, we analyzed 30,780 admixed genomes with African-like ancestry from the All of Us research program, in which Neanderthal alleles encountered novel genetic backgrounds during the last 15 generations. Observed amounts of Neanderthal DNA approximately match expectations based on ancestry proportions, suggesting neutral evolution. Nevertheless, we identified genomic regions that have significantly less or more Neanderthal ancestry than expected and are associated with spermatogenesis, innate immunity, and other biological processes. We also identified three novel introgression desert-like regions in recently admixed genomes, whose genetic features are compatible with hybrid incompatibilities and intrinsic negative selection. Overall, we find that much of the remaining Neanderthal DNA in human genomes is not under strong selection, and complex evolutionary dynamics have shaped introgression landscapes in our species.

Genomic and Cell-Specific Regulation of Benzylisoquinoline Alkaloid Biosynthesis in Opium Poppy.

Opium poppy is a crop of great commercial value as a source of several opium alkaloids for the pharmaceutical industries including morphine, codeine, thebaine, noscapine and papaverine. Most enzymes involved in benzylisoquinoline alkaloids (BIAs) biosynthesis in opium poppy have been functionally characterized, and opium poppy currently serves as a model system to study BIA metabolism in plants. BIA biosynthesis in opium poppy involves two biosynthetic gene clusters associated respectively with the morphine and noscapine branches. Recent reports have shown that genes in the same cluster are co-expressed, suggesting they might also be co-regulated. However, the transcriptional regulation of opium poppy BIA biosynthesis is not well studied. Opium poppy BIA biosynthesis involves three cell types associated with the phloem system: companion cells, sieve elements and laticifers. The transcripts and enzymes associated with BIA biosynthesis are distributed across cell types, requiring the translocation of key enzymes and pathway intermediates between cell types. Together, these suggest that the regulation of BIA biosynthesis in opium poppy is multilayered and complex, involving biochemical, genomic, and physiological mechanisms. In this review, we highlight recent advances in genome sequencing and single cell and spatial transcriptomics with a focus on how these efforts can improve our understanding of the genomic and cell-specific regulation of BIA biosynthesis. Such knowledge is vital for opium poppy genetic improvement and metabolic engineering efforts targeting the modulation of alkaloid yield and composition.