Genetic Drivers of Epigenetic and Transcriptional Variation in Human Immune Cells.

Characterizing the multifaceted contribution of genetic and epigenetic factors to disease phenotypes is a major challenge in human genetics and medicine. We carried out high-resolution genetic, epigenetic, and transcriptomic profiling in three major human immune cell types (CD14+ monocytes, CD16+ neutrophils, and naive CD4+ T cells) from up to 197 individuals. We assess, quantitatively, the relative contribution of cis-genetic and epigenetic factors to transcription and evaluate their impact as potential sources of confounding in epigenome-wide association studies. Further, we characterize highly coordinated genetic effects on gene expression, methylation, and histone variation through quantitative trait locus (QTL) mapping and allele-specific (AS) analyses. Finally, we demonstrate colocalization of molecular trait QTLs at 345 unique immune disease loci. This expansive, high-resolution atlas of multi-omics changes yields insights into cell-type-specific correlation between diverse genomic inputs, more generalizable correlations between these inputs, and defines molecular events that may underpin complex disease risk.

Population Variation and Genetic Control of Modular Chromatin Architecture in Humans.

Chromatin state variation at gene regulatory elements is abundant across individuals, yet we understand little about the genetic basis of this variability. Here, we profiled several histone modifications, the transcription factor (TF) PU.1, RNA polymerase II, and gene expression in lymphoblastoid cell lines from 47 whole-genome sequenced individuals. We observed that distinct cis-regulatory elements exhibit coordinated chromatin variation across individuals in the form of variable chromatin modules (VCMs) at sub-Mb scale. VCMs were associated with thousands of genes and preferentially cluster within chromosomal contact domains. We mapped strong proximal and weak, yet more ubiquitous, distal-acting chromatin quantitative trait loci (cQTL) that frequently explain this variation. cQTLs were associated with molecular activity at clusters of cis-regulatory elements and mapped preferentially within TF-bound regions. We propose that local, sequence-independent chromatin variation emerges as a result of genetic perturbations in cooperative interactions between cis-regulatory elements that are located within the same genomic domain.

A brief history of human disease genetics.

A primary goal of human genetics is to identify DNA sequence variants that influence biomedical traits, particularly those related to the onset and progression of human disease. Over the past 25 years, progress in realizing this objective has been transformed by advances in technology, foundational genomic resources and analytical tools, and by access to vast amounts of genotype and phenotype data. Genetic discoveries have substantially improved our understanding of the mechanisms responsible for many rare and common diseases and driven development of novel preventative and therapeutic strategies. Medical innovation will increasingly focus on delivering care tailored to individual patterns of genetic predisposition.

2021 Curt Stern Award: Studying the biology of "junk".

This article is based on the address given by the author at the 2021 virtual meeting of the American Society of Human Genetics (ASHG). The video of the original address can be found at the ASHG website.

Regulation of HLA class I expression by non-coding gene variations.

The Human Leukocyte Antigen (HLA) is a critical genetic system for different outcomes after solid organ and hematopoietic cell transplantation. Its polymorphism is usually determined by molecular technologies at the DNA level. A potential role of HLA allelic expression remains under investigation in the context of the allogenic immune response between donors and recipients. In this study, we quantified the allelic expression of all three HLA class I loci (HLA-A, B and C) by RNA sequencing and conducted an analysis of expression quantitative traits loci (eQTL) to investigate whether HLA expression regulation could be associated with non-coding gene variations. HLA-B alleles exhibited the highest expression levels followed by HLA-C and HLA-A alleles. The max fold expression variation was observed for HLA-C alleles. The expression of HLA class I loci of distinct individuals demonstrated a coordinated and paired expression of both alleles of the same locus. Expression of conserved HLA-A~B~C haplotypes differed in distinct PBMC's suggesting an individual regulated expression of both HLA class I alleles and haplotypes. Cytokines TNFα /IFNß, which induced a very similar upregulation of HLA class I RNA and cell surface expression across alleles did not modify the individually coordinated expression at the three HLA class I loci. By identifying cis eQTLs for the HLA class I genes, we show that the non-coding eQTLs explain 29%, 13%, and 31% of the respective HLA-A, B, C expression variance in unstimulated cells, and 9%, 23%, and 50% of the variance in cytokine-stimulated cells. The eQTLs have significantly higher effect sizes in stimulated cells compared to unstimulated cells for HLA-B and HLA-C genes expression. Our data also suggest that the identified eQTLs are independent from the coding variation which defines HLA alleles and thus may be influential on intra-allele expression variability although they might not represent the causal eQTLs.

