Trans-regulatory Effects Research Articles

Molecular characterization of the SP3a gene, a negative regulator of viral infection in the orange-spotted grouper, Epinephelus coioides

SP3 (specificity protein 3) is a transcription factor characterized by three conserved Cys2His2 zinc finger motifs that exert a transregulatory effect by binding to GC boxes, either upregulating or downregulating multiple genes or by co-regulating gene expression in coordination with other proteins. SP3 potentially regulates a series of processes, such as the cell cycle, growth, metabolic pathways, and apoptosis, and plays an important role in antiviral effect. The function of sp3 in fish is poorly understood. In this study, the Sp3a open reading frame was cloned from the orange-spotted grouper, Epinephelus coioides. The full-length open reading frame of Sp3a was 2034 bp, encoding 677 amino acids, with a predicted molecular weight of 72.34 kDa and an isoelectric point of 5.05. Phylogenetically, Sp3a in Epinephelus coioides was the most closely related to Sp3a in the Malabar grouper, Epinephelus malabaricus. RT-qPCR revealed ubiquitous expression of Sp3a in all examined grouper tissues, with no significant differences in expression levels among tissues. A eukaryotic expression vector, pEGFP-Sp3a, was constructed and transfected into grouper spleen (GS) cells. Subcellular localization of Sp3a was observed using an inverted fluorescence microscope. When Spa3 was overexpressed in GS cells, the expression of orange-spotted grouper nerve necrosis virus (RGNNV) genes (CP and RdRp) decreased significantly, indicating that Sp3a significantly inhibited RGNNV replication. siRNA inhibition of Sp3a accelerated the intracellular replication of RGNNV, implying the antiviral effect of Sp3a. Conclusively, our findings contribute to further research on the antiviral capabilities of Sp3a in grouper and other fish. Therefore, our research has potential implications on the development of the aquaculture industry.

Variation and Interaction of Distinct Subgenomes Contribute to Growth Diversity in Intergeneric Hybrid Fish.

Intergeneric hybridization greatly reshapes regulatory interactions among allelic and non-allelic genes. However, their effects on growth diversity remain poorly understood in animals. In this study, we conducted whole-genome sequencing and RNA sequencing (RNA-seq) analyses in diverse hybrid varieties resulting from the intergeneric hybridization of goldfish (Carassius auratus red var.) and common carp (Cyprinus carpio). These hybrid individuals were characterized by distinct mitochondrial genomes and copy number variations. Through a weighted gene correlation network analysis, we identified 3693 genes as candidate growth-regulated genes. Among them, the expression of 3672 genes in subgenome R (originating from goldfish) displayed negative correlations with growth rate, whereas 20 genes in subgenome C (originating from common carp) exhibited positive correlations. Notably, we observed intriguing patterns in the expression of slc2a12 in subgenome C, showing opposite correlations with body weight that changed with water temperatures, suggesting differential interactions between feeding activity and weight gain in response to seasonal changes for hybrid animals. In 40.31% of alleles, we observed dominant trans-regulatory effects in the regulatory interaction between distinct alleles from subgenomes R and C. Integrating analyses of allelic-specific expression and DNA methylation data revealed that the influence of DNA methylation on both subgenomes shapes the relative contribution of allelic expression to the growth rate. These findings provide novel insights into the interaction of distinct subgenomes that underlie heterosis in growth traits and contribute to a better understanding of multiple allele traits in animals.

Transcriptome analysis revealed the characteristics and functions of long non-coding RNAs in the hypothalamus during sexual maturation in goats

