Research progress in allele-specific expression and its regulatory mechanisms

Evolution and genetics of accessory gland transcriptome divergence between Drosophila melanogaster and D. simulans

Studies of allele-specific expression in interspecific hybrids have provided important insights into gene-regulatory divergence and hybrid incompatibilities. Many such investigations in Drosophila have used transcriptome data from complex mixtures of many tissues or from gonads, however, regulatory divergence may vary widely among species, sexes, and tissues. Thus, we lack sufficiently broad sampling to be confident about the general biological principles of regulatory divergence. Here, we seek to fill some of these gaps in the literature by characterizing regulatory evolution and hybrid misexpression in a somatic male sex organ, the accessory gland, in F1 hybrids between Drosophila melanogaster and D. simulans. The accessory gland produces seminal fluid proteins, which play an important role in male and female fertility and may be subject to adaptive divergence due to male-male or male-female interactions. We find that trans differences are relatively more abundant than cis, in contrast to most of the interspecific hybrid literature, though large effect-size trans differences are rare. Seminal fluid protein genes have significantly elevated levels of expression divergence and tend to be regulated through both cis and trans divergence. We find limited misexpression (over- or underexpression relative to both parents) in this organ compared to most other Drosophila studies. As in previous studies, male-biased genes are overrepresented among misexpressed genes and are much more likely to be underexpressed. ATAC-Seq data show that chromatin accessibility is correlated with expression differences among species and hybrid allele-specific expression. This work identifies unique regulatory evolution and hybrid misexpression properties of the accessory gland and suggests the importance of tissue-specific allele-specific expression studies.

Systematic evaluation of single-cell RNA-seq analyses performance based on long-read sequencing platforms

INTRODUCTION

The rapid development of next-generation sequencing (NGS)-based single-cell RNA sequencing (scRNA-seq) allows for detecting and quantifying gene expression in a high-throughput manner, providing a powerful tool for comprehensively understanding cellular function in various biological processes. However, the NGS-based scRNA-seq only quantifies gene expression and cannot reveal the exact transcript structures (isoforms) of each gene due to the limited read length. On the other hand, the long read length of third-generation sequencing (TGS) technologies, including Oxford Nanopore Technologies (ONT) and Pacific Biosciences (PacBio), enable direct reading of intact cDNA molecules.

OBJECTIVES

Both ONT and PacBio have been used in conjunction with scRNA-seq, but their performance in single-cell analyses has not been systematically evaluated.

METHODS

To address this, we generated ONT and PacBio data from the same single-cell cDNA libraries containing different amount of cells.

RESULTS

Using NGS as a control, we assessed the performance of each platform in cell type identification. Additionally, the reliability in identifying novel isoforms and allele-specific gene/isoform expression by both platforms was verified, providing a systematic evaluation to design the sequencing strategies in single-cell transcriptome studies.

CONCLUSION

Beyond gene expression analysis, which the NGS-based scRNA-seq only affords, TGS-based scRNA-seq achieved gene splicing analyses, identifying novel isoforms. Attribute to higher sequencing quality of PacBio, it outperforms ONT in accuracy of novel transcripts identification and allele-specific gene/isoform expression.

High-quality genome assembly enables prediction of allele-specific gene expression in hybrid poplar

Poplar (Populus) is a well-established model system for tree genomics and molecular breeding, and hybrid poplar is widely used in forest plantations. However, distinguishing its diploid homologous chromosomes is difficult, complicating advanced functional studies on specific alleles. In this study, we applied a trio-binning design and PacBio high-fidelity long-read sequencing to obtain haplotype-phased telomere-to-telomere genome assemblies for the 2 parents of the well-studied F1 hybrid "84K" (Populus alba × Populus tremula var. glandulosa). Almost all chromosomes, including the telomeres and centromeres, were completely assembled for each haplotype subgenome apart from 2 small gaps on one chromosome. By incorporating information from these haplotype assemblies and extensive RNA-seq data, we analyzed gene expression patterns between the 2 subgenomes and alleles. Transcription bias at the subgenome level was not uncovered, but extensive-expression differences were detected between alleles. We developed machine-learning (ML) models to predict allele-specific expression (ASE) with high accuracy and identified underlying genome features most highly influencing ASE. One of our models with 15 predictor variables achieved 77% accuracy on the training set and 74% accuracy on the testing set. ML models identified gene body CHG methylation, sequence divergence, and transposon occupancy both upstream and downstream of alleles as important factors for ASE. Our haplotype-phased genome assemblies and ML strategy highlight an avenue for functional studies in Populus and provide additional tools for studying ASE and heterosis in hybrids.

Patterns of Gene Expression, Splicing, and Allele-Specific Expression Vary among Macular Tissues and Clinical Stages of Age-Related Macular Degeneration

Age-related macular degeneration (AMD) is a leading cause of blindness, and elucidating its underlying disease mechanisms is vital to the development of appropriate therapeutics. We identified differentially expressed genes (DEGs) and differentially spliced genes (DSGs) across the clinical stages of AMD in disease-affected tissue, the macular retina pigment epithelium (RPE)/choroid and the macular neural retina within the same eye. We utilized 27 deeply phenotyped donor eyes (recovered within a 6 h postmortem interval time) from Caucasian donors (60-94 years) using a standardized published protocol. Significant findings were then validated in an independent set of well-characterized donor eyes (n = 85). There was limited overlap between DEGs and DSGs, suggesting distinct mechanisms at play in AMD pathophysiology. A greater number of previously reported AMD loci overlapped with DSGs compared to DEGs between disease states, and no DEG overlap with previously reported loci was found in the macular retina between disease states. Additionally, we explored allele-specific expression (ASE) in coding regions of previously reported AMD risk loci, uncovering a significant imbalance in C3 rs2230199 and CFH rs1061170 in the macular RPE/choroid for normal eyes and intermediate AMD (iAMD), and for CFH rs1061147 in the macular RPE/choroid for normal eyes and iAMD, and separately neovascular AMD (NEO). Only significant DEGs/DSGs from the macular RPE/choroid were found to overlap between disease states. STAT1, validated between the iAMD vs. normal comparison, and AGTPBP1, BBS5, CERKL, FGFBP2, KIFC3, RORα, and ZNF292, validated between the NEO vs. normal comparison, revealed an intricate regulatory network with transcription factors and miRNAs identifying potential upstream and downstream regulators. Findings regarding the complement genes C3 and CFH suggest that coding variants at these loci may influence AMD development via an imbalance of gene expression in a tissue-specific manner. Our study provides crucial insights into the multifaceted genomic underpinnings of AMD (i.e., tissue-specific gene expression changes, potential splice variation, and allelic imbalance), which may open new avenues for AMD diagnostics and therapies specific to iAMD and NEO.

Haplotype-Resolution Transcriptome Analysis Reveals Important Responsive Gene Modules and Allele-Specific Expression Contributions under Continuous Salt and Drought in Camellia sinensis

The tea plant, Camellia sinensis (L.) O. Kuntze, is one of the most important beverage crops with significant economic and cultural value. Global climate change and population growth have led to increased salt and drought stress, negatively affecting tea yield and quality. The response mechanism of tea plants to these stresses remains poorly understood due to the lack of reference genome-based transcriptional descriptions. This study presents a high-quality genome-based transcriptome dynamic analysis of C. sinensis' response to salt and drought stress. A total of 2244 upregulated and 2164 downregulated genes were identified under salt and drought stress compared to the control sample. Most of the differentially expression genes (DEGs) were found to involve divergent regulation processes at different time points under stress. Some shared up- and downregulated DEGs related to secondary metabolic and photosynthetic processes, respectively. Weighted gene co-expression network analysis (WGCNA) revealed six co-expression modules significantly positively correlated with C. sinensis' response to salt or drought stress. The MEpurple module indicated crosstalk between the two stresses related to ubiquitination and the phenylpropanoid metabolic regulation process. We identified 1969 salt-responsive and 1887 drought-responsive allele-specific expression (ASE) genes in C. sinensis. Further comparison between these ASE genes and tea plant heterosis-related genes suggests that heterosis likely contributes to the adversity and stress resistance of C. sinensis. This work offers new insight into the underlying mechanisms of C. sinensis' response to salt and drought stress and supports the improved breeding of tea plants with enhanced salt and drought tolerance.

Gene expression variations and allele-specific expression of two rice and their hybrid in caryopses at single-nucleus resolution

