Variants In Region Research Articles

Datasets-Based IMPDH1 Revisited: Heterozygous Missense Variants for Dominant Retinitis Pigmentosa While Truncation Variants Are Likely Non-Pathogenic

Purpose Heterozygous variants of IMPDH1 are associated with autosomal dominant retinitis pigmentosa (adRP). The current study aims to investigate the characteristics of the adRP-associated variants. Methods IMPDH1 variants from our exome sequencing dataset were retrieved and systemically evaluated through multiple online prediction tools, comparative genomics (in-house dataset, HGMD, and gnomAD), and phenotypic association. Potential pathogenic variants (PPVs) were further confirmed by Sanger sequencing and segregation analysis. Results In total, seven heterozygous PPVs (six missenses and one inframe) were identified in 10 families with RP, in which six of the seven might be classified as pathogenic or likely pathogenic while one others as variants of uncertain significance. IMPDH1 variants contributed to 0.7% (10/1519) of RP families in our cohort, ranking the top four genes implicated in adRP. These adRP-associated variants were located in exons 8–10, a region within or downstream of the CBS domain. All these variants were predicted to be damaged by at least three of the six online prediction tools. Two truncation variants were considered non-pathogenic. Hitherto, 41 heterozygous variants of IMPDH1 were detected in 110 families in published literature, including 33 missenses, two inframes, and six truncations (including a synonymous variant affecting splicing). Of the 35 missense and inframe variants, most were clustered in exons 8–10 (77.1%, 27/35), including 18 (51.4%, 18/35) in exon 10 accounting for 70.9% (78/110) of the families. However, truncation variants were enriched in the general population with a pLI value of 0 (tolerated), and the reported variants in patients with RP did not cluster in specific region. Conclusions Our data together with comprehensive analysis of existing datasets suggest that causative variants of IMPDH1 are usually missense and mostly clustered in exons 8–10. Conversely, most missense variants outside this region and truncation variants should be interpreted with great care in clinical gene test.

A regulatory variant impacting TBX1 expression contributes to basicranial morphology in Homo sapiens

Changes in gene regulatory elements play critical roles in human phenotypic divergence. However, identifying the base-pair changes responsible for the distinctive morphology of Homo sapiens remains challenging. Here, we report a noncoding single-nucleotide polymorphism (SNP), rs41298798, as a potential causal variant contributing to the morphology of the skull base and vertebral structures found in Homo sapiens. Screening for differentially regulated genes between Homo sapiens and extinct relatives revealed 13 candidate genes associated with basicranial development, with TBX1, implicated in DiGeorge syndrome, playing a pivotal role. Epigenetic markers and in silico analyses prioritized rs41298798 within a TBX1 intron for functional validation. CRISPR editing revealed that the 41-base-pair region surrounding rs41298798 modulates gene expression at 22q11.21. The derived allele of rs41298798 acts as an allele-specific enhancer mediated by E2F1, resulting in increased TBX1 expression levels compared to the ancestral allele. Tbx1-knockout mice exhibited skull base and vertebral abnormalities similar to those seen in DiGeorge syndrome. Phenotypic differences associated with TBX1 deficiency are observed between Homo sapiens and Neanderthals (Homo neanderthalensis). In conclusion, the regulatory divergence of TBX1 contributes to the formation of skull base and vertebral structures found in Homo sapiens.

The quantitative genetics of gene expression in Mimulus guttatus.

Gene expression can be influenced by genetic variants that are closely linked to the expressed gene (cis eQTLs) and variants in other parts of the genome (trans eQTLs). We created a multiparental mapping population by sampling genotypes from a single natural population of Mimulus guttatus and scored gene expression in the leaves of 1,588 plants. We find that nearly every measured gene exhibits cis regulatory variation (91% have FDR < 0.05). cis eQTLs are usually allelic series with three or more functionally distinct alleles. The cis locus explains about two thirds of the standing genetic variance (on average) but varies among genes and tends to be greatest when there is high indel variation in the upstream regulatory region and high nucleotide diversity in the coding sequence. Despite mapping over 10,000 trans eQTL / affected gene pairs, most of the genetic variance generated by trans acting loci remains unexplained. This implies a large reservoir of trans acting genes with subtle or diffuse effects. Mapped trans eQTLs show lower allelic diversity but much higher genetic dominance than cis eQTLs. Several analyses also indicate that trans eQTLs make a substantial contribution to the genetic correlations in expression among different genes. They may thus be essential determinants of "gene expression modules," which has important implications for the evolution of gene expression and how it is studied by geneticists.