The hypothalamus is an essential neuroendocrine area in animals that regulates sexual development. Long non-coding RNAs (lncRNAs) are hypothesized to regulate physiological processes related to animal reproduction. However, the regulatory mechanism by which lncRNAs participate in sexual maturity in goats is poorly known, particularly from birth to sexual maturation. In this study, RNAseq analysis was conducted on the hypothalamus of four developmental stages (1day (D1, n = 5), 2 months (M2, n = 5), 4 months (M4, n = 5), and 6 months (M6, n = 5)) of Jining grey goats. The results showed that a total of 237 differentially expressed lncRNAs (DELs) were identified in the hypothalamus. Among these, 221 DELs exhibited cis-regulatory effects on 693 target genes, while 24 DELs demonstrated trans-regulatory effects on 63 target genes. The target genes of these DELs are mainly involved in biological processes related to energy metabolism, signal transduction and hormone secretion, such as sphingolipid signaling pathway, adipocytokine signaling pathway, neurotrophic signaling pathway, glutamatergic synapse, P53 signaling pathway and GnRH signaling pathway. In addition, XR_001918477.1, TCONS_00077463, XR_001918760.1, and TCONS_00029048 and their potential target genes may play a crucial role in the process of goat sexual maturation. This study advances our understanding of lncRNA in hypothalamic tissue during sexual maturation in goats and will give a theoretical foundation for improving goat reproductive features.

Paralog transcriptional differentiation in the D. melanogaster-specific gene family Sdic across populations and spermatogenesis stages

How recently originated gene copies become stable genomic components remains uncertain as high sequence similarity of young duplicates precludes their functional characterization. The tandem multigene family Sdic is specific to Drosophila melanogaster and has been annotated across multiple reference-quality genome assemblies. Here we show the existence of a positive correlation between Sdic copy number and total expression, plus vast intrastrain differences in mRNA abundance among paralogs, using RNA-sequencing from testis of four strains with variable paralog composition. Single cell and nucleus RNA-sequencing data expose paralog expression differentiation in meiotic cell types within testis from third instar larva and adults. Additional RNA-sequencing across synthetic strains only differing in their Y chromosomes reveal a tissue-dependent trans-regulatory effect on Sdic: upregulation in testis and downregulation in male accessory gland. By leveraging paralog-specific expression information from tissue- and cell-specific data, our results elucidate the intraspecific functional diversification of a recently expanded tandem gene family.

Acute Myeloid Leukemia with Isolated Trisomy 4 Is Characterized By Frequent Mutations of TET2 and a Distinct DNA Methylation Signature