Seeds are an indispensable part of the flowering plant life cycle and a critical determinant of agricultural production. Distinct differences in the anatomy and morphology of seeds separate monocots and dicots. Although some progress has been made with respect to understanding seed development in Arabidopsis, the transcriptomic features of monocotyledon seeds at the cellular level are much less understood. Since most important cereal crops, such as rice, maize, and wheat, are monocots, it is essential to study transcriptional differentiation and heterogeneity during seed development at a finer scale. Here, we present single-nucleus RNA sequencing (snRNA-seq) results of over three thousand nuclei from caryopses of the rice cultivars Nipponbare and 9311 and their intersubspecies F₁ hybrid. A transcriptomics atlas that covers most of the cell types present during the early developmental stage of rice caryopses was successfully constructed. Additionally, novel specific marker genes were identified for each nuclear cluster in the rice caryopsis. Moreover, with a focus on rice endosperm, the differentiation trajectory of endosperm subclusters was reconstructed to reveal the developmental process. Allele-specific expression (ASE) profiling in endosperm revealed 345 genes with ASE (ASEGs). Further pairwise comparisons of the differentially expressed genes (DEGs) in each endosperm cluster among the three rice samples demonstrated transcriptional divergence. Our research reveals differentiation in rice caryopsis from the single-nucleus perspective and provides valuable resources to facilitate clarification of the molecular mechanism underlying caryopsis development in rice and other monocots.

Genetic analysis and allele-specific expression of SMAD7 3'UTR variants in human colorectal cancer reveal a novel somatic variant exhibiting allelic imbalance

BACKGROUND

Considering the impact of SMAD7 deregulation in colorectal cancer (CRC) progression and the significance of single nucleotide variant (SNV)-mediated disruptions of microRNA (miRNA)-dependent regulation for cancer susceptibility, our study aimed to analyze genetic variation in the SMAD7 3' untranslated region ( 3'UTR) in CRC, measure differences in allelic mRNA expression, and evaluate its interference with miRNA-mediated post-transcriptional regulation.

PATIENTS AND METHODS

This study included 80 patients with different CRC stages and six human colon cancer cell lines of various histological origins. SMAD7 3'UTR was analyzed by direct sequencing, followed by the relative quantification of differential allelic expression of detected variants by allele-specific qRT-PCR. In silico tools were employed for predictions of regulatory consequences of detected variants.

RESULTS

A total of four different SNVs in one cell line and nine patients were found, among which were a novel somatic point variant and three already known germline variants (rs16950113, rs1050799536, and rs1043778717). All evaluated SNVs exhibited variable extents of allelic imbalance in expression. In silico analysis predicted significant effects of SNVs on miRNA binding efficiency, with each SNV disrupting existing and creating new target sites for one or more miRNAs.

CONCLUSION

Imbalance observed in the expression of SNV alleles altering miRNA binding suggests that all investigated SNVs are potential contributing factors impacting SMAD7 expression regulation in CRC that further studies should investigate.

Genome-wide analysis of NBS-LRR genes revealed contribution of disease resistance from Saccharum spontaneum to modern sugarcane cultivar

INTRODUCTION

During plant evolution, nucleotide-binding sites (NBS) and leucine-rich repeat (LRR) genes have made significant contributions to plant disease resistance. With many high-quality plant genomes sequenced, identification and comprehensive analyses of NBS-LRR genes at whole genome level are of great importance to understand and utilize them.

METHODS

In this study, we identified the NBS-LRR genes of 23 representative species at whole genome level, and researches on NBS-LRR genes of four monocotyledonous grass species, Saccharum spontaneum, Saccharum officinarum, Sorghum bicolor and Miscanthus sinensis, were focused.

RESULTS AND DISCUSSION

We found that whole genome duplication, gene expansion, and allele loss could be factors affecting the number of NBS-LRR genes in the species, and whole genome duplication is likely to be the main cause of the number of NBS-LRR genes in sugarcane. Meanwhile, we also found a progressive trend of positive selection on NBS-LRR genes. These studies further elucidated the evolutionary pattern of NBS-LRR genes in plants. Transcriptome data from multiple sugarcane diseases revealed that more differentially expressed NBS-LRR genes were derived from S. spontaneum than from S. officinarum in modern sugarcane cultivars, and the proportion was significantly higher than the expected. This finding reveals that S. spontaneum has a greater contribution to disease resistance for modern sugarcane cultivars. In addition, we observed allelespecific expression of seven NBS-LRR genes under leaf scald, and 125 NBS-LRR genes responding to multiple diseases were identified. Finally, we built a plant NBS-LRR gene database to facilitate subsequent analysis and use of NBSLRR genes obtained here. In conclusion, this study complemented and completed the research of plant NBS-LRR genes, and discussed how NBS-LRR genes responding to sugarcane diseases, which provided a guide and genetic resources for further research and utilization of NBS-LRR genes.

Transcriptome analysis of sugarcane reveals rapid defense response of SES208 to Xanthomonas albilineans in early infection

BACKGROUND

Diseases are the major factor affecting the quality and yield of sugarcane during its growth and development. However, our knowledge about the factors regulating disease responses remain limited. The present study focuses on identifying genes regulating transcriptional mechanisms responsible for resistance to leaf scald caused by Xanthomonas albilineans in S. spontaneum and S. officinarum.

RESULTS

After inoculation of the two sugarcane varieties SES208 (S. spontaneum) and LA Purple (S. officinarum) with Xanthomonas albilineans, SES208 exhibited significantly greater resistance to leaf scald caused by X. albilineans than did LA Purple. Using transcriptome analysis, we identified a total of 4323 and 1755 differentially expressed genes (DEGs) in inoculated samples of SES208 and LA Purple, respectively. Significantly, 262 DEGs were specifically identified in SES208 that were enriched for KEGG pathway terms such as plant-pathogen interaction, MAPK signaling pathway, and plant hormone signal transduction. Furthermore, we built a transcriptional regulatory co-expression network that specifically identified 16 and 25 hub genes in SES208 that were enriched for putative functions in plant-pathogen interactions, MAPK signaling, and plant hormone signal transduction. All of these essential genes might be significantly involved in resistance-regulating responses in SES208 after X. albilineans inoculation. In addition, we found allele-specific expression in SES208 that was associated with the resistance phenotype of SES208 when infected by X. albilineans. After infection with X. albilineans, a great number of DEGs associated with the KEGG pathways 'phenylpropanoid biosynthesis' and 'flavonoid biosynthesis' exhibited significant expression changes in SES208 compared to LA Purple that might contribute to superior leaf scald resistance in SES208.

CONCLUSIONS

We provided the first systematical transcriptome map that the higher resistance of SES208 is associated with and elicited by the rapid activation of multiple clusters of defense response genes after infection by X. albilineans and not merely due to changes in the expression of genes generically associated with stress resistance. These results will serve as the foundation for further understanding of the molecular mechanisms of resistance against X. albilineans in S. spontaneum.

Allele-specific expression and chromatin accessibility contribute to heterosis in tea plants (Camellia sinensis)

Heterosis is extensively used to improve crop productivity, yet its allelic and chromatin regulation remains unclear. Based on our resolved genomes of the maternal TGY and paternal HD, we analyzed the contribution of allele-specific expression (ASE) and chromatin accessibility of JGY and HGY, the artificial hybrids of oolong tea with the largest cultivated area in China. The ASE genes (ASEGs) of tea hybrids with maternal-biased were mainly related to the energy and terpenoid metabolism pathways, whereas the ASEGs with paternal-biased tend to be enriched in glutathione metabolism, and these parental bias of hybrids may coordinate and lead to the acquisition of heterosis in more biological pathways. ATAC-seq results showed that hybrids have significantly higher accessible chromatin regions (ACRs) compared with their parents, which may confer broader and stronger transcriptional activity of genes in hybrids. The number of ACRs with significantly increased accessibility in hybrids was much greater than decreased, and the associated alleles were also affected by differential ACRs across different parents, suggesting enhanced positive chromatin regulation and potential genetic effects in hybrids. Core ASEGs of terpene and purine alkaloid metabolism pathways with significant positive heterosis have greater chromatin accessibility in hybrids, and were potentially regulated by several members of the MYB, DOF and TRB families. The binding motif of CsMYB85 in the promoter ACR of the rate-limiting enzyme CsDXS was verified by DAP-seq. These results suggest that higher numbers and more accessible ACRs in hybrids contribute to the regulation of ASEGs, thereby affecting the formation of heterotic metabolites.