Coordination among frequent genetic variants imparts substance use susceptibility and pathogenesis.

Determining the key genetic variants is a crucial step to comprehensively understand substance use disorders (SUDs). In this study, utilizing whole exome sequences of five multi-generational pedigrees with SUDs, we used an integrative omics-based approach to uncover candidate genetic variants that impart susceptibility to SUDs and influence addition traits. We identified several SNPs and rare, protein-function altering variants in genes, GRIA3, NCOR1, and SHANK1; compound heterozygous variants in LNPEP, LRP1, and TBX2, that play a significant role in the neurotransmitter-neuropeptide axis, specifically in the dopaminergic circuits. We also noted a greater frequency of heterozygous and recessive variants in genes involved in the structural and functional integrity of synapse receptors, CHRNA4, CNR2, GABBR1, DRD4, NPAS4, ADH1B, ADH1C, OPRM1, and GABBR2. Variant analysis in upstream promoter regions revealed regulatory variants in NEK9, PRRX1, PRPF4B, CELA2A, RABGEF1, and CRBN, crucial for dopamine regulation. Using family-and pedigree-based data, we identified heterozygous recessive alleles in LNPEP, LRP1 (4 frameshift deletions), and TBX2 (2 frameshift deletions) linked to SUDs. GWAS overlap identified several SNPs associated with SUD susceptibility, including rs324420 and rs1229984. Furthermore, miRNA variant analysis revealed notable variants in mir-548 U and mir-532. Pathway studies identified the presence of extensive coordination among these genetic variants to impart substance use susceptibility and pathogenesis. This study identified variants that were found to be overrepresented among genes of dopaminergic circuits participating in the neurotransmitter-neuropeptide axis, suggesting pleiotropic influences in the development and sustenance of chronic substance use. The presence of a diverse set of haploinsufficient variants in varying frequencies demonstrates the existence of extraordinary coordination among them in attributing risk and modulating severity to SUDs.

Allele-Specific Regulation of the Candidate Autism Liability Gene RAI1 by the Enhancer Variant rs4925102 (C/G).

Retinoic acid-induced 1 (RAI1) is a dosage-sensitive gene that causes autistic phenotypes when deleted or duplicated. Observations from clinical cases and animal models also suggest that changes of RAI1 expression levels contribute to autism. Previously, we used a bioinformatic approach to identify several single nucleotide polymorphisms (SNPs) located within the 5'-region of RAI1 that correlate with RAI1 mRNA expression in the human brain. In particular, the SNP rs4925102 was identified as a candidate cis-acting regulatory variant, the genotype of which may affect the binding of transcription factors that influence RAI1 mRNA expression. In this study, we provide experimental evidence based on reporter gene, chromatin immunoprecipitation (ChIP), and chromatin conformation capture (3C) assays that rs4925102 regulates RAI1 mRNA expression in an allele-specific manner in human cell lines, including the neuroblastoma-derived cell line SH-SY5Y. We also describe a statistically significant association between rs4925102 genotype and autism spectrum disorder (ASD) diagnosis in a case-control study and near-statistically significant association in an Autism Genome Project (AGP) transmission disequilibrium (TDT) study using Caucasian subjects.

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.

Exploring non-coding genetic variability in ACE2: Functional annotation and in vitro validation of regulatory variants

The surge in human whole-genome sequencing data has facilitated the study of non-coding region variations, yet understanding their biological significance remains a challenge. We used a computational workflow to assess the regulatory potential of non-coding variants, with a particular focus on the Angiotensin Converting Enzyme 2 (ACE2) gene. This gene is crucial in physiological processes and serves as the entry point for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus causing coronavirus disease 19 (COVID-19). In our analysis, using data from the gnomAD population database and functional annotation, we identified 17 significant Single Nucleotide Variants (SNVs) in ACE2, particularly in its enhancers, promoters, and 3′ untranslated regions (UTRs). We found preliminary evidence supporting the regulatory impact of some of these variants on ACE2 expression. Our detailed examination of two SNVs, rs147718775 and rs140394675, in the ACE2 promoter revealed that these co-occurring SNVs, when mutated, significantly enhance promoter activity, suggesting a possible increase in specific ACE2 isoform expression. This method proves effective in identifying and interpreting impactful non-coding variants, aiding in further studies and enhancing understanding of molecular bases of monogenic and complex traits.