Background: Trisomy 4 (+4) is a rare cytogenetic abnormality found in less than 1% of acute myeloid leukemia (AML) either as sole abnormality (isolated AML+4, iAML+4) or associated with other genetic lesions. However, iAML+4 is inadequately studied due to its rarity and little is known about its molecular profile. Aim: The aim of this study was to characterize the molecular profile of iAML+4, including gene mutations, gene expression and DNA methylation changes. Methods: For the molecular characterization, we have performed targeted amplicon sequencing (n=68 target genes), gene expression profiling (SCRB-seq) and DNA methylation profiling (Illumina Methylation EPIC arrays) on 11 iAML+4 patients from German AML study groups (AML-CG 1999 and AML-CG 2008) and 105 patients with cytogenetically normal AML (CN-AML). Results: By targeted sequencing, 35 mutations were detected in 11 genes with a median of 3 mutations per patient (range 1-7, Figure 1A). Recurrently mutated genes were TET2 (n=9 patients), NPM1 (n=5) and FLT3-ITD (n=4), DNMT3A (n=2), and WT1 (n=2). Of note, in 8/9 patients with TET2 mutations, a truncating mutation has been detected. Thus, the frequency of truncating TET2 mutations was exceptionally high (8/11; 73%) in iAML+4, as compared to CN-AML with 12/105 patients mutated (11%, P=2.4e-5, Fisher exact test), which is in line with previous findings in the literature. Moreover, 3 patients harbored multiple TET2 mutations. We have further performed targeted amplicon sequencing in 8 AML patients with non-isolated +4, where trisomy 4 was present along with other cytogenetic abnormalities. In non-isolated +4 AML, TET2 mutations were less frequently found, as compared to iAML+4 (2/8, 25%, P=0.07). Gene expression profiling revealed a distinct gene expression pattern in iAML+4 as compared to CN-AML patients. Across all chromosomes, differentially expressed genes were detected; with the highest numbers of overexpressed genes located on chromosome 4, followed by chromosomes 2, 12, 17 and 18, as well as a substantial number of downregulated genes, particularly on chromosomes 11 and 19, indicating a trans-regulatory effect of +4 (Figure 1B). DNA methylation analysis revealed a specific hypermethylation signature among iAML+4 patients, as compared to CN-AML (Figure 1C). A global increase of CpG methylation levels in iAML+4 may be likely due to the high frequency of truncating TET2 mutations, as TET2 is known to be involved in hydroxymethylation of the DNA. However, the hypermethylation signature was confirmed even when comparing iAML+4 to CN-AML within TET2 mutated patient samples, indicating DNA hypermethylation beyond the effect of truncating TET2 mutations. Hierarchical clustering of iAML+4 and CN-AML samples resulted in a separate cluster with a distinct hypermethylation signature, which is strongly enriched for iAML+4 patients. Of note, besides TET2 truncating variants, no additional mutation was found to be enriched in that cluster, indicating that +4 is the major (cyto) genetic driver mediating this pattern. Conclusion: Our comprehensive multi-omics analysis revealed a homogenous molecular profile of iAML+4 with high frequency of truncating TET2 mutations indicating oncogenic cooperativity, a global pattern of differential gene expression as well as a uniform DNA hypermethylation signature, which is generally not associated with truncating TET2 mutations but a specific feature of trisomy 4 in AML. Figure 1View largeDownload PPTFigure 1View largeDownload PPT Close modal

Symmetric subgenomes and balanced homoeolog expression stabilize the establishment of allopolyploidy in cyprinid fish

BackgroundInterspecific postzygotic reproduction isolation results from large genetic divergence between the subgenomes of established hybrids. Polyploidization immediately after hybridization may reset patterns of homologous chromosome pairing and ameliorate deleterious genomic incompatibility between the subgenomes of distinct parental species in plants and animals. However, the observation that polyploidy is less common in vertebrates raises the question of which factors restrict its emergence. Here, we perform analyses of the genome, epigenome, and gene expression in the nascent allotetraploid lineage (2.95 Gb) derived from the intergeneric hybridization of female goldfish (Carassius auratus, 1.49 Gb) and male common carp (Cyprinus carpio, 1.42 Gb), to shed light on the changes leading to the stabilization of hybrids.ResultsWe firstly identify the two subgenomes derived from the parental lineages of goldfish and common carp. We find variable unequal homoeologous recombination in somatic and germ cells of the intergeneric F1 and allotetraploid (F22 and F24) populations, reflecting high plasticity between the subgenomes, and rapidly varying copy numbers between the homoeolog genes. We also find dynamic changes in transposable elements accompanied by genome merger and duplication in the allotetraploid lineage. Finally, we observe the gradual decreases in cis-regulatory effects and increases in trans-regulatory effects along with the allotetraploidization, which contribute to increases in the symmetrical homoeologous expression in different tissues and developmental stages, especially in early embryogenesis.ConclusionsOur results reveal a series of changes in transposable elements, unequal homoeologous recombination, cis- and trans-regulations (e.g. DNA methylation), and homoeologous expression, suggesting their potential roles in mediating adaptive stabilization of regulatory systems of the nascent allotetraploid lineage. The symmetrical subgenomes and homoeologous expression provide a novel way of balancing genetic incompatibilities, providing a new insight into the early stages of allopolyploidization in vertebrate evolution.

Mutational sources of trans-regulatory variation affecting gene expression in Saccharomyces cerevisiae.