Genetic Landscape of the ACE2 Coronavirus Receptor.

BACKGROUND: SARS-CoV-2, the causal agent of COVID-19, enters human cells using the ACE2 (angiotensin-converting enzyme 2) protein as a receptor. ACE2 is thus key to the infection and treatment of the coronavirus. ACE2 is highly expressed in the heart and respiratory and gastrointestinal tracts, playing important regulatory roles in the cardiovascular and other biological systems. However, the genetic basis of the ACE2 protein levels is not well understood. METHODS: We have conducted the largest genome-wide association meta-analysis of plasma ACE2 levels in >28 000 individuals of the SCALLOP Consortium (Systematic and Combined Analysis of Olink Proteins). We summarize the cross-sectional epidemiological correlates of circulating ACE2. Using the summary statistics-based high-definition likelihood method, we estimate relevant genetic correlations with cardiometabolic phenotypes, COVID-19, and other human complex traits and diseases. We perform causal inference of soluble ACE2 on vascular disease outcomes and COVID-19 severity using mendelian randomization. We also perform in silico functional analysis by integrating with other types of omics data. RESULTS: We identified 10 loci, including 8 novel, capturing 30% of the heritability of the protein. We detected that plasma ACE2 was genetically correlated with vascular diseases, severe COVID-19, and a wide range of human complex diseases and medications. An X-chromosome cis-protein quantitative trait loci-based mendelian randomization analysis suggested a causal effect of elevated ACE2 levels on COVID-19 severity (odds ratio, 1.63 [95% CI, 1.10-2.42]; P=0.01), hospitalization (odds ratio, 1.52 [95% CI, 1.05-2.21]; P=0.03), and infection (odds ratio, 1.60 [95% CI, 1.08-2.37]; P=0.02). Tissue- and cell type-specific transcriptomic and epigenomic analysis revealed that the ACE2 regulatory variants were enriched for DNA methylation sites in blood immune cells. CONCLUSIONS: Human plasma ACE2 shares a genetic basis with cardiovascular disease, COVID-19, and other related diseases. The genetic architecture of the ACE2 protein is mapped, providing a useful resource for further biological and clinical studies on this coronavirus receptor.

Identifying novel regulatory effects for clinically relevant genes through the study of the Greek population.

BACKGROUND: Expression quantitative trait loci (eQTL) studies provide insights into regulatory mechanisms underlying disease risk. Expanding studies of gene regulation to underexplored populations and to medically relevant tissues offers potential to reveal yet unknown regulatory variants and to better understand disease mechanisms. Here, we performed eQTL mapping in subcutaneous (S) and visceral (V) adipose tissue from 106 Greek individuals (Greek Metabolic study, GM) and compared our findings to those from the Genotype-Tissue Expression (GTEx) resource. RESULTS: We identified 1,930 and 1,515 eGenes in S and V respectively, over 13% of which are not observed in GTEx adipose tissue, and that do not arise due to different ancestry. We report additional context-specific regulatory effects in genes of clinical interest (e.g. oncogene ST7) and in genes regulating responses to environmental stimuli (e.g. MIR21, SNX33). We suggest that a fraction of the reported differences across populations is due to environmental effects on gene expression, driving context-specific eQTLs, and suggest that environmental effects can determine the penetrance of disease variants thus shaping disease risk. We report that over half of GM eQTLs colocalize with GWAS SNPs and of these colocalizations 41% are not detected in GTEx. We also highlight the clinical relevance of S adipose tissue by revealing that inflammatory processes are upregulated in individuals with obesity, not only in V, but also in S tissue. CONCLUSIONS: By focusing on an understudied population, our results provide further candidate genes for investigation regarding their role in adipose tissue biology and their contribution to disease risk and pathogenesis.

Proteomic analysis of 92 circulating proteins and their effects in cardiometabolic diseases.

