The Revised Reference Genome of the Leopard Gecko (\u3cem\u3eEublepharis macularius\u3c/em\u3e) Provides Insight into the Considerations of Genome Phasing and Assembly

Genomic resources across squamate reptiles (lizards and snakes) have lagged behind other vertebrate systems and high-quality reference genomes remain scarce. Of the 23 chromosome-scale reference genomes across the order, only 12 of the ~60 squamate families are represented. Within geckos (infraorder Gekkota), a species-rich clade of lizards, chromosome-level genomes are exceptionally sparse representing only two of the seven extant families. Using the latest advances in genome sequencing and assembly methods, we generated one of the highest-quality squamate genomes to date for the leopard gecko, Eublepharis macularius (Eublepharidae). We compared this assembly to the previous, short-read only, E. macularius reference genome published in 2016 and examined potential factors within the assembly influencing contiguity of genome assemblies using PacBio HiFi data. Briefly, the read N50 of the PacBio HiFi reads generated for this study was equal to the contig N50 of the previous E. macularius reference genome at 20.4 kilobases. The HiFi reads were assembled into a total of 132 contigs, which was further scaffolded using HiC data into 75 total sequences representing all 19 chromosomes. We identified 9 of the 19 chromosomal scaffolds were assembled as a near-single contig, whereas the other 10 chromosomes were each scaffolded together from multiple contigs. We qualitatively identified that the percent repeat content within a chromosome broadly affects its assembly contiguity prior to scaffolding. This genome assembly signifies a new age for squamate genomics where high-quality reference genomes rivaling some of the best vertebrate genome assemblies can be generated for a fraction of previous cost estimates. This new E. macularius reference assembly is available on NCBI at JAOPLA010000000

The Impact of Donor and Recipient Genetic Variation on Outcomes After Solid Organ Transplantation:a Scoping Review and Future Perspectives

At the outset of solid organ transplantation, genetic variation between donors and recipients was recognized as a major player in mechanisms such as allograft tolerance and rejection. Genome-wide association studies have been very successful in identifying novel variant-trait associations, but have been difficult to perform in the field of solid organ transplantation due to complex covariates, era effects, and poor statistical power for detecting donor-recipient interactions. To overcome a lack of statistical power, consortia such as the International Genetics and Translational Research in Transplantation Network have been established. Studies have focused on the consequences of genetic dissimilarities between donors and recipients and have reported associations between polymorphisms in candidate genes or their regulatory regions with transplantation outcomes. However, knowledge on the exact influence of genetic variation is limited due to a lack of comprehensive characterization and harmonization of recipients' or donors' phenotypes and validation using an experimental approach. Causal research in genetics has evolved from agnostic discovery in genome-wide association studies to functional annotation and clarification of underlying molecular mechanisms in translational studies. In this overview, we summarize how the recent advances and progresses in the field of genetics and genomics have improved the understanding of outcomes after solid organ transplantation

Gene-centric meta-analyses of 108,912 individuals confirm known body mass index loci and reveal three novel signals

Recent genetic association studies have made progress in uncovering components of the genetic architecture of body mass index (BMI). We used the ITMAT-Broad-CARe (IBC) array comprising up to 49,320 single nucleotide polymorphisms (SNPs) across ~2,100 metabolic and cardiovascular-related loci to genotype up to 108,912 individuals of European ancestry (EA), African Americans, Hispanics, and East Asians, from 46 studies, to provide additional insight into SNPs underpinning BMI. We used a five-phase study design: Phase I focused on meta-analysis of EA studies providing individual level genotype data; Phase II performed a replication of cohorts providing summary level EA data; Phase III meta-analyzed results from the first two phases; associated SNPs from Phase III were used for replication in Phase IV; finally in Phase V, a multi-ethnic meta-analysis of all samples from four ethnicities was performed. At an array-wide significance (P<2.40E-06), we identify novel BMI associations in loci TOMM40-APOE-APOC1 (rs2075650, P=2.95E-10), SREBF2 (a sterol regulatory element binding transcription factor gene, rs5996074, P=9.43E-07) and NTRK2 (a BDNF receptor, rs1211166, P=1.04E-06) in the Phase IV meta-analysis. Of ten loci with previous evidence for BMI association represented on IBC array, eight were replicated, with the remaining two showing nominal significance. Conditional analyses revealed two independent BMI associated signals in BDNF and MC4R regions. Of the 11 array-wide significant SNPs, three are associated with gene expression levels in both primary B-cells and monocytes; with rs4788099 in SH2B1 notably being associated with the expression of multiple genes in cis. These multi-ethnic meta-analyses expand our knowledge of BMI genetics

IBC CARe Microarray Allelic Population Prevalences in an American Indian Population