Comparative transcriptomic analysis of maize ear heterosis during the inflorescence meristem differentiation stage

BACKGROUND

Heterosis is widely used in many crops and is important for global food safety, and maize is one of the most successful crops to take advantage of heterosis. Gene expression patterns control the development of the maize ear, but the mechanisms by which heterosis affects transcriptional-level control are not fully understood.

RESULTS

In this study, we sampled ear inflorescence meristems (IMs) from the single-segment substitution maize (Zea mays) line lx9801^hlEW2b, which contains the heterotic locus hlEW2b associated with ear width, as well as the receptor parent lx9801, the test parent Zheng58, and their corresponding hybrids Zheng58 × lx9801^hlEW2b (HY) and Zheng58 × lx9801 (CK). After RNA sequencing and transcriptomic analysis, 2531 unique differentially expressed genes (DEGs) were identified between the two hybrids (HY vs. CK). Our results showed that approximately 64% and 48% of DEGs exhibited additive expression in HY and CK, whereas the other genes displayed a non-additive expression pattern. The DEGs were significantly enriched in GO functional categories of multiple metabolic processes, plant organ morphogenesis, and hormone regulation. These essential processes are potentially associated with heterosis performance during the maize ear developmental stage. In particular, 125 and 100 DEGs from hybrids with allele-specific expression (ASE) were specifically identified in HY and CK, respectively. Comparison between the two hybrids suggested that ASE genes were involved in different development-related processes that may lead to the hybrid vigor phenotype during maize ear development. In addition, several critical genes involved in auxin metabolism and IM development were differentially expressed between the hybrids and showed various expression patterns (additive, non-additive, and ASE). Changes in the expression levels of these genes may lead to differences in auxin homeostasis in the IM, affecting the transcription of core genes such as WUS that control IM development.

CONCLUSIONS

Our research suggests that additive, non-additive, and allele-specific expression patterns may fine-tune the expression of crucial DEGs that modulate carbohydrate and protein metabolic processes, nitrogen assimilation, and auxin metabolism to optimal levels, and these transcriptional changes may play important roles in maize ear heterosis. The results provide new information that increases our understanding of the relationship between transcriptional variation and heterosis during maize ear development, which may be helpful for clarifying the genetic and molecular mechanisms of heterosis.

Transposon insertions regulate genome-wide allele-specific expression and underpin flower colour variations in apple (Malus spp.)

Allele-specific expression (ASE) can lead to phenotypic diversity and evolution. However, the mechanisms regulating ASE are not well understood, particularly in woody perennial plants. In this study, we investigated ASE genes in the apple cultivar 'Royal Gala' (RG). A high quality chromosome-level genome was assembled using a homozygous tetra-haploid RG plant, derived from anther cultures. Using RNA-sequencing (RNA-seq) data from RG flower and fruit tissues, we identified 2091 ASE genes. Compared with the haploid genome of 'Golden Delicious' (GD), a parent of RG, we distinguished the genomic sequences between the two alleles of 817 ASE genes, and further identified allele-specific presence of a transposable element (TE) in the upstream region of 354 ASE genes. These included MYB110a that encodes a transcription factor regulating anthocyanin biosynthesis. Interestingly, another ASE gene, MYB10 also showed an allele-specific TE insertion and was identified using genome data of other apple cultivars. The presence of the TE insertion in both MYB genes was positively associated with ASE and anthocyanin accumulation in apple petals through analysis of 231 apple accessions, and thus underpins apple flower colour evolution. Our study demonstrated the importance of TEs in regulating ASE on a genome-wide scale and presents a novel method for rapid identification of ASE genes and their regulatory elements in plants.

Drosophila as a Model System for Studying of the Evolution and Functional Specialization of the Y Chromosome

The Y chromosome is one of the sex chromosomes found in males of animals of different taxa, including insects and mammals. Among all chromosomes, the Y chromosome is characterized by a unique chromatin landscape undergoing dynamic evolutionary change. Being entirely heterochromatic, the Y chromosome as a rule preserves few functional genes, but is enriched in tandem repeats and transposons. Due to difficulties in the assembly of the highly repetitive Y chromosome sequence, deep analyses of Y chromosome evolution, structure, and functions are limited to a few species, one of them being Drosophila melanogaster. Despite Y chromosomes exhibiting high structural divergence between even closely related species, Y-linked genes have evolved convergently and are mainly associated with spermatogenesis-related activities. This indicates that male-specific selection is a dominant force shaping evolution of Y chromosomes across species. This review presents our analysis of current knowledge concerning Y chromosome functions, focusing on recent findings in Drosophila. Here we dissect the experimental and bioinformatics data about the Y chromosome accumulated to date in Drosophila species, providing comparative analysis with mammals, and discussing the relevance of our analysis to a wide range of eukaryotic organisms, including humans.

The multi-omics basis of potato heterosis

Heterosis is a fundamental biological phenomenon characterized by the superior performance of hybrids over their parents. Although tremendous progress has been reported in seed crops, the molecular mechanisms underlying heterosis in clonally propagated crops are largely unknown. Potato (Solanum tuberosum L.) is the most important tuber crop and an ongoing revolution is transforming potato from a clonally propagated tetraploid crop into a seed-propagated diploid hybrid potato. In our previous study, we developed the first generation of highly homozygous inbred lines of potato and hybrids with strong heterosis. Here, we integrated transcriptome, metabolome, and DNA methylation data to explore the genetic and molecular basis of potato heterosis at three developmental stages. We found that the initial establishment of heterosis in diploid potato was mainly due to dominant complementation. Flower color, male fertility, and starch and sucrose metabolism showed obvious gene dominant complementation in hybrids, and hybrids devoted more energy to primary metabolism for rapid growth. In addition, we identified ~2 700 allele-specific expression genes at each stage, which likely function in potato heterosis and might be regulated by CHH allele-specific methylation level. Our multi-omics analysis provides insight into heterosis in potato and facilitates the exploitation of heterosis in potato breeding.

Chemerin Impact on Alternative mRNA Transcription in the Porcine Luteal Cells

Chemerin participates in the regulation of processes related to physiological and disorder mechanisms in mammals, including metabolism, obesity, inflammation, and reproduction. In this study, we have investigated chemerin influence on alternative mRNA transcription within the porcine luteal cell transcriptome, such as differential expression of long non-coding RNAs (DELs) and their interactions with differentially expressed genes (DEGs), differences in alternative splicing of transcripts (DASs), and allele-specific expression (ASEs) related to the single nucleotide variants (SNVs) frequency. Luteal cells were collected from gilts during the mid-luteal phase of the oestrous cycle. After in vitro culture of cells un-/treated with chemerin, the total RNA was isolated and sequenced using the high-throughput method. The in silico analyses revealed 24 DELs cis interacting with 6 DEGs and trans-correlated with 300 DEGs, 137 DASs events, and 18 ASEs. The results enabled us to analyse metabolic and signalling pathways in detail, providing new insights into the effects of chemerin on the corpus luteum functions related to inflammatory response, leukocyte infiltration, the occurrence of luteotropic and luteolytic signals (leading to apoptosis and/or necroptosis). Validation of the results using qPCR confirmed the predicted expression changes. Chemerin at physiological concentrations significantly modifies the transcription processes in the porcine luteal cells.

Gene age shapes the transcriptional landscape of sexual morphogenesis in mushroom forming fungi (Agaricomycetes)

Multicellularity has been one of the most important innovations in the history of life. The role of gene regulatory changes in driving transitions to multicellularity is being increasingly recognized; however, factors influencing gene expression patterns are poorly known in many clades. Here, we compared the developmental transcriptomes of complex multicellular fruiting bodies of eight Agaricomycetes and Cryptococcus neoformans, a closely related human pathogen with a simple morphology. In-depth analysis in Pleurotus ostreatus revealed that allele-specific expression, natural antisense transcripts, and developmental gene expression, but not RNA editing or a ‘developmental hourglass,’ act in concert to shape its transcriptome during fruiting body development. We found that transcriptional patterns of genes strongly depend on their evolutionary ages. Young genes showed more developmental and allele-specific expression variation, possibly because of weaker evolutionary constraint, suggestive of nonadaptive expression variance in fruiting bodies. These results prompted us to define a set of conserved genes specifically regulated only during complex morphogenesis by excluding young genes and accounting for deeply conserved ones shared with species showing simple sexual development. Analysis of the resulting gene set revealed evolutionary and functional associations with complex multicellularity, which allowed us to speculate they are involved in complex multicellular morphogenesis of mushroom fruiting bodies.