Heritable variation in a gene's expression arises from mutations impacting cis- and trans-acting components of its regulatory network. Here, we investigate how trans-regulatory mutations are distributed within the genome and within a gene regulatory network by identifying and characterizing 69 mutations with trans-regulatory effects on expression of the same focal gene in Saccharomyces cerevisiae. Relative to 1766 mutations without effects on expression of this focal gene, we found that these trans-regulatory mutations were enriched in coding sequences of transcription factors previously predicted to regulate expression of the focal gene. However, over 90% of the trans-regulatory mutations identified mapped to other types of genes involved in diverse biological processes including chromatin state, metabolism, and signal transduction. These data show how genetic changes in diverse types of genes can impact a gene's expression in trans, revealing properties of trans-regulatory mutations that provide the raw material for trans-regulatory variation segregating within natural populations.

Identification of TBX15 as an adipose master trans regulator of abdominal obesity genes

BackgroundObesity predisposes individuals to multiple cardiometabolic disorders, including type 2 diabetes (T2D). As body mass index (BMI) cannot reliably differentiate fat from lean mass, the metabolically detrimental abdominal obesity has been estimated using waist-hip ratio (WHR). Waist-hip ratio adjusted for body mass index (WHRadjBMI) in turn is a well-established sex-specific marker for abdominal fat and adiposity, and a predictor of adverse metabolic outcomes, such as T2D. However, the underlying genes and regulatory mechanisms orchestrating the sex differences in obesity and body fat distribution in humans are not well understood.MethodsWe searched for genetic master regulators of WHRadjBMI by employing integrative genomics approaches on human subcutaneous adipose RNA sequencing (RNA-seq) data (n ~ 1400) and WHRadjBMI GWAS data (n ~ 700,000) from the WHRadjBMI GWAS cohorts and the UK Biobank (UKB), using co-expression network, transcriptome-wide association study (TWAS), and polygenic risk score (PRS) approaches. Finally, we functionally verified our genomic results using gene knockdown experiments in a human primary cell type that is critical for adipose tissue function.ResultsHere, we identified an adipose gene co-expression network that contains 35 obesity GWAS genes and explains a significant amount of polygenic risk for abdominal obesity and T2D in the UKB (n = 392,551) in a sex-dependent way. We showed that this network is preserved in the adipose tissue data from the Finnish Kuopio Obesity Study and Mexican Obesity Study. The network is controlled by a novel adipose master transcription factor (TF), TBX15, a WHRadjBMI GWAS gene that regulates the network in trans. Knockdown of TBX15 in human primary preadipocytes resulted in changes in expression of 130 network genes, including the key adipose TFs, PPARG and KLF15, which were significantly impacted (FDR < 0.05), thus functionally verifying the trans regulatory effect of TBX15 on the WHRadjBMI co-expression network.ConclusionsOur study discovers a novel key function for the TBX15 TF in trans regulating an adipose co-expression network of 347 adipose, mitochondrial, and metabolically important genes, including PPARG, KLF15, PPARA, ADIPOQ, and 35 obesity GWAS genes. Thus, based on our converging genomic, transcriptional, and functional evidence, we interpret the role of TBX15 to be a main transcriptional regulator in the adipose tissue and discover its importance in human abdominal obesity.

Human-chimpanzee fused cells reveal cis-regulatory divergence underlying skeletal evolution.

Gene regulatory divergence is thought to play a central role in determining human-specific traits. However, our ability to link divergent regulation to divergent phenotypes is limited. Here, we utilized human-chimpanzee hybrid induced pluripotent stem cells to study gene expression separating these species. The tetraploid hybrid cells allowed us to separate cis- from trans-regulatory effects, and to control for non-genetic confounding factors. We differentiated these cells into cranial neural crest cells (CNCCs), the primary cell type giving rise to the face. We discovered evidence of lineage-specific selection on the hedgehog signaling pathway, including a human-specific 6-fold down-regulation of EVC2 (LIMBIN), a key hedgehog gene. Inducing a similar down-regulation of EVC2 substantially reduced hedgehog signaling output. Mice and humans lacking functional EVC2 show striking phenotypic parallels to human-chimpanzee craniofacial differences, suggesting that the regulatory divergence of hedgehog signaling may have contributed to the unique craniofacial morphology of humans.