BACKGROUND: Human plasma contains a wide variety of circulating proteins. These proteins can be important clinical biomarkers in disease and also possible drug targets. Large scale genomics studies of circulating proteins can identify genetic variants that lead to relative protein abundance. METHODS: We conducted a meta-analysis on genome-wide association studies of autosomal chromosomes in 22,997 individuals of primarily European ancestry across 12 cohorts to identify protein quantitative trait loci (pQTL) for 92 cardiometabolic associated plasma proteins. RESULTS: We identified 503 (337 cis and 166 trans) conditionally independent pQTLs, including several novel variants not reported in the literature. We conducted a sex-stratified analysis and found that 118 (23.5%) of pQTLs demonstrated heterogeneity between sexes. The direction of effect was preserved but there were differences in effect size and significance. Additionally, we annotate trans-pQTLs with nearest genes and report plausible biological relationships. Using Mendelian randomization, we identified causal associations for 18 proteins across 19 phenotypes, of which 10 have additional genetic colocalization evidence. We highlight proteins associated with a constellation of cardiometabolic traits including angiopoietin-related protein 7 (ANGPTL7) and Semaphorin 3F (SEMA3F). CONCLUSION: Through large-scale analysis of protein quantitative trait loci, we provide a comprehensive overview of common variants associated with plasma proteins. We highlight possible biological relationships which may serve as a basis for further investigation into possible causal roles in cardiometabolic diseases.

Leveraging interindividual variability of regulatory activity for refining genetic regulation of gene expression in schizophrenia.

Schizophrenia is a polygenic psychiatric disorder with limited understanding about the mechanistic changes in gene expression regulation. To elucidate on this, we integrate interindividual variability of regulatory activity (ChIP-sequencing for H3K27ac histone mark) with gene expression and genotype data captured from the prefrontal cortex of 272 cases and controls. By measuring interindividual correlation among proximal chromatin peaks, we show that regulatory element activity is structured into 10,936 and 10,376 cis-regulatory domains in cases and controls, respectively. The schizophrenia-specific cis-regulatory domains are enriched for fetal-specific (p = 0.0014, OR = 1.52) and depleted of adult-specific regulatory activity (p = 3.04 × 10-50, OR = 0.57) and are enriched for SCZ heritability (p = 0.001). By studying the interplay among genetic variants, gene expression, and cis-regulatory domains, we ascertain that changes in coordinated regulatory activity tag alterations in gene expression levels (p = 3.43 × 10-5, OR = 1.65), unveil case-specific QTL effects, and identify regulatory machinery changes for genes affecting synaptic function and dendritic spine morphology in schizophrenia. Altogether, we show that accounting for coordinated regulatory activity provides a novel mechanistic approach to reduce the search space for unveiling genetically perturbed regulation of gene expression in schizophrenia.

Contribution of schizophrenia polygenic burden to longitudinal phenotypic variance in 22q11.2 deletion syndrome.

While the recurrent 22q11.2 deletion is one of the strongest genetic risk factors for schizophrenia (SCZ), variability of its associated neuropsychiatric endophenotypes reflects its incomplete penetrance for psychosis development. To assess whether this phenotypic variability is linked to common variants associated with SCZ, we studied the association between SCZ polygenic risk score (PRS) and longitudinally acquired phenotypic information of the Swiss 22q11.2DS cohort (n = 97, 50% females, mean age 17.7 yr, mean visit interval 3.8 yr). The SCZ PRS with the best predictive performance was ascertained in the Estonian Biobank (n = 201,146) with LDpred. The infinitesimal SCZ PRS model showed the strongest capacity in discriminating SCZ cases from controls with one SD difference in SCZ PRS corresponding to an odds ratio (OR) of 1.73 (95% CI 1.57-1.90, P = 1.47 × 10-29). In 22q11.2 patients, random-effects ordinal regression modelling using longitudinal data showed SCZ PRS to have the strongest effect on social anhedonia (OR = 2.09, P = 0.0002), and occupational functioning (OR = 1.82, P = 0.0003) within the negative symptoms course, and dysphoric mood (OR = 2.00, P = 0.002) and stress intolerance (OR = 1.76, P = 0.0002) within the general symptoms course. Genetic liability for SCZ was additionally associated with full scale cognitive decline (ß = -0.25, P = 0.02) and with longitudinal volumetric reduction of the right and left hippocampi (ß = -0.28, P = 0.005; ß = -0.23, P = 0.02, respectively). Our results indicate that the polygenic contribution to SCZ acts upon the threshold-lowering first hit (i.e., the deletion). It modifies the endophenotypes of 22q11.2DS and augments the derailment of developmental trajectories of negative and general symptoms, cognition, and hippocampal volume.

Efficacy and safety of fremanezumab for migraine prophylaxis in patients with at least three previous preventive failures: Prospective, multicenter, real-world data from a Greek registry.