Limited allele-specific gene expression in highly polyploid sugarcane

Polyploidy is widespread in plants, allowing the different copies of genes to be expressed differently in a tissue-specific or developmentally specific way. This allele-specific expression (ASE) has been widely reported, but the proportion and nature of genes showing this characteristic have not been well defined. We now report an analysis of the frequency and patterns of ASE at the whole-genome level in the highly polyploid sugarcane genome. Very high depth whole-genome sequencing and RNA sequencing revealed strong correlations between allelic proportions in the genome and in expressed sequences. This level of sequencing allowed discrimination of each of the possible allele doses in this 12-ploid genome. Most genes were expressed in direct proportion to the frequency of the allele in the genome with examples of polymorphisms being found with every possible discrete level of dose from 1:11 for single-copy alleles to 12:0 for monomorphic sites. The rarer cases of ASE were more frequent in the expression of defense-response genes, as well as in some processes related to the biosynthesis of cell walls. ASE was more common in genes with variants that resulted in significant disruption of function. The low level of ASE may reflect the recent origin of polyploid hybrid sugarcane. Much of the ASE present can be attributed to strong selection for resistance to diseases in both nature and domestication.

Variation burst during dedifferentiation and increased CHH-type DNA methylation after 30 years of in vitro culture of sweet orange

Somaclonal variation arising from tissue culture may provide a valuable resource for the selection of new germplasm, but may not preserve true-to-type characteristics, which is a major concern for germplasm conservation or genome editing. The genomic changes associated with dedifferentiation and somaclonal variation during long-term in vitro culture are largely unknown. Sweet orange was one of the earliest plant species to be cultured in vitro and induced via somatic embryogenesis. We compared four sweet orange callus lines after 30 years of constant tissue culture with newly induced calli by comprehensively determining the single-nucleotide polymorphisms, copy number variations, transposable element insertions, methylomic and transcriptomic changes. We identified a burst of variation during early dedifferentiation, including a retrotransposon outbreak, followed by a variation purge during long-term in vitro culture. Notably, CHH methylation showed a dynamic pattern, initially disappearing during dedifferentiation and then more than recovering after 30 years of in vitro culture. We also analyzed the effects of somaclonal variation on transcriptional reprogramming, and indicated subgenome dominance was evident in the tetraploid callus. We identified a retrotransposon insertion and DNA modification alternations in the potential regeneration-related gene CLAVATA3/EMBRYO SURROUNDING REGION-RELATED 16. This study provides the foundation to harness in vitro variation and offers a deeper understanding of the variation introduced by tissue culture during germplasm conservation, somatic embryogenesis, gene editing, and breeding programs.

From Cave Dragons to Genomics: Advancements in the Study of Subterranean Tetrapods

Throughout most of the kingdom Animalia, evolutionary transitions from surface life to a life permanently bound to caves and other subterranean habitats have occurred innumerous times. Not so in tetrapods, where a mere 14 cave-obligate species—all plethodontid and proteid salamanders—are known. We discuss why cave tetrapods are so exceptional and why only salamanders have made the transition. Their evolution follows predictable and convergent, albeit independent pathways. Among the many known changes associated with transitions to subterranean life, eye degeneration, starvation resistance, and longevity are especially relevant to human biomedical research. Recently, sequences of salamander genomes have become available opening up genomic research for cave tetrapods. We discuss new genomic methods that can spur our understanding of the evolutionary mechanisms behind convergent phenotypic change, the relative roles of selective and neutral evolution, cryptic species diversity, and data relevant for conservation such as effective population size and demography.

The integrative biology of genetic dominance

Dominance is a basic property of inheritance systems describing the link between a diploid genotype at a single locus and the resulting phenotype. Models for the evolution of dominance have long been framed as an opposition between the irreconcilable views of Fisher in 1928 supporting the role of largely elusive dominance modifiers and Wright in 1929, who viewed dominance as an emerging property of the structure of enzymatic pathways. Recent theoretical and empirical advances however suggest that these opposing views can be reconciled, notably using models investigating the regulation of gene expression and developmental processes. In this more comprehensive framework, phenotypic dominance emerges from departures from linearity between any levels of integration in the genotype‐to‐phenotype map. Here, we review how these different models illuminate the emergence and evolution of dominance. We then detail recent empirical studies shedding new light on the diversity of molecular and physiological mechanisms underlying dominance and its evolution. By reconciling population genetics and functional biology, we hope our review will facilitate cross‐talk among research fields in the integrative study of dominance evolution.

Allele Specific Expression (ASE) analysis between Bos Taurus and Bos Indicus cows using RNA-Seq data at SNP level and gene level

In the current study, allele specific expression analysis was performed in two subspecies cows (Bos taurus and Bos indicus) at SNP and gene levels. RNA-Seq data of 21,078,477 and 20940063 paired end reads from pooling of whole blood samples (Leukocyte) from 40 US Holstein (Bos Taurus) and 45 Cholistani cows (Bos indicus) obtained from SRA database in NCBI. Quality control and trimming of row RNA-Seq data were processed by FASTQC and Trimmomatic softwares. The transcriptome was assembled by TopHat2 software in two cow's population by aligning and mapping the RNA-Seq reads on bovine reference genome. The SNPs were discovered by Samtools software and ASE analysis was performed by Chi-square test. Results showed that 50183 and 137954 SNPs were discovered on the assembled transcriptome of Holstein and Cholistani cow samples, respectively, and 15308 SNPs were common in both breeds. 10158 SNPs from 50183 (20%) in Holstein and 31523 SNPs from 137954 (23%) in Cholistani cows were identified as ASE-SNPs. Reference allele and alternative allele count in Holstein and Cholistani cows were 3041 and 7155, respectively. Among 131 discovered SNPs in 41 genes with different expression in Holstein and Cholistani cows, 31 ASE-SNPs (5 in Holstein; 26 in Cholistani cows) were discovered.

Biased allelic expression in tissues of F1 hybrids between tropical and temperate lotus (Nelumbo nuicfera)

The genome-wide allele-specific expression in F1 hybrids from the cross of tropical and temperate lotus unveils how cis-regulatory divergences affect genes in key pathways related to ecotypic divergence. Genetic variation, particularly cis-regulatory variation, plays a crucial role in phenotypic variation and adaptive evolution in plants. Temperate and tropical lotus, the two ecotypes of Nelumbo nucifera, show distinction in the degree of rhizome enlargement, which is associated with winter dormancy. To understand the roles of genome-wide cis-regulatory divergences on adaptive evolution of temperate and tropical lotus (Nelumbo nucifera), here we performed allele-specific expression (ASE) analyses on the tissues including flowers, leaves and rhizome from F1 hybrids of tropical and temperate lotus. For all investigated tissues in F1s, about 36% of genes showed ASE and about 3% of genes showed strong consistent ASE. Most of ASEs were biased towards the tropical parent in all surveyed samples, indicating that the tropical genome might be dominant over the temperate genome in gene expression of tissues from their F1 hybrids. We found that promoter sequences with similar allelic expression are more conserved than genes with significant or conditional ASE, suggesting the cis-regulatory sequence divergence underlie the allelic expression bias. We further uncovered biased genes being related to phenotypic differentiation between two lotus ecotypes, especially metabolic and phytohormone-related pathways in the rhizome. Overall, our study provides a global landscape of cis-regulatory variations between two lotus ecotypes and highlights their roles in rhizome growth variation for the climatic adaptation.

Exploring allele specific methylation in drug dependence susceptibility