Detection of Novel BEST1 Variations in Autosomal Recessive Bestrophinopathy Using Third-generation Sequencing.

Autosomal recessive bestrophinopathy (ARB), a retinal degenerative disease, is characterized by central visual loss, yellowish multifocal diffuse subretinal deposits, and a dramatic decrease in the light peak on electrooculogram. The potential pathogenic mechanism involves mutations in the BEST1 gene, which encodes Ca2+-activated Cl- channels in the retinal pigment epithelium (RPE), resulting in degeneration of RPE and photoreceptor. In this study, the complete clinical characteristics of two Chinese ARB families were summarized. Pacific Biosciences (PacBio) single-molecule real-time (SMRT) sequencing was performed on the probands to screen for disease-causing gene mutations, and Sanger sequencing was applied to validate variants in the patients and their family members. Two novel mutations, c.202T>C (chr11:61722628, p.Y68H) and c.867+97G>A, in the BEST1 gene were identified in the two Chinese ARB families. The novel missense mutation BEST1 c.202T>C (p.Y68H) resulted in the substitution of tyrosine with histidine in the N-terminal region of transmembrane domain 2 of bestrophin-1. Another novel variant, BEST1 c.867+97G>A (chr11:61725867), located in intron 7, might be considered a regulatory variant that changes allele-specific binding affinity based on motifs of important transcriptional regulators. Our findings represent the first use of third-generation sequencing (TGS) to identify novel BEST1 mutations in patients with ARB, indicating that TGS can be a more accurate and efficient tool for identifying mutations in specific genes. The novel variants identified further broaden the mutation spectrum of BEST1 in the Chinese population.

Insights from the SNP analysis of TYMP gene linking MNGIE.

TYMP gene, which codes for thymidine phosphorylase (TP) is also known as platelet-derived endothelial cell growth factor (PD-ECGF). TP plays crucial roles in nucleotide metabolism and angiogenesis. Mutations in the TYMP gene can lead to Mitochondrial Neurogastrointestinal Encephalopathy (MNGIE) syndrome, a rare genetic disorder. Our main objective was to evaluate the impact of detrimental non-synonymous single nucleotide polymorphisms (nsSNPs) on TP protein structure and predict harmful variants in untranslated regions (UTR). We employed a combination of predictive algorithms to identify nsSNPs with potential deleterious effects, followed by molecular modeling analysis to understand their effects on protein structure and function. Using 13 algorithms, we identified 119 potentially deleterious nsSNPs, with 82 located in highly conserved regions. Of these, 53 nsSNPs were functional and exposed, while 79 nsSNPs reduced TP protein stability. Further analysis of 18 nsSNPs through 3D protein structure analysis revealed alterations in amino acid interactions, indicating their potential impact on protein function. This will help in the development of faster and more efficient genetic tests for detecting TYMP gene mutations.

Omics-based construction of regulatory variants can be applied to help decipher pig liver-related traits

Genetic variants can influence complex traits by altering gene expression through changes to regulatory elements. However, the genetic variants that affect the activity of regulatory elements in pigs are largely unknown, and the extent to which these variants influence gene expression and contribute to the understanding of complex phenotypes remains unclear. Here, we annotate 90,991 high-quality regulatory elements using acetylation of histone H3 on lysine 27 (H3K27ac) ChIP-seq of 292 pig livers. Combined with genome resequencing and RNA-seq data, we identify 28,425 H3K27ac quantitative trait loci (acQTLs) and 12,250 expression quantitative trait loci (eQTLs). Through the allelic imbalance analysis, we validate two causative acQTL variants in independent datasets. We observe substantial sharing of genetic controls between gene expression and H3K27ac, particularly within promoters. We infer that 46% of H3K27ac exhibit a concomitant rather than causative relationship with gene expression. By integrating GWAS, eQTLs, acQTLs, and transcription factor binding prediction, we further demonstrate their application, through metabolites dulcitol, phosphatidylcholine (PC) (16:0/16:0) and published phenotypes, in identifying likely causal variants and genes, and discovering sub-threshold GWAS loci. We provide insight into the relationship between regulatory elements and gene expression, and the genetic foundation for dissecting the molecular mechanism of phenotypes.