OBJECTIVE: To prospectively assess the efficacy and safety of fremanezumab for migraine prophylaxis in patients with failure of at least three previous preventive treatments. Changes in disability as quality-of-life outcomes after fremanezumab treatment were also examined. METHODS: Two hundred and four patients with either high-frequency EM (HFEM) or chronic migraine (CM), who attained at least three consecutive monthly sessions with fremanezumab 225 mg and otherwise met the inclusion criteria, were included in the study. The crude response (at least 50% reduction in monthly headache days [MHD]) rates to fremanezumab were assessed. Scores in the following efficacy outcomes were then compared from baseline to the last efficacy evaluation follow-up: (i) MHD, (ii) monthly days with moderate/severe peak headache intensity, and (iii) monthly days with intake of abortive medication. The disability was evaluated with the Migraine Disability Assessment; the quality of life (QOL) status was assessed with the Headache Impact-6 Test, and the EQ-5D questionnaire. RESULTS: In the majority of HFEM cases (n = 81/97; 83.5%) and CM patients (n = 67/107; 62.6%), fremanezumab proved effective in reducing the MHDs by at least 50% and was associated with clinically meaningful improvement in all other efficacy variables. The migraine-related disability experienced by our patients decreased and their QOL increased. We recorded just 36 cases reporting mild adverse events, including pain, rash or pruritus (n = 26), flu-like symptoms (n = 8), and hair loss (n = 2). CONCLUSION: With our prospective results, we provide further real-world data to support the favorable benefit/risk profile of fremanezumab in the prophylaxis of both HFEM and CM.

A Survey of Greek Primary Care Physicians on Their Likeability in Treating Migraines and Other Common Diseases.

Background and objectives: Migraine is considered the most clinically important primary headache due to its high prevalence and significant burden. Although globally categorized as one of the leading causes of disability, it is still largely underdiagnosed and undertreated. Worldwide, migraine care is in most cases provided by primary care physicians. The aim of our study was to assess the attitudes of Greek primary care physicians toward treating migraine compared to other common neurological and general medical disorders. Methods: We surveyed 182 primary care physicians with the use of a 5-point questionnaire regarding their preference in treating ten common medical conditions, including migraine, hypertension, hyperlipidemia, upper respiratory tract infections, diabetes mellitus, lower back pain, dizziness, transient ischemic attack, diabetic peripheral neuropathy, and fibromyalgia. Results: Overall, with regards to preference to treat, migraine scored very low (3.6 ± 1.0), next to diabetic peripheral neuropathy (3.6 ± 1.0), and third from the bottom to fibromyalgia (3.25 ± 1.06). In contrast, physicians reported a much higher preference to treat hypertension (4.66 ± 0.60) and hyperlipidemia (4.6 ± 1.0). Conclusions: Our results indicate that Greek primary care physicians dislike treating migraines but also other neurological diseases. Topics for further investigation include the reasons for this dislike, any associations with poor patient satisfaction, treatment results, or both.

Cross-species transcriptome analysis for early detection and specific therapeutic targeting of human lupus nephritis.

OBJECTIVES: Patients with lupus nephritis (LN) are in urgent need for early diagnosis and therapeutic interventions targeting aberrant molecular pathways enriched in affected kidneys. METHODS: We used mRNA-sequencing in effector (spleen) and target (kidneys, brain) tissues from lupus and control mice at sequential time points, and in the blood from 367 individuals (261 systemic lupus erythematosus (SLE) patients and 106 healthy individuals). Comparative cross-tissue and cross-species analyses were performed. The human dataset was split into training and validation sets and machine learning was applied to build LN predictive models. RESULTS: In murine SLE, we defined a kidney-specific molecular signature, as well as a molecular signature that underlies transition from preclinical to overt disease and encompasses pathways linked to metabolism, innate immune system and neutrophil degranulation. The murine kidney transcriptome partially mirrors the blood transcriptome of patients with LN with 11 key transcription factors regulating the cross-species active LN molecular signature. Integrated protein-to-protein interaction and drug prediction analyses identified the kinases TRRAP, AKT2, CDK16 and SCYL1 as putative targets of these factors and capable of reversing the LN signature. Using murine kidney-specific genes as disease predictors and machine-learning training of the human RNA-sequencing dataset, we developed and validated a peripheral blood-based algorithm that discriminates LN patients from normal individuals (based on 18 genes) and non-LN SLE patients (based on 20 genes) with excellent sensitivity and specificity (area under the curve range from 0.80 to 0.99). CONCLUSIONS: Machine-learning analysis of a large whole blood RNA-sequencing dataset of SLE patients using human orthologs of mouse kidney-specific genes can be used for early, non-invasive diagnosis and therapeutic targeting of LN. The kidney-specific gene predictors may facilitate prevention and early intervention trials.