Drug dependence is a neuropsychiatric condition that involves genetic, epigenetic and environmental factors. Allele-specific methylation (ASM) is a common and stable epigenetic mechanism that involves genetic variants correlating with differential levels of methylation at CpG sites. We selected 182 single-nucleotide polymorphisms (SNPs) described to influence cis ASM in human brain regions to evaluate their possible contribution to drug dependence susceptibility. We performed a case-control association study in a discovery sample of 578 drug-dependent patients (including 428 cocaine-dependent subjects) and 656 controls from Spain, and then, we followed-up the significant associations in an independent sample of 1119 cases (including 589 cocaine-dependent subjects) and 1092 controls. In the discovery sample, we identified five nominal associations, one of them replicated in the follow-up sample (rs6020251). The pooled analysis revealed an association between drug dependence and rs6020251 but also rs11585570, both overcoming the Bonferroni correction for multiple testing. We performed the same analysis considering only cocaine-dependent patients and obtained similar results. The rs6020251 variant correlates with differential methylation levels of cg17974185 and lies in the first intron of the CTNNBL1 gene, in a genomic region with multiple histone marks related to enhancer and promoter regions in brain. Rs11585570 is an eQTL in brain and blood for the SCP2 and ECHDC2 genes and correlates with differential methylation of cg27535305 and cg13461509, located in the promoter regions of both genes. To conclude, using an approach that combines genetic and epigenetic data, we highlighted the CTNNBL1, SCP2 and ECHDC2 genes as potential contributors to drug dependence susceptibility.

Global transcriptional profiling between inbred parents and hybrids provides comprehensive insights into ear-length heterosis of maize (Zea mays)

BACKGROUND

Maize (Zea mays) ear length, which is an important yield component, exhibits strong heterosis. Understanding the potential molecular mechanisms of ear-length heterosis is critical for efficient yield-related breeding.

RESULTS

Here, a joint netted pattern, including six parent-hybrid triplets, was designed on the basis of two maize lines harboring long (T121 line) and short (T126 line) ears. Global transcriptional profiling of young ears (containing meristem) was performed. Multiple comparative analyses revealed that 874 differentially expressed genes are mainly responsible for the ear-length variation between T121 and T126 lines. Among them, four key genes, Zm00001d049958, Zm00001d027359, Zm00001d048502 and Zm00001d052138, were identified as being related to meristem development, which corroborated their roles in the superior additive genetic effects on ear length in T121 line. Non-additive expression patterns were used to identify candidate genes related to ear-length heterosis. A non-additively expressed gene (Zm00001d050649) was associated with the timing of meristematic phase transition and was determined to be the homolog of tomato SELF PRUNING, which assists SINGLE FLOWER TRUSS in driving yield-related heterosis, indicating that Zm00001d050649 is a potential contributor to drive heterotic effect on ear length.

CONCLUSION

Our results suggest that inbred parents provide genetic and heterotic effects on the ear lengths of their corresponding F₁ hybrids through two independent pathways. These findings provide comprehensive insights into the transcriptional regulation of ear length and improve the understanding of ear-length heterosis in maize.

Inducible uniparental chromosome disomy to probe genomic imprinting at single-cell level in brain and beyond

Genomic imprinting is an epigenetic mechanism that results in parental allele-specific expression of ~1% of all genes in mouse and human. Imprinted genes are key developmental regulators and play pivotal roles in many biological processes such as nutrient transfer from the mother to offspring and neuronal development. Imprinted genes are also involved in human disease, including neurodevelopmental disorders, and often occur in clusters that are regulated by a common imprint control region (ICR). In extra-embryonic tissues ICRs can act over large distances, with the largest surrounding Igf2r spanning over 10 million base-pairs. Besides classical imprinted expression that shows near exclusive maternal or paternal expression, widespread biased imprinted expression has been identified mainly in brain. In this review we discuss recent developments mapping cell type specific imprinted expression in extra-embryonic tissues and neocortex in the mouse. We highlight the advantages of using an inducible uniparental chromosome disomy (UPD) system to generate cells carrying either two maternal or two paternal copies of a specific chromosome to analyze the functional consequences of genomic imprinting. Mosaic Analysis with Double Markers (MADM) allows fluorescent labeling and concomitant induction of UPD sparsely in specific cell types, and thus to over-express or suppress all imprinted genes on that chromosome. To illustrate the utility of this technique, we explain how MADM-induced UPD revealed new insights about the function of the well-studied Cdkn1c imprinted gene, and how MADM-induced UPDs led to identification of highly cell type specific phenotypes related to perturbed imprinted expression in the mouse neocortex. Finally, we give an outlook on how MADM could be used to probe cell type specific imprinted expression in other tissues in mouse, particularly in extra-embryonic tissues.

Allele-specific expression of Parkinson's disease susceptibility genes in human brain

Genome-wide association studies have identified genetic variation in genomic loci associated with susceptibility to Parkinson’s disease (PD), the most common neurodegenerative movement disorder worldwide. We used allelic expression profiling of genes located within PD-associated loci to identify cis-regulatory variation affecting gene expression. DNA and RNA were extracted from post-mortem superior frontal gyrus tissue and whole blood samples from PD patients and controls. The relative allelic expression of transcribed SNPs in 12 GWAS risk genes was analysed by real-time qPCR. Allele-specific expression was identified for 9 out of 12 genes tested (GBA, TMEM175, RAB7L1, NUCKS1, MCCC1, BCKDK, ZNF646, LZTS3, and WDHD1) in brain tissue samples. Three genes (GPNMB, STK39 and SIPA1L2) did not show significant allele-specific effects. Allele-specific effects were confirmed in whole blood for three genes (BCKDK, LZTS3 and MCCC1), whereas two genes (RAB7L1 and NUCKS1) showed brain-specific allelic expression. Our study supports the hypothesis that changes to the cis-regulation of gene expression is a major mechanism behind a large proportion of genetic associations in PD. Interestingly, allele-specific expression was also observed for coding variants believed to be causal variants (GBA and TMEM175), indicating that splicing and other regulatory mechanisms may be involved in disease development.

Natural Selection on Exonic SNPs Shapes Allelic Expression Imbalance (AEI) Adaptability in Lung Cancer Progression

Tumors are driven by a sequence of genetic and epigenetic alterations. Previous studies have mostly focused on the roles of somatic mutations in tumorigenesis, but how germline variants act is largely unknown. In this study, we hypothesized that allelic expression imbalance (AEI) participated in the process of germline variants on tumorigenesis. We screened single-nucleotide polymorphisms (SNPs) as representative germline variants. By using 127 patients’ RNA sequencing data from paired lung cancer and adjacent normal tissues from public databases, we analyzed the effects of the functional consequence of SNPs, function and conservativeness on genes with AEI. We found that natural selection can affect AEI. Functional adaptability of genes with a high frequency of AEI and a correlation of the incidence of AEI with conservativeness were observed in both adjacent tissues and tumor tissues. Moreover, we observed a higher incidence of AEI in genes with non-synonymous SNPs than in those with synonymous SNPs. However, we also found that AEI was affected by allele expression noise, especially in tumor tissues, which led to an increased proportion of AEI, weakened the effect of natural selection and eliminated the influence of the functional consequence of SNPs on AEI. We unveiled a previously unknown adaptive regulatory mechanism in which the effect of natural selection on SNPs can be reflected in allelic expression, which provides insight into a better understanding of cancer evolution.

Unzipping haplotypes in diploid and polyploid genomes

Diploid genomes consist of two homologous copies of chromosomes with one from each parent while polyploid genomes contain more than two homologous sets of chromosomes. Most of the reference genome assemblies collapsed haplotypes that represent 'mosaic' sequences, ignoring allelic variants that may be involved in important cellular and biological functions. Unzipping haplotypes into distinct sets of sequences has been a growing trend in recent genome studies, as it is an essential tool towards resolving important clinical and biological questions, such as compound heterozygotes, heterosis, and evolution. Herein, we review existing methods for alignment-based and assembly-based haplotype phasing for heterozygous diploid and polyploid genomes, as well as recent advances of experimental approaches for improved genome phasing. We anticipate that full haplotype phasing could become a routine procedure in genome studies in the near future.

Embryo microinjection of the lecithotrophic sea urchin Heliocidaris erythrogramma

Microinjection is a common embryological technique used for many types of experiments, including lineage tracing, manipulating gene expression, or genome editing. Injectable reagents include mRNA overexpression, mis-expression, or dominant-negative experiments to examine a gene of interest, a morpholino antisense oligo to prevent translation of an mRNA or spliceoform of interest and CRISPR-Cas9 reagents. Thus, the technique is broadly useful for basic embryological studies, constructing gene regulatory networks, and directly testing hypotheses about cis-regulatory and coding sequence changes underlying the evolution of development. However, the methods for microinjection in typical planktotrophic marine invertebrates may not work well in the highly modified eggs and embryos of lecithotrophic species. This protocol is optimized for the lecithotrophic sea urchin Heliocidaris erythrogramma.