Molecular and evolutionary processes generating variation in gene expression.

Heritable variation in gene expression is common within and between species. This variation arises from mutations that alter the form or function of molecular gene regulatory networks that are then filtered by natural selection. High-throughput methods for introducing mutations and characterizing their cis- and trans-regulatory effects on gene expression (particularly, transcription) are revealing how different molecular mechanisms generate regulatory variation, and studies comparing these mutational effects with variation seen in the wild are teasing apart the role of neutral and non-neutral evolutionary processes. This integration of molecular and evolutionary biology allows us to understand how the variation in gene expression we see today came to be and to predict how it is most likely to evolve in the future.

Regulatory Potential of Long Non-Coding RNAs (lncRNAs) in Boar Spermatozoa with Good and Poor Freezability.

Long non-coding RNAs (lncRNAs) are suggested to play an important role in the sperm biological processes. We performed de novo transcriptome assembly to characterize lncRNAs in spermatozoa, and to investigate the role of the potential target genes of the differentially expressed lncRNAs (DElncRNAs) in sperm freezability. We detected approximately 4007 DElncRNAs, which were differentially expressed in spermatozoa from boars classified as having good and poor semen freezability (GSF and PSF, respectively). Most of the DElncRNAs were upregulated in boars of the PSF group and appeared to significantly affect the sperm’s response to the cryopreservation conditions. Furthermore, we predicted that the potential target genes were regulated by DElncRNAs in cis or trans. It was found that DElncRNAs of both freezability groups had potential cis- and trans-regulatory effects on different protein-coding genes, such as COX7A2L, TXNDC8 and SOX-7. Gene Ontology (GO) enrichment revealed that the DElncRNA target genes are associated with numerous biological processes, including signal transduction, response to stress, cell death (apoptosis), motility and embryo development. Significant differences in the de novo assembled transcriptome expression profiles of the DElncRNAs between the freezability groups were confirmed by quantitative real-time PCR analysis. This study reveals the potential effects of protein-coding genes of DElncRNAs on sperm functions, which could contribute to further research on their relevance in semen freezability.

Integrative genomics analysis of various omics data and networks identify risk genes and variants vulnerable to childhood-onset asthma

BackgroundChildhood-onset asthma is highly affected by genetic components. In recent years, many genome-wide association studies (GWAS) have reported a large group of genetic variants and susceptible genes associated with asthma-related phenotypes including childhood-onset asthma. However, the regulatory mechanisms of these genetic variants for childhood-onset asthma susceptibility remain largely unknown.MethodsIn the current investigation, we conducted a two-stage designed Sherlock-based integrative genomics analysis to explore the cis- and/or trans-regulatory effects of genome-wide SNPs on gene expression as well as childhood-onset asthma risk through incorporating a large-scale GWAS data (N = 314,633) and two independent expression quantitative trait loci (eQTL) datasets (N = 1890). Furthermore, we applied various bioinformatics analyses, including MAGMA gene-based analysis, pathway enrichment analysis, drug/disease-based enrichment analysis, computer-based permutation analysis, PPI network analysis, gene co-expression analysis and differential gene expression analysis, to prioritize susceptible genes associated with childhood-onset asthma.ResultsBased on comprehensive genomics analyses, we found 31 genes with multiple eSNPs to be convincing candidates for childhood-onset asthma risk; such as, PSMB9 (cis-rs4148882 and cis-rs2071534) and TAP2 (cis-rs9267798, cis-rs4148882, cis-rs241456, and trans-10,447,456). These 31 genes were functionally interacted with each other in our PPI network analysis. Our pathway enrichment analysis showed that numerous KEGG pathways including antigen processing and presentation, type I diabetes mellitus, and asthma were significantly enriched to involve in childhood-onset asthma risk. The co-expression patterns among 31 genes were remarkably altered according to asthma status, and 25 of 31 genes (25/31 = 80.65%) showed significantly or suggestively differential expression between asthma group and control group.ConclusionsWe provide strong evidence to highlight 31 candidate genes for childhood-onset asthma risk, and offer a new insight into the genetic pathogenesis of childhood-onset asthma.

