rqt: an R package for gene-level meta-analysis

Matrisome and Immune Pathways Contribute to Extreme Vascular Outcomes in Williams-Beuren Syndrome

BACKGROUND

Supravalvar aortic stenosis (SVAS) is a characteristic feature of Williams-Beuren syndrome (WBS). Its severity varies: ~20% of people with Williams-Beuren syndrome have SVAS requiring surgical intervention, whereas ~35% have no appreciable SVAS. The remaining individuals have SVAS of intermediate severity. Little is known about genetic modifiers that contribute to this variability.

METHODS AND RESULTS

We performed genome sequencing on 473 individuals with Williams-Beuren syndrome and developed strategies for modifier discovery in this rare disease population. Approaches include extreme phenotyping and nonsynonymous variant prioritization, followed by gene set enrichment and pathway-level association tests. We next used GTEx v8 and proteomic data sets to verify expression of candidate modifiers in relevant tissues. Finally, we evaluated overlap between the genes/pathways identified here and those ascertained through larger aortic disease/trait genome-wide association studies. We show that SVAS severity in Williams-Beuren syndrome is associated with increased frequency of common and rarer variants in matrisome and immune pathways. Two implicated matrisome genes (ACAN and LTBP4) were uniquely expressed in the aorta. Many genes in the identified pathways were previously reported in genome-wide association studies for aneurysm, bicuspid aortic valve, or aortic size.

CONCLUSIONS

Smaller sample sizes in rare disease studies necessitate new approaches to detect modifiers. Our strategies identified variation in matrisome and immune pathways that are associated with SVAS severity. These findings suggest that, like other aortopathies, SVAS may be influenced by the balance of synthesis and degradation of matrisome proteins. Leveraging multiomic data and results from larger aorta-focused genome-wide association studies may accelerate modifier discovery for rare aortopathies like SVAS.

Antimicrobial Peptide Arsenal Predicted from the Venom Gland Transcriptome of the Tropical Trap-Jaw Ant Odontomachus chelifer

With about 13,000 known species, ants are the most abundant venomous insects. Their venom consists of polypeptides, enzymes, alkaloids, biogenic amines, formic acid, and hydrocarbons. In this study, we investigated, using in silico techniques, the peptides composing a putative antimicrobial arsenal from the venom gland of the neotropical trap-jaw ant Odontomachus chelifer. Focusing on transcripts from the body and venom gland of this insect, it was possible to determine the gland secretome, which contained about 1022 peptides with putative signal peptides. The majority of these peptides (75.5%) were unknown, not matching any reference database, motivating us to extract functional insights via machine learning-based techniques. With several complementary methodologies, we investigated the existence of antimicrobial peptides (AMPs) in the venom gland of O. chelifer, finding 112 non-redundant candidates. Candidate AMPs were predicted to be more globular and hemolytic than the remaining peptides in the secretome. There is evidence of transcription for 97% of AMP candidates across the same ant genus, with one of them also verified as translated, thus supporting our findings. Most of these potential antimicrobial sequences (94.8%) matched transcripts from the ant's body, indicating their role not solely as venom toxins.

Transcriptomic and biochemical analysis from the venom gland of the neotropical ant Odontomachus chelifer

The genus Odontomachus is widely distributed in neotropical areas throughout Central and South America. It is a stinging ant that subdues its prey (insects) by injecting them a cocktail of toxic molecules (venom). Ant venoms are generally composed of formic acid, alkaloids, hydrocarbons, amines, peptides, and proteins. Odontomachus chelifer is an ant that inhabits neotropical regions from Mexico to Argentina. Unlike the venom of other animals such as scorpions, spiders and snakes, this ant venom has seldom been analyzed comprehensively, and their compositions are not yet completely known. In the present study, we performed a partial investigation of enzymatic and functional activities of O. chelifer ant venom, and we provide a global insight on the transcripts expressed in the venom gland to better understand their properties. The crude venom showed phospholipase A₂ and antiparasitic activities. RNA sequencing (Illumina platform) of the venom gland of O. chelifer generated 61, 422, 898 reads and de novo assembly Trinity generated 50,220 contigs. BUSCO analysis against Arthropoda_db10 showed that 92.89% of the BUSCO groups have complete gene representation (single-copy or duplicated), while 4.05% are only partially recovered, and 3.06% are missing. The 30 most expressed genes in O. chelifer venom gland transcriptome included important transcripts involved in venom function such as U-poneritoxin (01)-Om1a-like (pilosulin), chitinase 2, venom allergen 3, chymotrypsin 1 and 2 and glutathione S-transferase. Analysis of the molecular function revealed that the largest number of transcripts were related to catalytic activity, including phospholipases. These data emphasize the potential of O. chelifer venom for prospection of molecules with biotechnological application.