Exploring genetic variation that influences brain methylation in attention-deficit/hyperactivity disorder

Attention-deficit/hyperactivity disorder (ADHD) is a neurodevelopmental disorder caused by an interplay of genetic and environmental factors. Epigenetics is crucial to lasting changes in gene expression in the brain. Recent studies suggest a role for DNA methylation in ADHD. We explored the contribution to ADHD of allele-specific methylation (ASM), an epigenetic mechanism that involves SNPs correlating with differential levels of DNA methylation at CpG sites. We selected 3896 tagSNPs reported to influence methylation in human brain regions and performed a case-control association study using the summary statistics from the largest GWAS meta-analysis of ADHD, comprising 20,183 cases and 35,191 controls. We observed that genetic risk variants for ADHD are enriched in ASM SNPs and identified associations with eight tagSNPs that were significant at a 5% false discovery rate (FDR). These SNPs correlated with methylation of CpG sites lying in the promoter regions of six genes. Since methylation may affect gene expression, we inspected these ASM SNPs together with 52 ASM SNPs in high LD with them for eQTLs in brain tissues and observed that the expression of three of those genes was affected by them. ADHD risk alleles correlated with increased expression (and decreased methylation) of ARTN and PIDD1 and with a decreased expression (and increased methylation) of C2orf82. Furthermore, these three genes were predicted to have altered expression in ADHD, and genetic variants in C2orf82 correlated with brain volumes. In summary, we followed a systematic approach to identify risk variants for ADHD that correlated with differential cis-methylation, identifying three novel genes contributing to the disorder.

Analysis of allele-specific expression using RNA-seq of the Korean native pig and Landrace reciprocal cross

Objective

We tried to analyze allele-specific expression in the pig neocortex using bioinformatic analysis of high-throughput sequencing results from the parental genomes and offspring transcriptomes from reciprocal crosses between Korean Native and Landrace pigs.

Methods

We carried out sequencing of parental genomes and offspring transcriptomes using next generation sequencing. We subsequently carried out genome scale identification of single nucleotide polymorphisms (SNPs) in two different ways using either individual genome mapping or joint genome mapping of the same breed parents that were used for the reciprocal crosses. Using parent-specific SNPs, allele-specifically expressed genes were analyzed.

Results

Because of the low genome coverage (~4×) of the sequencing results, most SNPs were non-informative for parental lineage determination of the expressed alleles in the offspring and were thus excluded from our analysis. Consequently, 436 SNPs covering 336 genes were applicable to measure the imbalanced expression of paternal alleles in the offspring. By calculating the read ratios of parental alleles in the offspring, we identified seven genes showing allele-biased expression (p<0.05) including three previously reported and four newly identified genes in this study.

Conclusion

The newly identified allele-specifically expressing genes in the neocortex of pigs should contribute to improving our knowledge on genomic imprinting in pigs. To our knowledge, this is the first study of allelic imbalance using high throughput analysis of both parental genomes and offspring transcriptomes of the reciprocal cross in outbred animals. Our study also showed the effect of the number of informative animals on the genome level investigation of allele-specific expression using RNA-seq analysis in livestock species.

Patterns of genome-wide allele-specific expression in hybrid rice and the implications on the genetic basis of heterosis

Utilization of heterosis has greatly increased productivity of many crops globally. Allele-specific expression (ASE) has been suggested as a mechanism for causing heterosis. We performed a genome-wide analysis of ASE in three tissues of an elite rice hybrid grown under four conditions. The analysis identified 3,270 genes showing various patterns of ASE in response to developmental and environmental cues, which provides a glimpse of the ASE landscape in the hybrid genome. We showed that the ASE patterns may have distinct implications in the genetic basis of heterosis, especially in light of the classical dominance and overdominance hypotheses. The genes showing ASE provide the candidates for future studies of the genetic and molecular mechanism of heterosis.

Utilization of heterosis has greatly increased the productivity of many crops worldwide. Although tremendous progress has been made in characterizing the genetic basis of heterosis using genomic technologies, molecular mechanisms underlying the genetic components are much less understood. Allele-specific expression (ASE), or imbalance between the expression levels of two parental alleles in the hybrid, has been suggested as a mechanism of heterosis. Here, we performed a genome-wide analysis of ASE by comparing the read ratios of the parental alleles in RNA-sequencing data of an elite rice hybrid and its parents using three tissues from plants grown under four conditions. The analysis identified a total of 3,270 genes showing ASE (ASEGs) in various ways, which can be classified into two patterns: consistent ASEGs such that the ASE was biased toward one parental allele in all tissues/conditions, and inconsistent ASEGs such that ASE was found in some but not all tissues/conditions, including direction-shifting ASEGs in which the ASE was biased toward one parental allele in some tissues/conditions while toward the other parental allele in other tissues/conditions. The results suggested that these patterns may have distinct implications in the genetic basis of heterosis: The consistent ASEGs may cause partial to full dominance effects on the traits that they regulate, and direction-shifting ASEGs may cause overdominance. We also showed that ASEGs were significantly enriched in genomic regions that were differentially selected during rice breeding. These ASEGs provide an index of the genes for future pursuit of the genetic and molecular mechanism of heterosis.

Analysis of allele-specific expression of seven candidate genes involved in lipid metabolism in pig skeletal muscle and fat tissues reveals allelic imbalance of ACACA, LEP, SCD, and TNF

Analysis of allele-specific expression may help to elucidate the genetic architecture of complex traits including fat deposition in pigs. Here, we used pyrosequencing to investigate the allele proportions of candidate genes (ACACA, ADIPOR1, FASN, LEP, ME1, SCD, and TNF) involved in regulation of lipid metabolism in two fat deposits (subcutaneous and visceral fat) and longissimus dorsi muscle of pigs representing Polish Large White, Polish Landrace, Duroc, and Pietrain breeds. We detected differential allelic expression of ACACA, LEP, SCD, and TNF in all tissues analyzed. To search for putative cis-regulatory elements involved in allele-specific expression, we quantified the methylation level within CpG islands located in 5′-flanking regions of ACACA and SCD. Comparison between samples showing markedly disproportionate allelic expression and control groups with similar levels of both alleles did not reveal significant differences. We also assessed the association of rs321308225 (c.*195C>A) an SNP located in the 3′UTR of ACACA with its allelic expression in Polish Landrace pigs, but it was not significant. We conclude that allelic imbalance occurs frequently in regard to genes involved in regulation of lipid deposition in pigs, and further studies are necessary to identify cis-regulatory elements affecting ACACA, LEP, SCD, and TNF expression in porcine fat tissues and skeletal muscle.

Investigation of allele-specific expression of genes involved in adipogenesis and lipid metabolism suggests complex regulatory mechanisms of PPARGC1A expression in porcine fat tissues

Background

The expression of genes involved in regulating adipogenesis and lipid metabolism may affect economically important fatness traits in pigs. Allele-specific expression (ASE) reflects imbalance between allelic transcript levels and can be used to identify underlying cis-regulatory elements. ASE has not yet been intensively studied in pigs. The aim of this investigation was to analyze the differential allelic expression of four genes, PPARA, PPARG, SREBF1, and PPARGC1A, which are involved in the regulation of fat deposition in porcine subcutaneous and visceral fat and longissimus dorsi muscle.

Results

Quantification of allelic proportions by pyrosequencing revealed that both alleles of PPARG and SREBF1 are expressed at similar levels. PPARGC1A showed the greatest ASE imbalance in fat deposits in Polish Large White (PLW), Polish Landrace and Pietrain pigs; and PPARA in PLW pigs. Significant deviations of mean PPARGC1A allelic transcript ratio between cDNA and genomic DNA were detected in all tissues, with the most pronounced difference (p < 0.001) in visceral fat of PLW pigs. To search for potential cis-regulatory elements affecting ASE in the PPARGC1A gene we analyzed the effects of four SNPs (rs337351686, rs340650517, rs336405906 and rs345224049) in the promoter region, but none were associated with ASE in the breeds studied. DNA methylation analysis revealed significant CpG methylation differences between samples showing balanced (allelic transcript ratio ≈1) and imbalanced allelic expression for CpG site at the genomic position in chromosome 8 (SSC8): 18527678 in visceral fat (p = 0.017) and two CpG sites (SSC8:18525215, p = 0.030; SSC8:18525237, p = 0.031) in subcutaneous fat.