Simple inheritance, complex regulation: Supergene-mediated fire ant queen polymorphism.

The fire ant Solenopsis invicta exists in two alternate social forms: monogyne nests contain a single reproductive queen and polygyne nests contain multiple reproductive queens. This colony-level social polymorphism corresponds with individual differences in queen physiology, queen dispersal patterns and worker discrimination behaviours, all evidently regulated by an inversion-based supergene that spans more than 13Mb of a "social chromosome," contains over 400 protein-coding genes and rarely undergoes recombination. The specific mechanisms by which this supergene influences expression of the many distinctive features that characterize the alternate forms remain almost wholly unknown. To advance our understanding of these mechanisms, we explore the effects of social chromosome genotype and natal colony social form on gene expression in queens sampled as they embarked on nuptial flights, using RNA-sequencing of brains and ovaries. We observe a large effect of natal social form, that is, of the social/developmental environment, on gene expression profiles, with similarly substantial effects of genotype, including: (a) supergene-associated gene upregulation, (b) allele-specific expression and (c) pronounced extra-supergene trans-regulatory effects. These findings, along with observed spatial variation in differential and allele-specific expression within the supergene region, highlight the complex gene regulatory landscape that emerged following divergence of the inversion-mediated Sb haplotype from its homologue, which presumably largely retained the ancestral gene order. The distinctive supergene-associated gene expression trajectories we document at the onset of a queen's reproductive life expand the known record of relevant molecular correlates of a complex social polymorphism and point to putative genetic factors underpinning the alternate social syndromes.

Population-scale proteome variation in human induced pluripotent stem cells.

Human disease phenotypes are driven primarily by alterations in protein expression and/or function. To date, relatively little is known about the variability of the human proteome in populations and how this relates to variability in mRNA expression and to disease loci. Here, we present the first comprehensive proteomic analysis of human induced pluripotent stem cells (iPSC), a key cell type for disease modelling, analysing 202 iPSC lines derived from 151 donors, with integrated transcriptome and genomic sequence data from the same lines. We characterised the major genetic and non-genetic determinants of proteome variation across iPSC lines and assessed key regulatory mechanisms affecting variation in protein abundance. We identified 654 protein quantitative trait loci (pQTLs) in iPSCs, including disease-linked variants in protein-coding sequences and variants with trans regulatory effects. These include pQTL linked to GWAS variants that cannot be detected at the mRNA level, highlighting the utility of dissecting pQTL at peptide level resolution.

Identification of eQTLs associated with lipid metabolism in Longissimus dorsi muscle of pigs with different genetic backgrounds

Intramuscular fat content and its fatty acid composition affect porcine meat quality and its nutritional value. The present work aimed to identify genomic variants regulating the expression in the porcine muscle (Longissimus dorsi) of 45 candidate genes for lipid metabolism and fatty acid composition in three experimental backcrosses based on the Iberian breed. Expression genome-wide association studies (eGWAS) were performed between the muscle gene expression values, measured by real-time quantitative PCR, and the genotypes of 38,426 SNPs distributed along all chromosomes. The eGWAS identified 186 eSNPs located in ten Sus scrofa regions and associated with the expression of ACSM5, ACSS2, ATF3, DGAT2, FOS and IGF2 (FDR < 0.05) genes. Two expression quantitative trait loci (eQTLs) for IGF2 and ACSM5 were classified as cis-acting eQTLs, suggesting a mutation in the same gene affecting its expression. Conversely, ten eQTLs showed trans-regulatory effects on gene expression. When the eGWAS was performed for each backcross independently, only three common trans-eQTL regions were observed, indicating different regulatory mechanisms or allelic frequencies among the breeds. In addition, hotspot regions regulating the expression of several genes were detected. Our results provide new data to better understand the functional regulatory mechanisms of lipid metabolism genes in muscle.