Background: The prevalence of variant alleles among single nucleotide polymorphisms (SNPs) is not well known for many minority populations. These population allele frequencies (PAFs) are necessary to guide genetic epidemiology studies and to understand the population specific contribution of these variants to disease risk. Large differences in PAF among certain functional groups of genes could also indicate possible selection pressure or founder effects of interest. The 50K SNP, custom genotyping microarray (CARe) was developed, focusing on about 2,000 candidate genes and pathways with demonstrated pathophysiologic influence on cardiovascular disease (CVD). Methods: The CARe microarray was used to genotype 216 unaffected controls in a study of pre-eclampsia among a Northern Plains, American Indian tribe. The allelic prevalences of 34,240 SNPs suitable for analysis, were determined and compared with corresponding HapMap prevalences for the Caucasian population. Further analysis was conducted to compare the frequency of statistically different prevalences among functionally related SNPs, as determined by the DAVID Bioinformatics Resource. Results: Of the SNPs with PAFs in both datasets, 9.8%,37.2% and 47.1% showed allele frequencies among the American Indian population greater than, less than and either greater or less than (respectively) the HapMap Caucasian population. The 2,547 genes were divided into 53 functional groups using the highest stringency criteria. While none of these groups reached the Bonferroni corrected p value of 0.00094, there were 7 of these 53 groups with significantly more or less differing PAFs, each with a probability of less than 0.05 and an overall probability of 0.0046. Conclusion: In comparison to the HapMap Caucasian population, there are substantial differences in the prevalence among an American Indian community of SNPs related to CVD. Certain functional groups of genes and related SNPs show possible evidence of selection pressure or founder effects

Harnessing publicly available genetic data to prioritize lipid modifying therapeutic targets for prevention of coronary heart disease based on dysglycemic risk

Therapeutic interventions that lower LDL-cholesterol effectively reduce the risk of coronary artery disease (CAD). However, statins, the most widely prescribed LDL-cholesterol lowering drugs, increase diabetes risk. We used genome-wide association study (GWAS) data in the public domain to investigate the relationship of LDL-C and diabetes and identify loci encoding potential drug targets for LDL-cholesterol modification without causing dysglycemia. We obtained summary-level GWAS data for LDL-C from GLGC, glycemic traits from MAGIC, diabetes from DIAGRAM and CAD from CARDIoGRAMplusC4D consortia. Mendelian randomization analyses identified a one standard deviation (SD) increase in LDL-C caused an increased risk of CAD (odds ratio [OR] 1.63 (95 % confidence interval [CI] 1.55, 1.71), which was not influenced by removing SNPs associated with diabetes. LDL-C/CAD-associated SNPs showed consistent effect directions (binomial P = 6.85 × 10−5). Conversely, a 1-SD increase in LDL-C was causally protective of diabetes (OR 0.86; 95 % CI 0.81, 0.91), however LDL-cholesterol/diabetes-associated SNPs did not show consistent effect directions (binomial P = 0.15). HMGCR, our positive control, associated with LDL-C, CAD and a glycemic composite (derived from GWAS meta-analysis of four glycemic traits and diabetes). In contrast, PCSK9, APOB, LPA, CETP, PLG, NPC1L1 and ALDH2 were identified as “druggable” loci that alter LDL-C and risk of CAD without displaying associations with dysglycemia. In conclusion, LDL-C increases the risk of CAD and the relationship is independent of any association of LDL-C with diabetes. Loci that encode targets of emerging LDL-C lowering drugs do not associate with dysglycemia, and this provides provisional evidence that new LDL-C lowering drugs (such as PCSK9 inhibitors) may not influence risk of diabetes

Unsupervised mRNA-seq classification of heart transplant endomyocardial biopsies.

BACKGROUND: Endomyocardial biopsy (EMB) is currently considered the gold standard for diagnosing cardiac allograft rejection. However, significant limitations related to histological interpretation variability are well-recognized. We sought to develop a methodology to evaluate EMB solely based on gene expression, without relying on histology interpretation. METHODS: Sixty-four EMBs were obtained from 47 post-heart transplant recipients, who were evaluated for allograft rejection. EMBs were subjected to mRNA sequencing, in which an unsupervised classification algorithm was used to identify the molecular signatures that best classified the EMBs. Cytokine and natriuretic peptide peripheral blood profiling was also performed. Subsequently, we performed gene network analysis to identify the gene modules and gene ontology to understand their biological relevance. We correlated our findings with the unsupervised and histological classifications. RESULTS: Our algorithm classifies EMBs into three categories based solely on clusters of gene expression: unsupervised classes 1, 2, and 3. Unsupervised and histological classifications were closely related, with stronger gene module-phenotype correlations for the unsupervised classes. Gene ontology enrichment analysis revealed processes impacting on the regulation of cardiac and mitochondrial function, immune response, and tissue injury response. Significant levels of cytokines and natriuretic peptides were detected following the unsupervised classification. CONCLUSION: We have developed an unsupervised algorithm that classifies EMBs into three distinct categories, without relying on histology interpretation. These categories were highly correlated with mitochondrial, immune, and tissue injury response. Significant cytokine and natriuretic peptide levels were detected within the unsupervised classification. If further validated, the unsupervised classification could offer a more objective EMB evaluation