Single-cell transcriptomics reveal temporal dynamics of critical regulators of germ cell fate during mouse sex determination.

Despite the importance of germ cell (GC) differentiation for sexual reproduction, the gene networks underlying their fate remain unclear. Here, we comprehensively characterize the gene expression dynamics during sex determination based on single-cell RNA sequencing of 14 914 XX and XY mouse GCs between embryonic days (E) 9.0 and 16.5. We found that XX and XY GCs diverge transcriptionally as early as E11.5 with upregulation of genes downstream of the bone morphogenic protein (BMP) and nodal/Activin pathways in XY and XX GCs, respectively. We also identified a sex-specific upregulation of genes associated with negative regulation of mRNA processing and an increase in intron retention consistent with a reduction in mRNA splicing in XY testicular GCs by E13.5. Using computational gene regulation network inference analysis, we identified sex-specific, sequential waves of putative key regulator genes during GC differentiation and revealed that the meiotic genes are regulated by positive and negative master modules acting in an antagonistic fashion. Finally, we found that rare adrenal GCs enter meiosis similarly to ovarian GCs but display altered expression of master genes controlling the female and male genetic programs, indicating that the somatic environment is important for GC function. Our data are available on a web platform and provide a molecular roadmap of GC sex determination at single-cell resolution, which will serve as a valuable resource for future studies of gonad development, function, and disease.

Results of a Web-based survey of 2105 Greek migraine patients in 2020: demographics, clinical characteristics, burden and the effects of the COVID-19 pandemic on the course of migraine.

The Greek Society of Migraine and Headache Patients (GSMHP), maintaining a strong commitment to research and information, conducted its second web-based online survey named "Migraine in Greece-2020", following its first one conducted in 2018. The 2020 study included 2,105 migraine patients who were called to answer 151 questions. The purposes of the current research were to record the demographic and clinical characteristics of migraine patients in Greece, including the severity and effects of migraine on respondents' quality of life, as well as to survey the effects of the coronavirus pandemic on the course of migraine. Our population, internet-based study provides data that will hopefully contribute to better comprehend the clinical phenotype and course of migraine during the COVID-19 pandemic.

Long-term adherence, safety, and efficacy of repeated onabotulinumtoxinA over five years in chronic migraine prophylaxis.

BACKGROUND: OnabotulinumtoxinA (BoNTA) demonstrated a positive benefit-risk in chronic migraine (CM) patients in PREEMPT I and II phase III trials and many subsequent real-world studies. We herein aimed at evaluating the adherence to repeated BoNTA over a period of five years, while secondary objectives included the assessment of its long-term safety/efficacy and patients' satisfaction to treatment. METHODS: We studied 56 CM patients who had successfully received consequent cycles of BoNTA over five years. Adherence was calculated as the percentage of patients actively choosing to follow with repeated BoNTA treatment, as instructed. Safety and efficacy data were collected throughout the study period. The overall patients' belief in and satisfaction by the efficacy of treatment was assessed at last follow-up, using the self-report 7-point measure patient global impression of change (PGIC). RESULTS: A total of 36 (64.3%) out of 56 patients remained adherent to BoNTA over five years. Long-term BoNTA exposure was safe and well-tolerated, without severe side-effects justifying treatment discontinuation. The mean monthly headache days and associated clinical efficacy outcomes remained consistent and quite low at last follow-up with evidence of continuous improvements in headache monthly frequency between year three and over five years of therapy. All patients who were able to maintain treatment over five years (n = 36), remained very satisfied and scored at least 5 in PGIC. CONCLUSION: Considerably high adherence, considerable satisfaction and sustained safety/efficacy were observed in patients followed up for five years, supporting a favorable benefit/risk profile for consistently delivering long-term BoNTA in CM.

Expression estimation and eQTL mapping for HLA genes with a personalized pipeline.