Mitochondrial and Nuclear DNA Variants in Amyotrophic Lateral Sclerosis: Enrichment in the Mitochondrial Control Region and Sirtuin Pathway Genes in Spinal Cord Tissue.

Amyotrophic Lateral Sclerosis (ALS) is a progressive disease with prevalent mitochondrial dysfunctions affecting both upper and lower motor neurons in the motor cortex, brainstem, and spinal cord. Despite mitochondria having their own genome (mtDNA), in humans, most mitochondrial genes are encoded by the nuclear genome (nDNA). Our study aimed to simultaneously screen for nDNA and mtDNA genomes to assess for specific variant enrichment in ALS compared to control tissues. Here, we analysed whole exome (WES) and whole genome (WGS) sequencing data from spinal cord tissues, respectively, of 6 and 12 human donors. A total of 31,257 and 301,241 variants in nuclear-encoded mitochondrial genes were identified from WES and WGS, respectively, while mtDNA reads accounted for 73 and 332 variants. Despite technical differences, both datasets consistently revealed a specific enrichment of variants in the mitochondrial Control Region (CR) and in several of these genes directly associated with mitochondrial dynamics or with Sirtuin pathway genes within ALS tissues. Overall, our data support the hypothesis of a variant burden in specific genes, highlighting potential actionable targets for therapeutic interventions in ALS.

The Association between Anatomical Variants of Musculoskeletal Structures and Nerve Compressions of the Lower Limb: A Systematic Review and Meta-Analysis.

Objective: The aim of this study was to describe the main anatomical variants and morphofunctional alterations in the lower limb that compress surrounding nervous structures in the gluteal region, thigh region, and leg and foot region. Methods: We searched the Medline, Scopus, Web of Science, Google Scholar, CINAHL, and LILACS databases from their inception up to October 2023. An assurance tool for anatomical studies (AQUA) was used to evaluate methodological quality, and the Joanna Briggs Institute assessment tool for case reports was also used. Forest plots were generated to assess the prevalence of variants of the gluteal region, thigh, and leg. Results: According to the forest plot of the gluteal region, the prevalence was 0.18 (0.14-0.23), with a heterogeneity of 93.52%. For the thigh region, the forest plot presented a prevalence of 0.10 (0.03-0.17) and a heterogeneity of 91.18%. The forest plot of the leg region was based on seven studies, which presented a prevalence of 0.01 (0.01-0.01) and a heterogeneity of 96.18%. Conclusions: This review and meta-analysis showed that, in studies that analyzed nerve compressions, the prevalence was low in the thigh and leg regions, while in the gluteal region, it was slightly higher. This is mainly due to the PM region and its different variants. We believe that it is important to analyze all the variant regions defined in this study and that surgeons treating the lower limb should be attentive to these possible scenarios so that they can anticipate possible surgical situations and thus avoid surgical complications.

A noncoding regulatory variant in IKZF1 increases acute lymphoblastic leukemia risk in Hispanic/Latino children

Hispanic/Latino children have the highest risk of acute lymphoblastic leukemia (ALL) in the US compared to other racial/ethnic groups, yet the basis of this remains incompletely understood. Through genetic fine-mapping analyses, we identified a new independent childhood ALL risk signal near IKZF1 in self-reported Hispanic/Latino individuals, but not in non-Hispanic White individuals, with an effect size of ∼1.44 (95% confidence interval= 1.33-1.55) and a risk allele frequency of ∼18% in Hispanic/Latino populations and <0.5% in European populations. This risk allele was positively associated with Indigenous American ancestry, showed evidence of selection in human history, and was associated with reduced IKZF1 expression. We identified a putative causal variant in a downstream enhancer that is most active in pro-B cells and interacts with the IKZF1 promoter. This variant disrupts IKZF1 autoregulation at this enhancer and results in reduced enhancer activity in B cell progenitors. Our study reveals a genetic basis for the increased ALL risk in Hispanic/Latino children.

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.

Low heterozygosity for rs3811050, a 5 prime untranslated region variant of the gene encoding interleukin-38 (IL1F10), is associated with a reduced risk of systemic lupus erythematosus