Conclusions

Our analysis of differential allelic expression suggests that PPARGC1A is subjected to cis-regulation in porcine fat tissues. Further studies are necessary to identify other regulatory elements localized outside the PPARGC1A proximal promoter region.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0696-6) contains supplementary material, which is available to authorized users.

Allelic imbalance and haploinsufficiency in MYBPC3-linked hypertrophic cardiomyopathy

Mutations in cardiac myosin binding protein C (MYBPC3) represent the most frequent cause of familial hypertrophic cardiomyopathy (HCM), making up approximately 50% of identified HCM mutations. MYBPC3 is distinct among other sarcomere genes associated with HCM in that truncating mutations make up the vast majority, whereas nontruncating mutations predominant in other sarcomere genes. Several studies using myocardial tissue from HCM patients have found reduced abundance of wild-type MYBPC3 compared to control hearts, suggesting haploinsufficiency of full-length MYBPC3. Further, decreased mutant versus wild-type mRNA and lack of truncated mutant MYBPC3 protein has been demonstrated, highlighting the presence of allelic imbalance. In this review, we will begin by introducing allelic imbalance and haploinsufficiency, highlighting the broad role each plays within the spectrum of human disease. We will subsequently focus on the roles allelic imbalance and haploinsufficiency play within MYBPC3-linked HCM. Finally, we will explore the implications of these findings on future directions of HCM research. An improved understanding of allelic imbalance and haploinsufficiency may help us better understand genotype-phenotype relationships in HCM and develop novel targeted therapies, providing exciting future research opportunities.

Bayesian Inference of Allele-Specific Gene Expression Indicates Abundant Cis-Regulatory Variation in Natural Flycatcher Populations

Polymorphism in cis-regulatory sequences can lead to different levels of expression for the two alleles of a gene, providing a starting point for the evolution of gene expression. Little is known about the genome-wide abundance of genetic variation in gene regulation in natural populations but analysis of allele-specific expression (ASE) provides a means for investigating such variation. We performed RNA-seq of multiple tissues from population samples of two closely related flycatcher species and developed a Bayesian algorithm that maximizes data usage by borrowing information from the whole data set and combines several SNPs per transcript to detect ASE. Of 2,576 transcripts analyzed in collared flycatcher, ASE was detected in 185 (7.2%) and a similar frequency was seen in the pied flycatcher. Transcripts with statistically significant ASE commonly showed the major allele in >90% of the reads, reflecting that power was highest when expression was heavily biased toward one of the alleles. This would suggest that the observed frequencies of ASE likely are underestimates. The proportion of ASE transcripts varied among tissues, being lowest in testis and highest in muscle. Individuals often showed ASE of particular transcripts in more than one tissue (73.4%), consistent with a genetic basis for regulation of gene expression. The results suggest that genetic variation in regulatory sequences commonly affects gene expression in natural populations and that it provides a seedbed for phenotypic evolution via divergence in gene expression.

Neighborhood characteristics influence DNA methylation of genes involved in stress response and inflammation: The Multi-Ethnic Study of Atherosclerosis

Living in a disadvantaged neighborhood is associated with poor health outcomes even after accounting for individual-level socioeconomic factors. The chronic stress of unfavorable neighborhood conditions may lead to dysregulation of the stress reactivity and inflammatory pathways, potentially mediated through epigenetic mechanisms such as DNA methylation. We used multi-level models to examine the relationship between 2 neighborhood conditions and methylation levels of 18 genes related to stress reactivity and inflammation in purified monocytes from 1,226 participants of the Multi-Ethnic Study of Atherosclerosis (MESA), a population-based sample of US adults. Neighborhood socioeconomic disadvantage, a summary of 16 census-based metrics, was associated with DNA methylation [False discovery rate (FDR) q-value ≤ 0.1] in 2 out of 7 stress-related genes evaluated (CRF, SLC6A4) and 2 out of 11 inflammation-related genes (F8, TLR1). Neighborhood social environment, a summary measure of aesthetic quality, safety, and social cohesion, was associated with methylation in 4 of the 7 stress-related genes (AVP, BDNF, FKBP5, SLC6A4) and 7 of the 11 inflammation-related genes (CCL1, CD1D, F8, KLRG1, NLRP12, SLAMF7, TLR1). High socioeconomic disadvantage and worse social environment were primarily associated with increased methylation. In 5 genes with significant associations between neighborhood and methylation (FKBP5, CD1D, F8, KLRG1, NLRP12), methylation was associated with gene expression of at least one transcript. These results demonstrate that multiple dimensions of neighborhood context may influence methylation levels and subsequent gene expression of stress- and inflammation-related genes, even after accounting for individual socioeconomic factors. Further elucidating the molecular mechanisms underlying these relationships will be important for understanding the etiology of health disparities.

Analysis of Allelic Imbalance in Rice Hybrids Under Water Stress and Association of Asymmetrically Expressed Genes with Drought-Response QTLs

BACKGROUND

Information on the effect of stress on the allele-specific expression (ASE) profile of rice hybrids is limited. More so, the association of allelically imbalanced genes to important traits is yet to be understood. Here we assessed allelic imbalance (AI) in the heterozygote state of rice under non- and water-stress treatments and determined association of asymmetrically expressed genes with grain yield (GY) under drought stress by in-silico co-localization analysis and selective genotyping. The genotypes IR64, Apo and their F1 hybrid (IR64 × Apo) were grown under normal and water-limiting conditions. We sequenced the total RNA transcripts for all genotypes then reconstructed the two chromosomes in the heterozygote.

RESULTS

We are able to estimate the transcript abundance of and the differential expression (DE) between the two parent-specific alleles in the rice hybrids. The magnitude and direction of AI are classified into two categories: (1) symmetrical or biallelic and (2) asymmetrical. The latter can be further classified as either IR64- or Apo-favoring gene. Analysis showed that in the hybrids grown under non-stress conditions, 179 and 183 favor Apo- and IR64-specific alleles, respectively. Hence, the number of IR64- and Apo-favoring genes is relatively equal. Under water-stress conditions, 179 and 255 favor Apo- and IR64-specific alleles, respectively, indicating that the number of allelically imbalanced genes is skewed towards IR64. This is nearly 40-60 % preference for Apo and IR64 alleles, respectively, to the hybrid transcriptome. We also observed genes which exhibit allele preference switching when exposed to water-stress conditions. Results of in-silico co-localization procedure and selective genotyping of Apo/IR64 F3:5 progenies revealed significant association of several asymmetrically expressed genes with GY under drought stress conditions.

CONCLUSION

Our data suggest that water stress skews AI on a genome-wide scale towards the IR64 allele, the cross-specific maternal allele. Several asymmetrically expressed genes are strongly associated with GY under drought stress which may shed hints that genes associated with important traits are allelically imbalanced. Our approach of integrating hybrid expression analysis and QTL mapping analysis may be an efficient strategy for shortlisting candidate genes for gene discovery.

A New Mechanism for Mendelian Dominance in Regulatory Genetic Pathways: Competitive Binding by Transcription Factors

We report a new mechanism for allelic dominance in regulatory genetic interactions that we call binding dominance. We investigated a biophysical model of gene regulation, where the fractional occupancy of a transcription factor (TF) on the cis-regulated promoter site it binds to is determined by binding energy (-ΔG) and TF dosage. Transcription and gene expression proceed when the TF is bound to the promoter. In diploids, individuals may be heterozygous at the cis-site, at the TF's coding region, or at the TF's own promoter, which determines allele-specific dosage. We find that when the TF's coding region is heterozygous, TF alleles compete for occupancy at the cis-sites and the tighter-binding TF is dominant in proportion to the difference in binding strength. When the TF's own promoter is heterozygous, the TF produced at the higher dosage is also dominant. Cis-site heterozygotes have additive expression and therefore codominant phenotypes. Binding dominance propagates to affect the expression of downstream loci and it is sensitive in both magnitude and direction to genetic background, but its detectability often attenuates. While binding dominance is inevitable at the molecular level, it is difficult to detect in the phenotype under some biophysical conditions, more so when TF dosage is high and allele-specific binding affinities are similar. A body of empirical research on the biophysics of TF binding demonstrates the plausibility of this mechanism of dominance, but studies of gene expression under competitive binding in heterozygotes in a diversity of genetic backgrounds are needed.