Variant phasing and haplotypic expression from long-read sequencing in maize

Haplotype phasing maize genetic variants is important for genome interpretation, population genetic analysis and functional analysis of allelic activity. We performed an isoform-level phasing study using two maize inbred lines and their reciprocal crosses, based on single-molecule, full-length cDNA sequencing. To phase and analyze transcripts between hybrids and parents, we developed IsoPhase. Using this tool, we validated the majority of SNPs called against matching short-read data from embryo, endosperm and root tissues, and identified allele-specific, gene-level and isoform-level differential expression between the inbred parental lines and hybrid offspring. After phasing 6907 genes in the reciprocal hybrids, we annotated the SNPs and identified large-effect genes. In addition, we identified parent-of-origin isoforms, distinct novel isoforms in maize parent and hybrid lines, and imprinted genes from different tissues. Finally, we characterized variation in cis- and trans-regulatory effects. Our study provides measures of haplotypic expression that could increase accuracy in studies of allelic expression.

CFIm25-regulated lncRNA acv3UTR promotes gastric tumorigenesis via miR-590-5p/YAP1 axis

Accumulating evidences indicate that 3ʹUTR of the coding gene can act as crucial regulators in gastric cancer (GC). However, the detailed mechanisms and responsive targets are not well established. Here, we found that acvr1b gene 3ʹUTR (acv3UTR) was elevated in GC tissue, the expression of which was significantly correlated with advanced pTNM-stage and poor outcome in clinical patients. Forced expression of acv3UTR promoted GC cells growth in vitro and in vivo. Mechanistically, our results suggested that acv3UTR functioned as an oncogenic competing endogenous RNA via sponging miR-590-5p and enhancing YAP1 level. Tumor suppressor miR-590-5p was a molecular module in acv3UTR regulatory axis, the forced expression of which led to impairing of oncogenic potential of acv3UTR. The positive correlation of acv3UTR and YAP1 expression, and the negative correlation of acv3UTR and miR-590-5p expression, were verified in GC patients. Moreover, CFIm25 was identified as a key regulator contributing to acv3UTR aberrant expression in GC binding to UGUA-264 motif. Overall, our finding defines a mechanism for understanding the potential role of acv3UTR transcription in GC tumorigenesis, and indicates a correlation between 3ʹUTR trans-regulatory effect and GC development.

Genome-wide identification of circulating-miRNA expression quantitative trait loci reveals the role of several miRNAs in the regulation of cardiometabolic phenotypes.

