Variants In Non-coding Regions Research Articles

Identification of Noncoding Functional Regulatory Variants of STIM1 Gene in Idiopathic Pulmonary Arterial Hypertension.

STIM1 (stromal interaction molecule 1) regulates store-operated calcium entry and is involved in pulmonary artery vasoconstriction and pulmonary artery smooth muscle cell proliferation, leading to pulmonary arterial hypertension (PAH). Bioinformatics analysis and a 2-stage matched case-control study were conducted to screen for noncoding variants that may potentially affect STIM1 transcriptional regulation in 242 patients with idiopathic PAH and 414 healthy controls. Luciferase reporter assay, real-time quantitative polymerase chain reaction, western blot, 5-ethynyl-2'-deoxyuridine (EdU) assay, and intracellular Ca2+ measurement were performed to study the mechanistic roles of those STIM1 noncoding variants in PAH. Five noncoding variants (rs3794050, rs7934581, rs3750996, rs1561876, and rs3750994) were identified and genotyped using Sanger sequencing. Rs3794050, rs7934581, and rs1561876 were associated with idiopathic PAH (recessive model, all P<0.05). Bioinformatics analysis showed that these 3 noncoding variants possibly affect the enhancer function of STIM1 or the microRNA (miRNA) binding to STIM1. Functional validation performed in HEK293 and pulmonary artery smooth muscle cells demonstrated that the noncoding variant rs1561876-G (STIM1 mutant) had significantly stronger transcriptional activity than the wild-type counterpart, rs1561876-A, by affecting the transcriptional regulatory function of both hsa-miRNA-3140-5p and hsa-miRNA-4766-5p. rs1561876-G enhanced intracellular Ca2+ signaling in human pulmonary artery smooth muscle cells secondary to calcium-sensing receptor activation and promoted proliferation of pulmonary artery smooth muscle cells under both normoxia and hypoxia conditions, suggesting a possible contribution to PAH development. The potential clinical implications of the 3 noncoding variants of STIM1, rs3794050, rs7934581, and rs1561876, are 2-fold, as they may help predict the risk and prognosis of idiopathic PAH and guide investigations on novel therapeutic pathway(s).

Single-nucleus multi-omics of Parkinson’s disease reveals a glutamatergic neuronal subtype susceptible to gene dysregulation via alteration of transcriptional networks

The genetic architecture of Parkinson’s disease (PD) is complex and multiple brain cell subtypes are involved in the neuropathological progression of the disease. Here we aimed to advance our understanding of PD genetic complexity at a cell subtype precision level. Using parallel single-nucleus (sn)RNA-seq and snATAC-seq analyses we simultaneously profiled the transcriptomic and chromatin accessibility landscapes in temporal cortex tissues from 12 PD compared to 12 control subjects at a granular single cell resolution. An integrative bioinformatic pipeline was developed and applied for the analyses of these snMulti-omics datasets. The results identified a subpopulation of cortical glutamatergic excitatory neurons with remarkably altered gene expression in PD, including differentially-expressed genes within PD risk loci identified in genome-wide association studies (GWAS). This was the only neuronal subtype showing significant and robust overexpression of SNCA. Further characterization of this neuronal-subpopulation showed upregulation of specific pathways related to axon guidance, neurite outgrowth and post-synaptic structure, and downregulated pathways involved in presynaptic organization and calcium response. Additionally, we characterized the roles of three molecular mechanisms in governing PD-associated cell subtype-specific dysregulation of gene expression: (1) changes in cis-regulatory element accessibility to transcriptional machinery; (2) changes in the abundance of master transcriptional regulators, including YY1, SP3, and KLF16; (3) candidate regulatory variants in high linkage disequilibrium with PD-GWAS genomic variants impacting transcription factor binding affinities. To our knowledge, this study is the first and the most comprehensive interrogation of the multi-omics landscape of PD at a cell-subtype resolution. Our findings provide new insights into a precise glutamatergic neuronal cell subtype, causal genes, and non-coding regulatory variants underlying the neuropathological progression of PD, paving the way for the development of cell- and gene-targeted therapeutics to halt disease progression as well as genetic biomarkers for early preclinical diagnosis.

Rare variation in non-coding regions with evolutionary signatures contributes to autism spectrum disorder risk.