From SNP to pathway-based GWAS meta-analysis: do current meta-analysis approaches resolve power and replication in genetic association studies?

Genome-wide association studies (GWAS) have benefited greatly from enhanced high-throughput technology in recent decades. GWAS meta-analysis has become increasingly popular to highlight the genetic architecture of complex traits, informing about the replicability and variability of effect estimations across human ancestries. A wealth of GWAS meta-analysis methodologies have been developed depending on the input data and the outcome information of interest. We present a survey of current approaches from SNP to pathway-based meta-analysis by acknowledging the range of resources and methodologies in the field, and we provide a comprehensive review of different categories of Genome-Wide Meta-analysis methods employed. These methods highlight different levels at which GWAS meta-analysis may be done, including Single Nucleotide Polymorphisms, Genes and Pathways, for which we describe their framework outline. We also discuss the strengths and pitfalls of each approach and make suggestions regarding each of them.

Dissecting Meta-Analysis in GWAS Era: Bayesian Framework for Gene/Subnetwork-Specific Meta-Analysis

Over the past decades, advanced high-throughput technologies have continuously contributed to genome-wide association studies (GWASs). GWAS meta-analysis has been increasingly adopted, has cross-ancestry replicability, and has power to illuminate the genetic architecture of complex traits, informing about the reliability of estimation effects and their variability across human ancestries. However, detecting genetic variants that have low disease risk still poses a challenge. Designing a meta-analysis approach that combines the effect of various SNPs within genes or genes within pathways from multiple independent population GWASs may be helpful in identifying associations with small effect sizes and increasing the association power. Here, we proposed ancMETA, a Bayesian graph-based framework, to perform the gene/pathway-specific meta-analysis by combining the effect size of multiple SNPs within genes, and genes within subnetwork/pathways across multiple independent population GWASs to deconvolute the interactions between genes underlying the pathogenesis of complex diseases across human populations. We assessed the proposed framework on simulated datasets, and the results show that the proposed model holds promise for increasing statistical power for meta-analysis of genetic variants underlying the pathogenesis of complex diseases. To illustrate the proposed meta-analysis framework, we leverage seven different European bipolar disorder (BD) cohorts, and we identify variants in the angiotensinogen (AGT) gene to be significantly associated with BD across all 7 studies. We detect a commonly significant BD-specific subnetwork with the ESR1 gene as the main hub of a subnetwork, associated with neurotrophin signaling (p = 4e⁻¹⁴) and myometrial relaxation and contraction (p = 3e⁻⁰⁸) pathways. ancMETA provides a new contribution to post-GWAS methodologies and holds promise for comprehensively examining interactions between genes underlying the pathogenesis of genetic diseases and also underlying ethnic differences.

Whole exome sequencing in patients with Williams-Beuren syndrome followed by disease modeling in mice points to four novel pathways that may modify stenosis risk

Supravalvular aortic stenosis (SVAS) is a narrowing of the aorta caused by elastin (ELN) haploinsufficiency. SVAS severity varies among patients with Williams-Beuren syndrome (WBS), a rare disorder that removes one copy of ELN and 25-27 other genes. Twenty percent of children with WBS require one or more invasive and often risky procedures to correct the defect while 30% have no appreciable stenosis, despite sharing the same basic genetic lesion. There is no known medical therapy. Consequently, identifying genes that modify SVAS offers the potential for novel modifier-based therapeutics. To improve statistical power in our rare-disease cohort (N = 104 exomes), we utilized extreme-phenotype cohorting, functional variant filtration and pathway-based analysis. Gene set enrichment analysis of exome-wide association data identified increased adaptive immune system variant burden among genes associated with SVAS severity. Additional enrichment, using only potentially pathogenic variants known to differ in frequency between the extreme phenotype subsets, identified significant association of SVAS severity with not only immune pathway genes, but also genes involved with the extracellular matrix, G protein-coupled receptor signaling and lipid metabolism using both SKAT-O and RQTest. Complementary studies in Eln+/-; Rag1-/- mice, which lack a functional adaptive immune system, showed improvement in cardiovascular features of ELN insufficiency. Similarly, studies in mixed background Eln+/- mice confirmed that variations in genes that increase elastic fiber deposition also had positive impact on aortic caliber. By using tools to improve statistical power in combination with orthogonal analyses in mice, we detected four main pathways that contribute to SVAS risk.