The HLA (Human Leukocyte Antigens) genes are well-documented targets of balancing selection, and variation at these loci is associated with many disease phenotypes. Variation in expression levels also influences disease susceptibility and resistance, but little information exists about the regulation and population-level patterns of expression. This results from the difficulty in mapping short reads originated from these highly polymorphic loci, and in accounting for the existence of several paralogues. We developed a computational pipeline to accurately estimate expression for HLA genes based on RNA-seq, improving both locus-level and allele-level estimates. First, reads are aligned to all known HLA sequences in order to infer HLA genotypes, then quantification of expression is carried out using a personalized index. We use simulations to show that expression estimates obtained in this way are not biased due to divergence from the reference genome. We applied our pipeline to the GEUVADIS dataset, and compared the quantifications to those obtained with reference transcriptome. Although the personalized pipeline recovers more reads, we found that using the reference transcriptome produces estimates similar to the personalized pipeline (r ≥ 0.87) with the exception of HLA-DQA1. We describe the impact of the HLA-personalized approach on downstream analyses for nine classical HLA loci (HLA-A, HLA-C, HLA-B, HLA-DRA, HLA-DRB1, HLA-DQA1, HLA-DQB1, HLA-DPA1, HLA-DPB1). Although the influence of the HLA-personalized approach is modest for eQTL mapping, the p-values and the causality of the eQTLs obtained are better than when the reference transcriptome is used. We investigate how the eQTLs we identified explain variation in expression among lineages of HLA alleles. Finally, we discuss possible causes underlying differences between expression estimates obtained using RNA-seq, antibody-based approaches and qPCR.

Results of a Web-Based Survey on 2105 Greek Migraine Patients-Second Part: Efficacy of Acute and Prophylactic Migraine Treatments and Corresponding Patients' Reported Satisfaction.

Background and Objectives: The Greek Society of Migraine and Headache Patients conducted, in 2020, its second online survey, titled "Migraine in Greece-2020", after publication of the first similar online survey conducted in 2018. To compare the current findings with the corresponding data obtained in 2018, we herein release the second part of results obtained from the 2020 survey on the efficacy of preventive and symptomatic anti-migraine medications and the patients' reported satisfaction with these treatments. Materials and Methods: We surveyed 2105 migraine patients from all over Greece with the use of a 151-questions specific migraine-focused questionnaire in Greek language, which was distributed through the online research software "SurveyMonkey". Results: Triptans were mostly used with efficacy for the symptomatic relief of migraine attacks. About 2 of 3 surveyed patients had received various prophylactic oral medications and the majority of them discontinued these prophylactic medications as a result of inefficacy/safety issues. BoNTA was reported to be effective only when administration was commenced by a trained neurologist/headache specialist, while our current findings are generally comparable to those obtained in our 2018 pre-COVID-19 survey and the pandemic has not imposed any significant attitudes on migraine therapies and corresponding patients' satisfaction. Conclusion: Although a market change is anticipated with the evolving widespread use of anti-CGRPs monoclonal antibodies or gepants in the symptomatic and prophylactic treatment of migraine, it is of great interest to review published results of larger longitudinal population-based studies to further ascertain the satisfaction of patients to migraine therapies.

Correction to: The role of physical activity in metabolic homeostasis before and after the onset of type 2 diabetes: an IMI DIRECT study.

Unfortunately, 'Present address' was omitted from one of the addresses provided for Mark I. McCarthy (#26).

Overview of PAX gene family: analysis of human tissue-specific variant expression and involvement in human disease.

Paired-box (PAX) genes encode a family of highly conserved transcription factors found in vertebrates and invertebrates. PAX proteins are defined by the presence of a paired domain that is evolutionarily conserved across phylogenies. Inclusion of a homeodomain and/or an octapeptide linker subdivides PAX proteins into four groups. Often termed "master regulators", PAX proteins orchestrate tissue and organ development throughout cell differentiation and lineage determination, and are essential for tissue structure and function through maintenance of cell identity. Mutations in PAX genes are associated with myriad human diseases (e.g., microphthalmia, anophthalmia, coloboma, hypothyroidism, acute lymphoblastic leukemia). Transcriptional regulation by PAX proteins is, in part, modulated by expression of alternatively spliced transcripts. Herein, we provide a genomics update on the nine human PAX family members and PAX homologs in 16 additional species. We also present a comprehensive summary of human tissue-specific PAX transcript variant expression and describe potential functional significance of PAX isoforms. While the functional roles of PAX proteins in developmental diseases and cancer are well characterized, much remains to be understood regarding the functional roles of PAX isoforms in human health. We anticipate the analysis of tissue-specific PAX transcript variant expression presented herein can serve as a starting point for such research endeavors.