Little is known about the role of non-coding regions in the etiology of autism spectrum disorder (ASD). We examined three classes of non-coding regions: human accelerated regions (HARs), which show signatures of positive selection in humans; experimentally validated neural VISTA enhancers (VEs); and conserved regions predicted to act as neural enhancers (CNEs). Targeted and whole-genome analysis of >16,600 samples and >4,900 ASD probands revealed that likely recessive, rare, inherited variants in HARs, VEs, and CNEs substantially contribute to ASD risk in probands whose parents share ancestry, which enriches for recessive contributions, but modestly contribute, if at all, in simplex family structures. We identified multiple patient variants in HARs near IL1RAPL1 and in VEs near OTX1 and SIM1 and showed that they change enhancer activity. Our results implicate both human-evolved and evolutionarily conserved non-coding regions in ASD risk and suggest potential mechanisms of how regulatory changes can modulate social behavior.

In silico generation and augmentation of regulatory variants from massively parallel reporter assay using conditional variational autoencoder.

Predicting the functional consequences of genetic variants in non-coding regions is a challenging problem. Massively parallel reporter assays (MPRAs), which are an in vitro high-throughput method, can simultaneously test thousands of variants by evaluating the existence of allele specific regulatory activity. Nevertheless, the identified labelled variants by MPRAs, which shows differential allelic regulatory effects on the gene expression are usually limited to the scale of hundreds, limiting their potential to be used as the training set for achieving a robust genome-wide prediction. To address the limitation, we propose a deep generative model, MpraVAE, to in silico generate and augment the training sample size of labelled variants. By benchmarking on several MPRA datasets, we demonstrate that MpraVAE significantly improves the prediction performance for MPRA regulatory variants compared to the baseline method, conventional data augmentation approaches as well as existing variant scoring methods. Taking autoimmune diseases as one example, we apply MpraVAE to perform a genome-wide prediction of regulatory variants and find that predicted regulatory variants are more enriched than background variants in enhancers, active histone marks, open chromatin regions in immune-related cell types, and chromatin states associated with promoter, enhancer activity and binding sites of cMyC and Pol II that regulate gene expression. Importantly, predicted regulatory variants are found to link immune-related genes by leveraging chromatin loop and accessible chromatin, demonstrating the importance of MpraVAE in genetic and gene discovery for complex traits.

Investigation of inherited noncoding genetic variation impacting the pharmacogenomics of childhood acute lymphoblastic leukemia treatment.

Defining genetic factors impacting chemotherapy failure can help to better predict response and identify drug resistance mechanisms. However, there is limited understanding of the contribution of inherited noncoding genetic variation on inter-individual differences in chemotherapy response in childhood acute lymphoblastic leukemia (ALL). Here we map inherited noncoding variants associated with treatment outcome and/or chemotherapeutic drug resistance to ALL cis-regulatory elements and investigate their gene regulatory potential and target gene connectivity using massively parallel reporter assays and three-dimensional chromatin looping assays, respectively. We identify 54 variants with transcriptional effects and high-confidence gene connectivity. Additionally, functional interrogation of the top variant, rs1247117, reveals changes in chromatin accessibility, PU.1 binding affinity and gene expression, and deletion of the genomic interval containing rs1247117 sensitizes cells to vincristine. Together, these data demonstrate that noncoding regulatory variants associated with diverse pharmacological traits harbor significant effects on allele-specific transcriptional activity and impact sensitivity to antileukemic agents.

Complete Chloroplast Genome of Krascheninnikovia ewersmanniana: Comparative and Phylogenetic Analysis.

Krascheninnikovia ewersmanniana is a dominant desert shrub in Xinjiang, China, with high economic and ecological value. However, molecular systematics research on K. ewersmanniana is lacking. To resolve the genetic composition of K. ewersmanniana within Amaranthaceae and its systematic relationship with related genera, we used a second-generation Illumina sequencing system to detect the chloroplast genome of K. ewersmanniana and analyze its assembly, annotation, and phylogenetics. Total length of the chloroplast genome of K. ewersmanniana reached 152,287 bp, with 84 protein-coding genes, 36 tRNAs, and eight rRNAs. Codon usage analysis showed the majority of codons ending with base A/U. Mononucleotide repeats were the most common (85.42%) of the four identified simple sequence repeats. A comparison with chloroplast genomes of six other Amaranthaceae species indicated contraction and expansion of the inverted repeat boundary region in K. ewersmanniana, with some genes (rps19, ndhF, ycf1) differing in length and distribution. Among the seven species, the variation in non-coding regions was greater. Phylogenetic analysis revealed Krascheninnikovia ceratoides, Dysphania ambrosioides, Dysphania pumilio, and Dysphania botrys to have a close monophyletic relationship. By sequencing the K. ewersmanniana chloroplast genome, this research resolves the relatedness among 35 Amaranthaceae species, providing molecular insights for germplasm utilization, and theoretical support for studying evolutionary relationships.