Inguinal Ring RNA Sequencing Reveals Downregulation of Muscular Genes Related to Scrotal Hernia in Pigs

Scrotal hernias (SH) are common congenital defects in commercial pigs, characterized by the presence of abdominal contents in the scrotal sac, leading to considerable production and animal welfare losses. Since the etiology of SH remains obscure, we aimed to identify the biological and genetic mechanisms involved in its occurrence through the whole transcriptome analysis of SH affected and unaffected pigs’ inguinal rings. From the 22,452 genes annotated in the pig reference genome, 13,498 were expressed in the inguinal canal tissue. Of those, 703 genes were differentially expressed (DE, FDR < 0.05) between the two groups analyzed being, respectively, 209 genes upregulated and 494 downregulated in the SH-affected group. Thirty-seven significantly overrepresented GO terms related to SH were enriched, and the most relevant biological processes were muscular system, cell differentiation, sarcome reorganization, and myofibril assembly. The calcium signaling, hypertrophic cardiomyopathy, dilated cardiomyopathy, and cardiac muscle contraction were the major pathways possibly involved in the occurrence of the scrotal hernias. The expression profile of the DE genes was associated with the reduction of smooth muscle differentiation, followed by low calcium content in the cell, which could lead to a decreased apoptosis ratio and diminished muscle contraction of the inguinal canal region. We have demonstrated that genes involved with musculature are closely linked to the physiological imbalance predisposing to scrotal hernia. According to our study, the genes MYBPC1, BOK, SLC25A4, SLC8A3, DES, TPM2, MAP1CL3C, and FGF1 were considered strong candidates for future evaluation.

Proximal femoral head transcriptome reveals novel candidate genes related to epiphysiolysis in broiler chickens

BACKGROUND

The proximal femoral head separation (FHS) or epiphysiolysis is a prevalent disorder affecting the chicken femur epiphysis, being considered a risk factor to infection which can cause bacterial chondronecrosis with osteomyelitis in broilers. To identify the genetic mechanisms involved in epiphysiolysis, differentially expressed (DE) genes in the femur of normal and FHS-affected broilers were identified using RNA-Seq technology. Femoral growth plate (GP) samples from 35-day-old commercial male broilers were collected from 4 healthy and 4 FHS-affected broilers. Sequencing was performed using an Illumina paired-end protocol. Differentially expressed genes were obtained using the edgeR package based on the False Discovery Rate (FDR < 0.05).

RESULTS

Approximately 16 million reads/sample were generated with 2 × 100 bp paired-end reads. After data quality control, approximately 12 million reads/sample were mapped to the reference chicken genome (Galgal5). A total of 12,645 genes were expressed in the femur GP. Out of those, 314 were DE between groups, being 154 upregulated and 160 downregulated in FHS-affected broilers. In the functional analyses, several biological processes (BP) were overrepresented. Among them, those related to cell adhesion, extracellular matrix (ECM), bone development, blood circulation and lipid metabolism, which are more related to chicken growth, are possibly involved with the onset of FHS. On the other hand, BP associated to apoptosis or cell death and immune response, which were also found in our study, could be related to the consequence of the FHS.

CONCLUSIONS

Genes with potential role in the epiphysiolysis were identified through the femur head transcriptome analysis, providing a better understanding of the mechanisms that regulate bone development in fast-growing chickens. In this study, we highlighted the importance of cell adhesion and extracellular matrix related genes in triggering FHS. Furthermore, we have shown new insights on the involvement of lipidemia and immune response/inflammation with FHS in broilers. Understanding the changes in the GP transcriptome might support breeding strategies to address poultry robustness and to obtain more resilient broilers.