To identify genetic variants that have a regulatory impact on circulating microRNAs (miRNAs) and to connect genetic risk to blood traits/biomarkers through the circulating miRNAs. Leveraging miRNA-Seq data and the 1000 Genomes imputed genotypes, we carried out genome-wide association analysis for SNPs that regulate the expression of circulating miRNAs in a sample of 710 unrelated subjects of European ancestry. Wherever possible, we used data from the Framingham and the Geuvadis studies to replicate our findings. We found at least one genome-wide significant (P < 5e-8) miRNA-eQTL (mirQTL) for 143 circulating miRNAs. Overall each mirQTL explained a small portion (<1%) of variation in miRNA levels; however, we identified a few mirQTLs that explained 4% to 20% of variation in miRNA levels in plasma. Unlike trans-mirQTLs (P = 0.7), cis-mirQTLs tend to be also associated with their counterpart mature miRNAs (P < 0.0001), this suggests trans-mirQTLs exert their effect through processes that affect the stability of mature miRNAs; whereas, cis-mirQTLs mainly regulate the expression of primary-miRNAs. Next, we used the identified mirQTLs to investigate the links between circulating miRNAs with blood traits/biomarkers through Mendelian randomization analysis. We found miR-1908-5p plays an important role in regulating low-density lipoprotein (LDL), total cholesterol (TC), fasting glucose, HbA1c, and several lipid-metabolites in blood, whereas, miR-10b-5p mediates the trans-regulatory effect of the ABO locus on several blood proteins, coronary artery disease, and TC. Moreover, we demonstrated that a higher plasma level of miR-199a is causally associated with lower levels of LDL and TC. Finally, we found miR-143-3p and miR-145-5p are functionally related and mediate the effect of ZFPM2 on a number of its protein targets in blood including VEGFA, SERPINE1, and PDGFs. This study identifies SNPs that have a regulatory impact on circulating miRNAs, and underlines the role of several circulating miRNAs in mediating the effect of a number of GWAS loci on cardiometabolic phenotypes.

Cis-regulated additively expressed genes play a fundamental role in the formation of triploid loquat (Eriobotrya japonica (Thunb.) Lindl.) Heterosis.

Triploid loquat (Eriobotrya japonica (Thunb.) Lindl.) has greater vigor than their respective diploid and tetraploid parents, but the molecular basis of this triploid heterosis remains unclear. Recent studies have suggested that DNA methylation is involved in heterosis, which is a recognized method of suppressing gene expression. However, our previous studies revealed a trend of increased DNA methylation in triploid loquat hybrids compared to their parents. To elucidate the mechanism of triploid loquat heterosis, we investigated the levels and regulation of relative gene expression between hybrid and parental lines using RNA-Seq technology. We found that gene expression in the hybrid lines was down-regulated and gene expression analysis revealed that approximately 94.56 and 86.97% were expressed additively in triploid-A and triploid-B, respectively. Analyses of the allele-specific gene expression in the hybrids revealed significantly more Longquan-1 alleles were preferentially expressed in the two hybrid lines. Further analysis of cis- and trans-regulatory effects showed that gene expression variation between parental alleles is largely attributable to cis-regulatory variation in triploid loquat and analyses of genes belonging to cis-regulatory variation showed that 88-90% of cis genes contributed to an additive expression pattern. Taken together, our results suggest that gene expression variation in triploid loquat fundamentally cis-regulated may play a dominant role in triploid loquat heterosis.

Expression analysis of candidate genes for fatty acid composition in adipose tissue and identification of regulatory regions

The aim of this work was to study the genetic basis of the backfat expression of lipid-related genes associated with meat quality traits in pigs. We performed a genome-wide association study with the backfat gene expression measured in 44 genes by qPCR and the PorcineSNP60 BeadChip genotypes in 115 Iberian x Landrace backcross animals. A total of 193 expression-associated SNPs located in 19 chromosomal regions were associated with expression levels of ACSM5, ELOVL6, FABP4, FADS2, and SLC27A4 genes. Three expression quantitative trail loci (eQTLs) corresponding to ACSM5, FABP4, and FADS2 were classified as cis-acting eQTLs, whereas the remaining 16 eQTLs have trans-regulatory effects. Remarkably, a SNP in the ACSM5 promoter region and a SNP in the 3′UTR region of FABP4 were the most associated polymorphisms with the ACSM5 and FABP4 expression levels, respectively. Moreover, relevant lipid-related genes mapped in the trans-eQTLs regions associated with the ACSM5, FABP4, FADS2, and SLC27A4 genes. Interestingly, a trans-eQTL hotspot on SSC13 regulating the gene expression of ELOVL6, ELOLV5, and SCD, three important genes implicated in the elongation and desaturation of fatty acids, was identified. These findings provide new data to further understand the functional regulatory mechanisms implicated in the variation of fatty acid composition in pigs.