BackgroundInterleukin-38 (IL-38), an inflammatory cytokine discovered in recent years, has been implicated in the pathogenesis of systemic lupus erythematosus (SLE). IL-38 is encoded by the IL1F10 (interleukin 1 family member 10) gene. Genetic variants of this gene have been associated with susceptibility to a number of autoimmune and inflammatory diseases, while their association with SLE risk has not been explored. In this case–control study, two novel variants of the 5 prime untranslated region (5′UTR) of the IL1F10 gene, rs3811050 C/T and rs3811051 T/G, were investigated in 120 women with SLE and 120 age-matched control women. The TaqMan allelic discrimination assay was used for genotyping of rs3811050 and rs3811051.ResultsThe frequency of the rs3811050 CT genotype was significantly lower in SLE patients compared to controls (30.8 vs. 50.0%; odds ratio = 0.49; 95% confidence interval = 0.28–0.86; corrected probability = 0.045). The rs3811051 genotype frequencies did not show significant differences between patients and controls. Rs3811050 and rs3811051 showed weak linkage disequilibrium (LD) as indicated by the estimated LD coefficient and correlation coefficient values (0.32 and 0.05, respectively), and two-locus haplotype analysis revealed no significant differences between patients and controls. The frequencies of the rs3811050 T allele (38.8 vs. 20.6%; probability = 0.029) and the rs3811051 G allele (56.3 vs. 38.2%; probability = 0.038) were significantly higher in patients with mild/moderate disease activity than in patients with high disease activity, but significance was not maintained after applying Bonferroni correction (corrected probability = 0.058 and 0.076, respectively). Serum IL-38 concentrations (median and interquartile range) were significantly decreased in patients compared with controls (69.5 [64.1–74.8] vs. 73.5 [66.1–82.9] pg/mL; probability = 0.03), but were not influenced by SNP genotypes.ConclusionsThe heterozygous genotype of rs3811050, a 5'UTR variant, of the IL-38 encoding gene, IL1F10, is associated with a reduced risk of SLE among women. Furthermore, the rs3811050 T and rs3811051 G alleles may influence disease activity. In addition, serum IL-38 concentrations were down-regulated in SLE patients but were not affected by the rs3811050 and rs3811051 genotypes.

Single nucleotide variants in nuclear pore complex disassembly pathway associated with poor survival in osteosarcoma.

The bone tumor, osteosarcoma, remains challenging to treat in children and young adults, especially when patients present with metastatic disease. Developing new therapies based on genomic data from sequencing projects has proven difficult given the lack of recurrent genetic lesions across tumors. MYC overexpression has been associated with poor outcomes in osteosarcoma. However, other genomic markers of disease severity are lacking. We utilized whole genome sequencing of 106 tumors and matched normal controls in order to define genomic characteristics that correlate with overall survival. Single nucleotide variants were overlaid onto annotated molecular pathways in order to define aberrant pathway signatures specific to aggressive osteosarcoma. Additionally, we calculated differential gene expression in a subsample of 71 tumors. Differentially expressed genes were then queried for known MYC-responsive genes. Molecular pathways specific to nuclear pore complex disassembly (NPCD) show significant correlation with poor overall survival in osteosarcoma when mutations were present. Genes involved in immune response and immune regulation are enriched in the differential expression analysis of samples with and without NPCD pathway aberrations. Furthermore, neither MYC nor MYC-responsive genes show differential expression between NPCD-aberrant and non-aberrant groups. The NPCD pathway mutations are dominated by regulatory region variants rather than protein-altering mutations, suggesting that dysregulation of genetic regulatory networks may be the underlying mechanism for their relation to osteosarcoma phenotype. Overall survival is significantly worse in patients whose tumors show aberrations in the NPCD pathway. Moreover, this difference in survival is not driven by MYC-overexpression, suggesting a novel mechanism for some aggressive osteosarcomas. These findings add light to the evolving understanding of the drivers of osteosarcoma and may aid in the search for new treatments based on patient-specific genetic data.

IFI16 Is Indispensable for Promoting HIF-1α-Mediated APOL1 Expression in Human Podocytes under Hypoxic Conditions.