Gene coexpression network analysis of oil biosynthesis in an interspecific backcross of oil palm

Global demand for vegetable oils is increasing at a dramatic rate, while our understanding of the regulation of oil biosynthesis in plants remains limited. To gain insights into the mechanisms that govern oil synthesis and fatty acid (FA) composition in the oil palm fruit, we used a multilevel approach combining gene coexpression analysis, quantification of allele-specific expression and joint multivariate analysis of transcriptomic and lipid data, in an interspecific backcross population between the African oil palm, Elaeis guineensis, and the American oil palm, Elaeis oleifera, which display contrasting oil contents and FA compositions. The gene coexpression network produced revealed tight transcriptional coordination of fatty acid synthesis (FAS) in the plastid with sugar sensing, plastidial glycolysis, transient starch storage and carbon recapture pathways. It also revealed a concerted regulation, along with FAS, of both the transfer of nascent FA to the endoplasmic reticulum, where triacylglycerol assembly occurs, and of the production of glycerol-3-phosphate, which provides the backbone of triacylglycerols. Plastid biogenesis and auxin transport were the two other biological processes most tightly connected to FAS in the network. In addition to WRINKLED1, a transcription factor (TF) known to activate FAS genes, two novel TFs, termed NF-YB-1 and ZFP-1, were found at the core of the FAS module. The saturated FA content of palm oil appeared to vary above all in relation to the level of transcripts of the gene coding for β-ketoacyl-acyl carrier protein synthase II. Our findings should facilitate the development of breeding and engineering strategies in this and other oil crops.

Variation in the Oxytocin Receptor Gene Predicts Brain Region-Specific Expression and Social Attachment

BACKGROUND

Oxytocin (OXT) modulates several aspects of social behavior. Intranasal OXT is a leading candidate for treating social deficits in patients with autism spectrum disorder, and common genetic variants in the human OXTR gene are associated with emotion recognition, relationship quality, and autism spectrum disorder. Animal models have revealed that individual differences in Oxtr expression in the brain drive social behavior variation. Our understanding of how genetic variation contributes to brain OXTR expression is very limited.

METHODS

We investigated Oxtr expression in monogamous prairie voles, which have a well-characterized OXT system. We quantified brain region-specific levels of Oxtr messenger RNA and oxytocin receptor protein with established neuroanatomic methods. We used pyrosequencing to investigate allelic imbalance of Oxtr mRNA, a molecular signature of polymorphic genetic regulatory elements. We performed next-generation sequencing to discover variants in and near the Oxtr gene. We investigated social attachment using the partner preference test.

RESULTS

Our allelic imbalance data demonstrate that genetic variants contribute to individual differences in Oxtr expression, but only in particular brain regions, including the nucleus accumbens, where oxytocin receptor signaling facilitates social attachment. Next-generation sequencing identified one polymorphism in the Oxtr intron, near a putative cis-regulatory element, explaining 74% of the variance in striatal Oxtr expression specifically. Males homozygous for the high expressing allele display enhanced social attachment.

CONCLUSIONS

Taken together, these findings provide convincing evidence for robust genetic influence on Oxtr expression and provide novel insights into how noncoding polymorphisms in OXTR might influence individual differences in human social cognition and behavior.

Identifying genomic changes associated with insecticide resistance in the dengue mosquito Aedes aegypti by deep targeted sequencing

The capacity of mosquitoes to resist insecticides threatens the control of diseases such as dengue and malaria. Until alternative control tools are implemented, characterizing resistance mechanisms is crucial for managing resistance in natural populations. Insecticide biodegradation by detoxification enzymes is a common resistance mechanism; however, the genomic changes underlying this mechanism have rarely been identified, precluding individual resistance genotyping. In particular, the role of copy number variations (CNVs) and polymorphisms of detoxification enzymes have never been investigated at the genome level, although they can represent robust markers of metabolic resistance. In this context, we combined target enrichment with high-throughput sequencing for conducting the first comprehensive screening of gene amplifications and polymorphisms associated with insecticide resistance in mosquitoes. More than 760 candidate genes were captured and deep sequenced in several populations of the dengue mosquito Ae. aegypti displaying distinct genetic backgrounds and contrasted resistance levels to the insecticide deltamethrin. CNV analysis identified 41 gene amplifications associated with resistance, most affecting cytochrome P450s overtranscribed in resistant populations. Polymorphism analysis detected more than 30,000 variants and strong selection footprints in specific genomic regions. Combining Bayesian and allele frequency filtering approaches identified 55 nonsynonymous variants strongly associated with resistance. Both CNVs and polymorphisms were conserved within regions but differed across continents, confirming that genomic changes underlying metabolic resistance to insecticides are not universal. By identifying novel DNA markers of insecticide resistance, this study opens the way for tracking down metabolic changes developed by mosquitoes to resist insecticides within and among populations.

Genome-wide identification of allele-specific expression in response to Streptococcus suis 2 infection in two differentially susceptible pig breeds

Although allele expression imbalance has been recognized in many species, and strongly linked to diseases, no whole transcriptome allele imbalance has been detected in pigs during pathogen infections. The pathogen Streptococcus suis 2 (SS2) causes serious zoonotic disease. Different pig breeds show differential susceptibility/resistance to pathogen infection, but the biological insight is little known. Here we analyzed allele-specific expression (ASE) using the spleen transcriptome of four pigs belonging to two phenotypically different breeds after SS2 infection. The comparative analysis of allele specific SNPs between control and infected animals revealed 882 and 1096 statistically significant differentially expressed allele SNPs (criteria: ratio ≧ 2 or ≦ 0.5) in Landrace and Enshi black pig, respectively. Twenty nine allelically imbalanced SNPs were further verified by Sanger sequencing, and later six SNPs were quantified by pyrosequencing assay. The pyrosequencing results are in agreement with the RNA-seq results, except two SNPs. Looking at the role of ASE in predisposition to diseases, the discovery of causative variants by ASE analysis might help the pig industry in long term to design breeding programs for improving SS2 resistance.

Comparative analysis of response to selection with three insecticides in the dengue mosquito Aedes aegypti using mRNA sequencing

Background

Mosquito control programmes using chemical insecticides are increasingly threatened by the development of resistance. Such resistance can be the consequence of changes in proteins targeted by insecticides (target site mediated resistance), increased insecticide biodegradation (metabolic resistance), altered transport, sequestration or other mechanisms. As opposed to target site resistance, other mechanisms are far from being fully understood. Indeed, insecticide selection often affects a large number of genes and various biological processes can hypothetically confer resistance. In this context, the aim of the present study was to use RNA sequencing (RNA-seq) for comparing transcription level and polymorphism variations associated with adaptation to chemical insecticides in the mosquito Aedes aegypti. Biological materials consisted of a parental susceptible strain together with three child strains selected across multiple generations with three insecticides from different classes: the pyrethroid permethrin, the neonicotinoid imidacloprid and the carbamate propoxur.

Results

After ten generations, insecticide-selected strains showed elevated resistance levels to the insecticides used for selection. RNA-seq data allowed detecting over 13,000 transcripts, of which 413 were differentially transcribed in insecticide-selected strains as compared to the susceptible strain. Among them, a significant enrichment of transcripts encoding cuticle proteins, transporters and enzymes was observed. Polymorphism analysis revealed over 2500 SNPs showing > 50% allele frequency variations in insecticide-selected strains as compared to the susceptible strain, affecting over 1000 transcripts. Comparing gene transcription and polymorphism patterns revealed marked differences among strains. While imidacloprid selection was linked to the over transcription of many genes, permethrin selection was rather linked to polymorphism variations. Focusing on detoxification enzymes revealed that permethrin selection strongly affected the polymorphism of several transcripts encoding cytochrome P450 monooxygenases likely involved in insecticide biodegradation.

Conclusions

The present study confirmed the power of RNA-seq for identifying concomitantly quantitative and qualitative transcriptome changes associated with insecticide resistance in mosquitoes. Our results suggest that transcriptome modifications can be selected rapidly by insecticides and affect multiple biological functions. Previously neglected by molecular screenings, polymorphism variations of detoxification enzymes may play an important role in the adaptive response of mosquitoes to insecticides.