Prioritizing genomic variants pathogenicity via DNA, RNA, and protein-level features based on extreme gradient boosting.

Genetic diseases are mostly implicated with genetic variants, including missense, synonymous, non-sense, and copy number variants. These different kinds of variants are indicated to affect phenotypes in various ways from previous studies. It remains essential but challenging to understand the functional consequences of these genetic variants, especially the noncoding ones, due to the lack of corresponding annotations. While many computational methods have been proposed to identify the risk variants. Most of them have only curated DNA-level and protein-level annotations to predict the pathogenicity of the variants, and others have been restricted to missense variants exclusively. In this study, we have curated DNA-, RNA-, and protein-level features to discriminate disease-causing variants in both coding and noncoding regions, where the features of protein sequences and protein structures have been shown essential for analyzing missense variants in coding regions while the features related to RNA-splicing and RBP binding are significant for variants in noncoding regions and synonymous variants in coding regions. Through the integration of these features, we have formulated the Multi-level feature Genomic Variants Predictor (ML-GVP) using the gradient boosting tree. The method has been trained on more than 400,000 variants in the Sherloc-training set from the 6th critical assessment of genome interpretation with superior performance. The method is one of the two best-performing predictors on the blind test in the Sherloc assessment, and is further confirmed by another independent test dataset of de novo variants.

Abstract 7356: 3 UTR variations in esophageal squamous cell carcinoma in African Americans