Genetic variants in the protein-coding regions of APOL1 are associated with an increased risk and progression of chronic kidney disease (CKD) in African Americans. Hypoxia exacerbates CKD progression by stabilizing HIF-1α, which induces APOL1 transcription in kidney podocytes. However, the contribution of additional mediators to regulating APOL1 expression under hypoxia in podocytes is unknown. Here, we report that a transient accumulation of HIF-1α in hypoxia is sufficient to upregulate APOL1 expression in podocytes through a cGAS/STING/IRF3-independent pathway. Notably, IFI16 ablation impedes hypoxia-driven APOL1 expression despite the nuclear accumulation of HIF-1α. Co-immunoprecipitation assays indicate no direct interaction between IFI16 and HIF-1α. Our studies identify hypoxia response elements (HREs) in the APOL1 gene enhancer/promoter region, showing increased HIF-1α binding to HREs located in the APOL1 gene enhancer. Luciferase reporter assays confirm the role of these HREs in transcriptional activation. Chromatin immunoprecipitation (ChIP)-qPCR assays demonstrate that IFI16 is not recruited to HREs, and IFI16 deletion reduces HIF-1α binding to APOL1 HREs. RT-qPCR analysis indicates that IFI16 selectively affects APOL1 expression, with a negligible impact on other hypoxia-responsive genes in podocytes. These findings highlight the unique contribution of IFI16 to hypoxia-driven APOL1 gene expression and suggest alternative IFI16-dependent mechanisms regulating APOL1 gene expression under hypoxic conditions.

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.

The axillary artery high bifurcation: coexisting variants and clinical significance.

Several branching patterns of the axillary artery (AA) have been described. Unusually, the brachial artery (BA) follows a course in front of the median nerve (MN), the so-called superficial brachial artery (SBA). The SBA may result in MN entrapment. The current cadaveric report highlights a high AA bifurcation, its continuation as SBA and the coexistence of muscular, neural, and vascular asymmetric aberrations. At the right side, the coracobrachialis muscle (CB) had a single head, and the ipsilateral musculocutaneous nerve (MCN) followed a medial course. The AA was highly divided into superficial and deep stems (SAS and DAS), at the 2nd rib lower border. Between two stems, the brachial plexus (BP) lateral and medial cords were identified. The MN originated from the BP lateral cord. The SAS, continued as SBA with a tortuous course. The DAS coursed posterior to the BP medial and lateral cords and gave off the subscapular artery. A bilateral 3rd head of the biceps brachii was identified. The MN atypically originated from the BP lateral cord. At the left side, the two-headed CB was typically penetrated by the MCN. A common trunk of the circumflex humeral arteries was identified in coexistence with an interconnection of the BP lateral cord with the MN medial root. The rare coexistence of muscular, neural, and arterial variants in axillary and brachial region is emphasized, taking into consideration the AA high division and related branching pattern. Documentation of such rare vascular variants is important in aneurysm and trauma surgery, and angiography, where all therapeutic manipulations must be accurately performed due to the possibility of complications.

Extreme restructuring of cis-regulatory regions controlling a deeply conserved plant stem cell regulator.

A striking paradox is that genes with conserved protein sequence, function and expression pattern over deep time often exhibit extremely divergent cis-regulatory sequences. It remains unclear how such drastic cis-regulatory evolution across species allows preservation of gene function, and to what extent these differences influence how cis-regulatory variation arising within species impacts phenotypic change. Here, we investigated these questions using a plant stem cell regulator conserved in expression pattern and function over ~125 million years. Using in-vivo genome editing in two distantly related models, Arabidopsis thaliana (Arabidopsis) and Solanum lycopersicum (tomato), we generated over 70 deletion alleles in the upstream and downstream regions of the stem cell repressor gene CLAVATA3 (CLV3) and compared their individual and combined effects on a shared phenotype, the number of carpels that make fruits. We found that sequences upstream of tomato CLV3 are highly sensitive to even small perturbations compared to its downstream region. In contrast, Arabidopsis CLV3 function is tolerant to severe disruptions both upstream and downstream of the coding sequence. Combining upstream and downstream deletions also revealed a different regulatory outcome. Whereas phenotypic enhancement from adding downstream mutations was predominantly weak and additive in tomato, mutating both regions of Arabidopsis CLV3 caused substantial and synergistic effects, demonstrating distinct distribution and redundancy of functional cis-regulatory sequences. Our results demonstrate remarkable malleability in cis-regulatory structural organization of a deeply conserved plant stem cell regulator and suggest that major reconfiguration of cis-regulatory sequence space is a common yet cryptic evolutionary force altering genotype-to-phenotype relationships from regulatory variation in conserved genes. Finally, our findings underscore the need for lineage-specific dissection of the spatial architecture of cis-regulation to effectively engineer trait variation from conserved productivity genes in crops.