Abstract Esophageal squamous cell carcinoma (ESCC) predominantly affects African Americans (AA). Esophageal Cancer (EC) consists of two major histological subtypes: esophageal adenocarcinoma and squamous cell carcinoma. ESCC ranks among the most lethal cancers and is particularly aggressive among AA, comprising about 75% of total ESCC patients. The 5-year survival rate for localized EC in AA patients is 25%; in contrast, the survival rate for EC patients of Caucasian origin is 48%. However, for metastatic disease, the survival rate decreases dramatically to 5% regardless of race. To conduct an extensive mutation screening, we performed whole exome sequencing (WES) in 10 paired AA ESCC and control tissues. Mutations included nonsynonymous mutations in multiple genes and copy number variations. Here, we report on further analysis of the WES data to identify 3 UTR noncoding region alterations, comprising the highest proportion of noncoding variations. Whole exome sequencing was performed using Agilent SureSelect XT Human All Exon V6 + UTR kit on matched normal-tumor samples from 10 AA patients (nine males and one female) with advanced-stage ESCC. The ages of patients ranged from 53 to 80 years. All patients, except for one, reported tobacco use and alcohol consumption. Our analysis yielded 5378 variants in noncoding regions in our cohort. The most variations were found in the 3UTR location, with 46% of all noncoding variations. In our dataset, variations in the 3UTR region of 48 genes were predicted by CScape to assert an oncogenic effect. The highest ten oncogenic predictive scores are for ETV1, FOXO3, MSI2, TCF7L2, KAT6A, FOXO1, IRF4, TRIM24, ELK4, and FGFR1 genes. Approximately 2500 genes having 3’UTR variations are over-presented in intracellular signaling and cell-cell communications, transcriptional regulation of pluripotent stem cells by OCT4, SOX2, NANOG, generic RNA polymerase II transcription, RHO GTPase cycle, and regulation of cytoskeletal remodeling. The combined effects of noncoding and coding region mutations may largely impact poor prognosis in AA with ESCC and yield targets for developing candidate therapeutics. Citation Format: Ceyda Erkizan, Robert Wadleigh. 3 UTR variations in esophageal squamous cell carcinoma in African Americans [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 7356.

Autozygome-guided exome-first study in a consanguineous cohort with early-onset retinal disease uncovers an isolated RIMS2 phenotype and a retina-enriched RIMS2 isoform.

Leber congenital amaurosis (LCA) and early-onset retinal degeneration (EORD) are inherited retinal diseases (IRD) characterized by early-onset vision impairment. Herein, we studied 15 Saudi families by whole exome sequencing (WES) and run-of-homozygosity (ROH) detection via AutoMap in 12/15 consanguineous families. This revealed (likely) pathogenic variants in 11/15 families (73%). A potential founder variant was found in RPGRIP1. Homozygous pathogenic variants were identified in known IRD genes (ATF6, CRB1, CABP4, RDH12, RIMS2, RPGRIP1, SPATA7). We established genotype-driven clinical reclassifications for ATF6, CABP4, and RIMS2. Specifically, we observed isolated IRD in the individual with the novel RIMS2 variant, and we found a retina-enriched RIMS2 isoform conserved but not annotated in mouse. The latter illustrates potential different phenotypic consequences of pathogenic variants depending on the particular tissue/cell-type specific isoforms they affect. Lastly, a compound heterozygous genotype in GUCY2D in one non-consanguineous family was demonstrated, and homozygous variants in novel candidate genes ATG2B and RUFY3 were found in the two remaining consanguineous families. Reporting these genes will allow to validate them in other IRD cohorts. Finally, the missing heritability of the two unsolved IRD cases may be attributed to variants in non-coding regions or structural variants that remained undetected, warranting future WGS studies.

Upstream open reading frame-introducing variants in patients with primary familial brain calcification.

More than 50% of patients with primary familial brain calcification (PFBC), a rare neurological disorder, remain genetically unexplained. While some causative genes are yet to be identified, variants in non-coding regions of known genes may represent a source of missed diagnoses. We hypothesized that 5'-Untranslated Region (UTR) variants introducing an AUG codon may initiate mRNA translation and result in a loss of function in some of the PFBC genes. After reannotation of exome sequencing data of 113 unrelated PFBC probands, we identified two upstream AUG-introducing variants in the 5'UTR of PDGFB. One, NM_002608.4:c.-373C>G, segregated with PFBC in the family. It was predicted to create an upstream open reading frame (ORF). The other one, NM_002608.4:c.-318C>T, was found in a simplex case. It was predicted to result in an ORF overlapping the natural ORF with a frameshift. In a GFP reporter assay, both variants were associated with a dramatic decrease in GFP levels, and, after restoring the reading frame with the GFP sequence, the c.-318C>T variant was associated with a strong initiation of translation as measured by western blotting. Overall, we found upstream AUG-introducing variants in the 5'UTR of PDGFB in 2/113 (1.7%) undiagnosed PFBC cases. Such variants thus represent a source of putative pathogenic variants.

NCAD v1.0: A database for non-coding variant annotation and interpretation

The application of whole genome sequencing is expanding in clinical diagnostics across various genetic disorders, and the significance of non-coding variants in penetrant diseases is increasingly being demonstrated. Therefore, it is urgent to improve the diagnostic yield by exploring the pathogenic mechanisms of variants in non-coding regions. However, the interpretation of non-coding variants remains a significant challenge, due to the complex functional regulatory mechanisms of non-coding regions and the current limitations of available databases and tools. Hence, we develop the non-coding variant annotation database (NCAD, http://www.ncawdb.net/), encompassing comprehensive insights into 665,679,194 variants, regulatory elements, and element interaction details. Integrating data from 96 sources, spanning both GRCh37 and GRCh38 versions, NCAD v1.0 provides vital information to support the genetic diagnosis of non-coding variants, including allele frequencies of 12 diverse populations, with a particular focus on the population frequency information for 230,235,698 variants in 20,964 Chinese individuals. Moreover, it offers prediction scores for variant functionality, five categories of regulatory elements, and four types of non-coding RNAs. With its rich data and comprehensive coverage, NCAD serves as a valuable platform, empowering researchers and clinicians with profound insights into non-coding regulatory mechanisms while facilitating the interpretation of non-coding variants.

Functional analysis of recurrent CDC20 promoter variants in human melanoma

Small nucleotide variants in non-coding regions of the genome can alter transcriptional regulation, leading to changes in gene expression which can activate oncogenic gene regulatory networks. Melanoma is heavily burdened by non-coding variants, representing over 99% of total genetic variation, including the well-characterized TERT promoter mutation. However, the compendium of regulatory non-coding variants is likely still functionally under-characterized. We developed a pipeline to identify hotspots, i.e. recurrently mutated regions, in melanoma containing putatively functional non-coding somatic variants that are located within predicted melanoma-specific regulatory regions. We identified hundreds of statistically significant hotspots, including the hotspot containing the TERT promoter variants, and focused on a hotspot in the promoter of CDC20. We found that variants in the promoter of CDC20, which putatively disrupt an ETS motif, lead to lower transcriptional activity in reporter assays. Using CRISPR/Cas9, we generated an indel in the CDC20 promoter in human A375 melanoma cell lines and observed decreased expression of CDC20, changes in migration capabilities, increased growth of xenografts, and an altered transcriptional state previously associated with a more proliferative and less migratory state. Overall, our analysis prioritized several recurrent functional non-coding variants that, through downregulation of CDC20, led to perturbation of key melanoma phenotypes.

LncRNA involvement in immune-related diseases - from SNP association to implication in pathogenesis and therapeutic potential

Development of new high throughput array-based techniques and, more recently, next-generation sequencing (NGS) technologies have revolutionized our capability to accurately characterize single nucleotide polymorphisms (SNPs) throughout the genome. These advances have facilitated large-scale genome-wide association studies (GWAS), which have served as fundamental elements in establishing links between SNPs and the susceptibility to several complex diseases, including those related to the immune system. Nevertheless, the molecular mechanisms underlying the development of most of these disorders are still poorly defined. Decoding the functionality of SNPs becomes increasingly challenging due to the predominant presence of these risk variants in non-coding regions of the genome. Among them, long non-coding RNAs (lncRNAs) are enriched in disease-associated SNPs. lncRNAs are involved in governing the control of gene expression both during transcription and at the post-transcriptional level. The existence of SNPs within the sequences of lncRNAs has the potential to alter their expression, structure, or function. This, in turn, can influence their regulatory roles and consequently contribute to the onset or progression of various diseases. In this review, we describe the implication of SNPs located in lncRNAs in the development of different immune-related diseases and highlight the potential of these molecules in the development of emerging RNA-based therapies.

Functional assay for assessment of pathogenicity of BAP1 variants.

Pathogenic germline variants in BRCA1-Associated Protein 1 (BAP1) cause BAP1 tumor predisposition syndrome (BAP1-TPDS). Carriers run especially a risk of uveal (UM) and cutaneous melanoma, malignant mesothelioma, and clear cell renal carcinoma. Approximately half of increasingly reported BAP1 variants lack accurate classification. Correct interpretation of pathogenicity can improve prognosis of the patients through tumor screening with better understanding of BAP1-TPDS. We edited five rare BAP1 variants with differing functional characteristics identified from patients with UM in HAP1 cells using CRISPR-Cas9 and assayed their effect on cell adhesion/spreading (at 4h) and proliferation (at 48h), measured as cell index (CI), using xCELLigence real-time analysis system. In BAP1 knockout HAP1 cultures, cell number was half of wild type (WT) cultures at 48h (p = 0.00021), reaching confluence later, and CI was 78% reduced (p < 0.0001). BAP1-TPDS-associated null variants c.67+1G>T and c.1780_1781insT, and a likely pathogenic missense variant c.281A>G reduced adhesion (all p ≤ 0.015) and proliferation by 74%-83% (all p ≤ 0.032). Another likely pathogenic missense variant c.680G>A reduced both by at least 50% (all p ≤ 0.032), whereas cells edited with likely benign one c.1526C>T grew similarly to WT. BAP1 is essential for optimal fitness of HAP1 cells. Pathogenic and likely pathogenic BAP1 variants reduced cell fitness, reflected in adhesion/spreading and proliferation properties. Further, moderate effects were quantifiable. Variant modelling in HAP1 with CRISPR-Cas9 enabled functional analysis of coding and non-coding region variants in an endogenous expression system.

Polygenic modelling and machine learning approaches in pharmacogenomics: importance in downstream analysis of genome-wide association study data.

Genome-wide association studies (GWAS) have identified genetic variations associated with adverse drug effects in pharmacogenomics (PGx) research. Yet, interpreting the biological implications of these associations remains a challenge. This review highlights two promising post-GWAS methods for PGx. First, we discuss the polygenic architecture of the PGx traits, especially for drug-induced liver injury (DILI). Their experimental modelling using multiple donors' human primary hepatocytes and human liver organoids demonstrated the polygenic architecture of DILI susceptibility and found biological vulnerability in genetically high-risk tissue donors. Second, we discuss the challenges of interpreting the roles of variants in non-coding regions. Beyond methods involving expression quantitative trait locus analysis and massively parallel reporter assays, we suggest the use of in silico mutagenesis through machine learning methods to understand the roles of variants in transcriptional regulation. This review underscores the importance of these post-GWAS methods in providing critical insights into PGx, potentially facilitating drug development and personalized treatment. (149 words).

Comprehensive genetic analysis reveals the mutational landscape of ABCA4-associated retinal dystrophy in a Chinese cohort

PurposeTo depict the variant profiles of the ABCA4 gene in a large Chinese cohort of patients with ABCA4-associated retinal dystrophy (ABCA4-RD). MethodsWe recruited 290 unrelated Chinese patients with ABCA4-RD and did ABCA4 mutational screening by a combination of Sanger sequencing, targeted exome sequencing, entire ABCA4 locus sequencing, and whole genome sequencing (WGS). The pathogenicity of variants was assessed using in silico tools or in vitro splicing assays following the American College of Medical Genetics and Genomics guidelines. ResultsTwo hundred sixty-eight distinct pathogenic variants were identified, and 57 were novel. In 580 alleles, 22 noncoding region variants outside canonical splice sites and 4 structural variations were found in 44 alleles accounting for 7.6% of all alleles. Bioinformatics analysis showed the complex mechanism of aberrant splicing productsnatural splice site disruption, branch point destruction, and cryptic splice site activation. Correspondingly, minigene assays validated the various abnormal splicing products, including exon skipping, exon elongation, partial exon deletion, and pseudoexon insertion. WGS identified the first inversion variation in ABCA4. ConclusionsThis study systematically depicted the variant profiles of ABCA4 and revealed the missing alleles of patients with ABCA4-RD in a large Chinese cohort. Our findings demonstrated the complexity of molecular diagnosis of Mendelian diseases and the efficiency of WGS for detecting structural variants.

Complete chloroplast of four Sanicula taxa (Apiaceae) endemic to China: lights into genome structure, comparative analysis, and phylogenetic relationships

BackgroundThe genus Sanicula comprises ca. 45 taxa, widely distributed from East Asia to North America, which is a taxonomically difficult genus with high medicinal value in Apiaceae. The systematic classification of the genus has been controversial for a long time due to varied characters in key morphological traits. China is one of the most important distributed centers, with ca. 18 species and two varieties. At present, chloroplast genomes are generally considered to be conservative and play an important role in evolutionary relationship study. To investigate the plastome evolution and phylogenetic relationships of Chinese Sanicula, we comprehensively analyzed the structural characteristics of 13 Chinese Sanicula chloroplasts and reconstructed their phylogenetic relationships.ResultsIn present study, four newly complete chloroplast genome of Sanicula taxa by using Illumina sequencing were reported, with the typical quadripartite structure and 155,396–155,757 bp in size. They encoded 126 genes, including 86 protein-coding genes, 32 tRNA genes and 8 rRNA genes. Genome structure, distributions of SDRs and SSRs, gene content, among Sanicula taxa, were similar. The nineteen intergenic spacers regions, including atpH-atpI, ndhC-trnM, petB-petD, petD-rpoA, petN-psbM, psaJ-rpl33, rbcL-accD, rpoB-trnC, rps16-trnQ, trnE-psbD, trnF-ndhJ, trnH-psbA, trnN-ndhF, trnS-psbZ, trnS-trnR, trnT-trnF, trnV-rps12, ycf3-trnS and ycf4-cemA, and one coding region (ycf1 gene) were the most variable. Results of maximum likelihood analysis based on 79 unique coding genes of 13 Chinese Sanicula samples and two Eryngium (Apiaceae-Saniculoideae) species as outgroup taxa revealed that they divided into four subclades belonged to two clades, and one subclade was consistent with previously traditional Sanicula section of its system. The current classification based on morphology at sect. Sanicla and Sect. Tuberculatae in Chinese Sanicula was not supported by analysis of cp genome phylogeny.ConclusionsThe chloroplast genome structure of Sanicula was similar to other angiosperms and possessed the typical quadripartite structure with the conserved genome arrangement and gene features. However, their size varied owing to expansion/contraction of IR/SC boundaries. The variation of non-coding regions was larger than coding regions of the chloroplast genome. Phylogenetic analysis within these Chinese Sanicula were determined using the 79 unique coding genes. These results could provide important data for systematic, phylogenomic and evolutionary research in the genus for the future studies.

Finding Needles in the Haystack: Strategies for Uncovering Noncoding Regulatory Variants.

Despite accumulating evidence implicating noncoding variants in human diseases, unraveling their functionality remains a significant challenge. Systematic annotations of the regulatory landscape and the growth of sequence variant data sets have fueled the development of tools and methods to identify causal noncoding variants and evaluate their regulatory effects. Here, we review the latest advances in the field and discuss potential future research avenues to gain a more in-depth understanding of noncoding regulatory variants.

Complex Regulation of Gamma-Hemolysin Expression Impacts Staphylococcus aureus Virulence.

Staphylococcus aureus gamma-hemolysin CB (HlgCB) is a core-genome-encoded pore-forming toxin that targets the C5a receptor, similar to the phage-encoded Panton-Valentine leucocidin (PVL). Absolute quantification by mass spectrometry of HlgCB in 39 community-acquired pneumonia (CAP) isolates showed considerable variations in the HlgC and HlgB yields between isolates. Moreover, although HlgC and HlgB are encoded on a single operon, their levels were dissociated in 10% of the clinical strains studied. To decipher the molecular basis for the variation in hlgCB expression and protein production among strains, different regulation levels were analyzed in representative clinical isolates and reference strains. Both the HlgCB level and the HlgC/HlgB ratio were found to depend on hlgC promoter activity and mRNA processing and translation. Strikingly, only one single nucleotide polymorphism (SNP) in the 5' untranslated region (UTR) of hlgCB mRNA strongly impaired hlgC translation in the USA300 strain, leading to a strong decrease in the level of HlgC but not in HlgB. Finally, we found that high levels of HlgCB synthesis led to mortality in a rabbit model of pneumonia, correlated with the implication of the role of HlgCB in severe S. aureus CAP. Taken together, this work illustrates the complexity of virulence factor expression in clinical strains and demonstrates a butterfly effect where subtle genomic variations have a major impact on phenotype and virulence. IMPORTANCE S. aureus virulence in pneumonia results in its ability to produce several virulence factors, including the leucocidin PVL. Here, we demonstrate that HlgCB, another leucocidin, which targets the same receptors as PVL, highly contributes to S. aureus virulence in pvl-negative strains. In addition, considerable variations in HlgCB quantities are observed among clinical isolates from patients with CAP. Biomolecular analyses have revealed that a few SNPs in the promoter sequences and only one SNP in the 5' UTR of hlgCB mRNA induce the differential expression of hlgCB, drastically impacting hlgC mRNA translation. This work illustrates the subtlety of regulatory mechanisms in bacteria, especially the sometimes major effects on phenotypes of single nucleotide variation in noncoding regions.

Blood cell traits’ GWAS loci colocalization with variation in PU.1 genomic occupancy prioritizes causal noncoding regulatory variants

Genome-wide association studies (GWASs) have uncovered numerous trait-associated loci across the human genome, most of which are located in noncoding regions, making interpretation difficult. Moreover, causal variants are hard to statistically fine-map at many loci because of widespread linkage disequilibrium. To address this challenge, we present a strategy utilizing transcription factor (TF) binding quantitative trait loci (bQTLs) for colocalization analysis to identify trait associations likely mediated by TF occupancy variation and to pinpoint likely causal variants using motif scores. We applied this approach to PU.1 bQTLs in lymphoblastoid cell lines and blood cell trait GWAS data. Colocalization analysis revealed 69 blood cell trait GWAS loci putatively driven by PU.1 occupancy variation. We nominate PU.1 motif-altering variants as the likely shared causal variants at 51 loci. Such integration of TF bQTL data with other GWAS data may reveal transcriptional regulatory mechanisms and causal noncoding variants underlying additional complex traits.