Differential Allelic Expression Research Articles

Pyrosequencing for accurate imprinted allele expression analysis.

Genomic imprinting is an epigenetic mechanism that restricts gene expression to one inherited allele. Improper maintenance of imprinting has been implicated in a number of human diseases and developmental syndromes. Assays are needed that can quantify the contribution of each paternal allele to a gene expression profile. We have developed a rapid, sensitive quantitative assay for the measurement of individual allelic ratios termed Pyrosequencing for Imprinted Expression (PIE). Advantages of PIE over other approaches include shorter experimental time, decreased labor, avoiding the need for restriction endonuclease enzymes at polymorphic sites, and prevent heteroduplex formation which is problematic in quantitative PCR-based methods. We demonstrate the improved sensitivity of PIE including the ability to detect differences in allelic expression down to 1%. The assay is capable of measuring genomic heterozygosity as well as imprinting in a single run. PIE is applied to determine the status of Insulin-like Growth Factor-2 (IGF2) imprinting in human and mouse tissues.

Recommendations for Accurate Resolution of Gene and Isoform Allele-Specific Expression in RNA-Seq Data.

Genetic variation modulates gene expression transcriptionally or post-transcriptionally, and can profoundly alter an individual’s phenotype. Measuring allelic differential expression at heterozygous loci within an individual, a phenomenon called allele-specific expression (ASE), can assist in identifying such factors. Massively parallel DNA and RNA sequencing and advances in bioinformatic methodologies provide an outstanding opportunity to measure ASE genome-wide. In this study, matched DNA and RNA sequencing, genotyping arrays and computationally phased haplotypes were integrated to comprehensively and conservatively quantify ASE in a single human brain and liver tissue sample. We describe a methodological evaluation and assessment of common bioinformatic steps for ASE quantification, and recommend a robust approach to accurately measure SNP, gene and isoform ASE through the use of personalized haplotype genome alignment, strict alignment quality control and intragenic SNP aggregation. Our results indicate that accurate ASE quantification requires careful bioinformatic analyses and is adversely affected by sample specific alignment confounders and random sampling even at moderate sequence depths. We identified multiple known and several novel ASE genes in liver, including WDR72, DSP and UBD, as well as genes that contained ASE SNPs with imbalance direction discordant with haplotype phase, explainable by annotated transcript structure, suggesting isoform derived ASE. The methods evaluated in this study will be of use to researchers performing highly conservative quantification of ASE, and the genes and isoforms identified as ASE of interest to researchers studying those loci.

Temperature stress mediates decanalization and dominance of gene expression in Drosophila melanogaster.

The regulatory architecture of gene expression remains an area of active research. Here, we studied how the interplay of genetic and environmental variation affects gene expression by exposing Drosophila melanogaster strains to four different developmental temperatures. At 18°C we observed almost complete canalization with only very few allelic effects on gene expression. In contrast, at the two temperature extremes, 13°C and 29°C a large number of allelic differences in gene expression were detected due to both cis- and trans-regulatory effects. Allelic differences in gene expression were mainly dominant, but for up to 62% of the genes the dominance swapped between 13 and 29°C. Our results are consistent with stabilizing selection causing buffering of allelic expression variation in non-stressful environments. We propose that decanalization of gene expression in stressful environments is not only central to adaptation, but may also contribute to genetic disorders in human populations.

Differential allelic expression of SOS1 and hyperexpression of the activating SOS1 c.755C variant in a Noonan syndrome family

Noonan syndrome (NS) is a genetic condition characterized by congenital heart defects, short stature and characteristic facial features. We here present the case of a girl with moderate learning disabilities, delayed language development, craniofacial features and skin anomalies reminiscent of NS. After a mutation screening of the known NS genes PTPN11, SOS1, RAF1, KRAS, GRB2, BRAF and SHOC2 we found the heterozygous c.755T>C variant in SOS1 causing the p.I252T amino-acid substitution, which was considered possibly pathogenetic by bioinformatic predictions. The same variant was present in the proband's mother, displaying some NS features, and maternal grandfather showing no NS traits, but also by a healthy subject in 1000 genomes project database without phenotype informations. The functional analysis revealed that SOS1 c.755C activated the RAS-ERK intracellular pathway, whereas no effects on RAC-JNK cascade have been detected. After a comparison between the sequence of SOS1 cDNA from peripheral blood and SOS1 genomic DNA, we showed for the first time a differential allelic expression of the SOS1 gene in healthy individuals, thus occurring as a physiologic condition. Interestingly, we found that the mutated allele C was 50% more expressed than the wild-type allele T in all familial carriers. The comparable amount of SOS1 mRNA between mutated individuals and the controls indicates that the variant does not affect SOS1 expression. The present study provides a first evidence of allelic imbalance of SOS1 and pinpoints this condition as a possible mechanism underlying a different penetrance of some SOS1-mutated alleles in unrelated carriers.

Control of Paternally Expressed Imprinted UPWARD CURLY LEAF1, a Gene Encoding an F-Box Protein That Regulates CURLY LEAF Polycomb Protein, in the Arabidopsis Endosperm.

Genomic imprinting, an epigenetic process in mammals and flowering plants, refers to the differential expression of alleles of the same genes in a parent-of-origin-specific manner. In Arabidopsis, imprinting occurs primarily in the endosperm, which nourishes the developing embryo. Recent high-throughput sequencing analyses revealed that more than 200 loci are imprinted in Arabidopsis; however, only a few of these imprinted genes and their imprinting mechanisms have been examined in detail. Whereas most imprinted loci characterized to date are maternally expressed imprinted genes (MEGs), PHERES1 (PHE1) and ADMETOS (ADM) are paternally expressed imprinted genes (PEGs). Here, we report that UPWARD CURLY LEAF1 (UCL1), a gene encoding an E3 ligase that degrades the CURLY LEAF (CLF) polycomb protein, is a PEG. After fertilization, paternally inherited UCL1 is expressed in the endosperm, but not in the embryo. The expression pattern of a β-glucuronidase (GUS) reporter gene driven by the UCL1 promoter suggests that the imprinting control region (ICR) of UCL1 is adjacent to a transposable element in the UCL1 5′-upstream region. Polycomb Repressive Complex 2 (PRC2) silences the maternal UCL1 allele in the central cell prior to fertilization and in the endosperm after fertilization. The UCL1 imprinting pattern was not affected in paternal PRC2 mutants. We found unexpectedly that the maternal UCL1 allele is reactivated in the endosperm of Arabidopsis lines with mutations in cytosine DNA METHYLTRANSFERASE 1 (MET1) or the DNA glycosylase DEMETER (DME), which antagonistically regulate CpG methylation of DNA. By contrast, maternal UCL1 silencing was not altered in mutants with defects in non-CpG methylation. Thus, silencing of the maternal UCL1 allele is regulated by both MET1 and DME as well as by PRC2, suggesting that divergent mechanisms for the regulation of PEGs evolved in Arabidopsis.

The clinical presentation of Marfan syndrome is modulated by expression of wild-type FBN1 allele.

Marfan syndrome is an autosomal dominant disorder mainly caused by mutations within FBN1 gene. The disease displays large variability in age of onset or severity and very poor phenotype/genotype correlations have been demonstrated. We investigated the hypothesis that phenotype severity could be related to the variable expression level of fibrillin-1 (FBN1) synthesized from the wild-type (WT) allele. Quantitative reverse-transcription and polymerase chain reaction was used to evaluate FBN1 levels in skin fibroblasts from 80 Marfan patients with premature termination codons and in skin fibroblasts from 80 controls. Results in controls showed a 3.9-fold variation in FBN1 mRNA synthesis level between subjects. A similar 4.4-fold variation was found in the Marfan population, but the mean level of FBN1 mRNA was a half of the control population. Differential allelic expression analysis in Marfan fibroblasts showed that over 90% of FBN1 mRNA was transcribed from the wild allele and the mutated allele was not detected. In the control population, independently of the expression level of FBN1, we observed steady-state equilibrium between the two allelic-mRNAs suggesting that FBN1 expression mainly depends on trans-acting regulators. Finally, we show that a low level of residual WT FBN1 mRNA accounts for a high risk of ectopia lentis and pectus abnormality and tends to increase the risk of aortic dilatation.

Negative feedback buffers effects of regulatory variants

Mechanisms conferring robustness against regulatory variants have been controversial. Previous studies suggested widespread buffering of RNA misexpression on protein levels during translation. We do not find evidence that translational buffering is common. Instead, we find extensive buffering at the level of RNA expression, exerted through negative feedback regulation acting in trans, which reduces the effect of regulatory variants on gene expression. Our approach is based on a novel experimental design in which allelic differential expression in a yeast hybrid strain is compared to allelic differential expression in a pool of its spores. Allelic differential expression in the hybrid is due to cis-regulatory differences only. Instead, in the pool of spores allelic differential expression is not only due to cis-regulatory differences but also due to local trans effects that include negative feedback. We found that buffering through such local trans regulation is widespread, typically compensating for about 15% of cis-regulatory effects on individual genes. Negative feedback is stronger not only for essential genes, indicating its functional relevance, but also for genes with low to middle levels of expression, for which tight regulation matters most. We suggest that negative feedback is one mechanism of Waddington's canalization, facilitating the accumulation of genetic variants that might give selective advantage in different environments.

Detecting Allelic Expression Imbalance at Candidate Genes Using 5' Exonuclease Genotyping Technology.

Genetic variation along the length of a chromosome can influence the transcription of a gene. In a heterozygous individual, this may lead to one chromosome producing different levels of RNA, compared to its paired chromosome, for a given gene. Allelic differences in gene expression can offer insight into the role of variation in transcription, and subsequently infer a route to conferring disease risk. This phenomenon is known as allele expression imbalance or AEI, which may be assayed using a PCR-based method that includes the quantification of the relative dosage of each allele (e.g., 5' exonuclease assays, TaqMan™). Importantly, in heterozygous individuals the resolution of expression imbalance is performed within a controlled system; the comparison of the alternate allele is reported relative to the wild-type, as the experiment can be performed within a single sample, controlled for background genetic information. Alternative methods for the detection of AEI include Primer-extension MALDI-TOF (Sequenom MassARRAY(®)), Next-Generation Sequencing, and SNP genotyping arrays. Here we present the methods used for the TaqMan™ approach and include a description of the SNP identification, allele-specific PCR, and analytic methods to convert allele amplification metrics to relative allele dosage.

Frequent Igh Fusion Transcripts with Clinical Impact in Multiple Myeloma

Background: Significant heterogeneity has been described in Multiple Myloma (MM), especially at the genomic level. Frequent gains and losses of DNA along with various mutations have been observed, and differential allelic expression is being characterized. Fusion proteins are common and maybe associated with cell transformation, growth and lethality of tumor cells. However, unlike in leukemia, no consistent fusion gene product has been consistently identified. As IgH-related translocations have an important role in myeloma, we have investigated fusion genes involving IgH, to understand their biology and explore a possible effect on survival in MM.Methods: We performed deep RNA-Seq on purified MM cells from 430 newly-diagnosed MM patients and analyzed gene expression profiles, isoform signatures and both novel and known fusion genes using two common algorithms: TopHat and MapSplice. We also correlated genomic data with patient data including cytogenetics and FISH, as well as survival.Results: We primarily focused on fusions involving the IGH gene (chromosome 14) and found that about one fourth of the patients (57 out of 430) presented an IGH fusion gene (97 patients according to TopHat, 303 according to Mapsplice, 57 according to both). These included the well described t(4,14) fusion involving the MMSET gene (found in 47 patients by both algorithms, 49 by Tophat, 54 by Mapsplice). Additionally we observed fusions involving chromosomes 14 and chromosomes 1, 4, 11, 12, and 16. The counterpart genes involved in the IgH fusions included PDE3A (chromosome 12 - 4 patients); HFM1 (chromosome 1, 2 patients); NFKB1, FGFR3, CIITA, WWOX and MRPL21 (chromosomes 4, 16, and 11, 1 patient). As RNA-Seq data allows the precise localization of the breakpoints, we were able to identify that out of the 47 t(4,14) patients, 62% were MB4-1, 9.5% were MB4-2 and 28.5% were MB4-3. Interestingly, we did not see fusion products involving IgH and other known parts on Chromosomes 8 and 20. We studied event-free survival in a subset of 265 patients with available survival data and found that, as predicted, patients with an IgH-MMSET fusion had significantly lower survival than others. However, patients with a fusion gene involving IGH and any other partner have a significantly better prognosis as a group. Moreover, the poor prognosis of IgH-MMSET fusion appears to be driven by MB4-3 patients. Importantly, the fusions identified using RNA-seq were also validated by FISH analysis. All t(4,14) fusions were characterized by a very high MMSET expression (FPKM greater than 20) while patients with other fusions presented a lower MMSET expression (FPKM lower than 10).Conclusion: Our study suggests that IgH-related translocations in myeloma may impact tumor biology by a number of mechanism, one of which is the generation of fusion proteins with functional consequences. It also highlights a possible clinical impact that requires validation in larger cohorts. DisclosuresNo relevant conflicts of interest to declare.

Differential and Limited Expression of Mutant Alleles in Multiple Myeloma

Multiple Myeloma (MM) is a heterogeneous disease but the hallmark genetic changes involve large numbers of genomic rearrangements. Recent studies have focused on attempts to identify individual driver mutations that might provide both prognostic information and unique therapeutic targets. Whole genome and exome sequencing of increasingly large numbers of patient samples have identified a number of commonly mutated genes in MM patients. However, none of these mutations are found in more than one quarter of patients and most are found in less than 10% of samples sequenced. We recently reported a large cohort of MM exome sequences involving 84 samples from 67 patients (Nat Commun. 2014;5:2997). We defined a diverse set of gene mutations with significant heterogeneity across our cohort with a median of 52 (range 21-488) mutations identified per sample. Although computational approaches can be used to prioritize mutations that are expected to alter protein structure and function, it is more challenging to determine which mutations are likely to be clinically meaningful. As a first step towards that understanding, here we report the frequency of expression of mutant alleles in Multiple Myeloma.In this study we report RNA-seq (100 million paired end reads on Illumina HiSeq) data on 14 samples from 10 MM patients for which we have previously performed exome sequencing and correlate allele-specific expression to the DNA mutant allele frequency. We find that a minority, average 27% (range 11-48%), of previously identified DNA mutations are expressed at detectable levels in MM patients. We also compared the allele frequency found in the RNA-seq to that from our exome sequencing to identify genes that demonstrate differential allelic expression and show that this is a common phenomenon in MM patients. We identified 42 such mutations in our analysis supported by at least 10 RNA-seq reads that showed a significant difference as determined by Bayesian hypothesis testing. For instance, the CCND1 mutant allele is expressed at a higher level than would be predicted based on exome-seq frequencies. Another gene showing a similar pattern of increased expression of the mutant allele in one patient was PARP4 (87% in RNA-seq vs 49% in exome-seq). Conversely, the mutant allele frequency of EIF1AX was lower than would be expected suggesting that the mutant allele may be suppressed in our patient (15% in RNA-seq vs 67% in exome-seq). Moreover, among a subset of genes previously identified as recurrently mutated within our patient samples we see that 8/11 (73%) express the mutant allele, providing further evidence that these genes may in fact be important in disease pathogenesis. Therefore, while a large number of mutations have been described in MM, only a small fraction of the mutant alleles have detectable expression and are likely to be biologically relevant. Unbalanced allelic expression of mutant alleles appears to be a relatively common occurrence in MM patients and may help explain why patients with the same identified mutation do not always behave in a similar fashion. This analysis for the first time highlights the important issue that DNA-based reporting of mutations may have significant limitations. It will be important in the future to study expression of mutant alleles in order to understand the biology, generate prognostic models and develop targeted therapies. DisclosuresNo relevant conflicts of interest to declare.

Differential and limited expression of mutant alleles in multiple myeloma

AbstractRecent work has delineated mutational profiles in multiple myeloma and reported a median of 52 mutations per patient, as well as a set of commonly mutated genes across multiple patients. In this study, we have used deep sequencing of RNA from a subset of these patients to evaluate the proportion of expressed mutations. We find that the majority of previously identified mutations occur within genes with very low or no detectable expression. On average, 27% (range, 11% to 47%) of mutated alleles are found to be expressed, and among mutated genes that are expressed, there often is allele-specific expression where either the mutant or wild-type allele is suppressed. Even in the absence of an overall change in gene expression, the presence of differential allelic expression within malignant cells highlights the important contribution of RNA-sequencing in identifying clinically significant mutational changes relevant to our understanding of myeloma biology and also for therapeutic applications.

Ecological genomics in full colour.

Colour patterns in animals have long offered an opportunity to observe adaptive traits in natural populations. Colour plays myriad roles in interactions within and among species, from reproductive signalling to predator avoidance, leading to multiple targets of natural and sexual selection and opportunities for diversification. Understanding the genetic and developmental underpinnings of variation in colour promises a fuller understanding of these evolutionary processes, but the path to unravelling these connections can be arduous. The advent of genomic techniques suitable for nonmodel organisms is now beginning to light the way. Two new studies in this issue of Molecular Ecology use genomic sequencing of laboratory crosses to map colour traits in cichlid fishes, a remarkably diverse group in which coloration has played a major role in diversification. They illustrate how genomic approaches, specifically RAD sequencing, can rapidly identify both simple and more complex genetic variation underlying ecologically important traits. In the first, Henning etal. () detect a single locus that appears to control in a Mendelian fashion the presence of horizontal stripes, a trait that has evolved in numerous cichlid lineages. In the second, Albertson etal. () identify several genes and epistatic interactions affecting multiple colour traits, as well as a novel metric describing integration across colour traits. Albertson etal. () go further, by quantifying differential expression of parental alleles at a candidate locus and by relating differentiation among natural populations at mapped loci to trait divergence. Herein lies the promise of ecological genomics - efficiently integrating genetic mapping of phenotypes with population genomic data to both identify functional genes and unravel their evolutionary history. These studies offer guidance on how genomic techniques can be tailored to a research question or study system, and they also add to the growing body of empirical examples addressing basic questions about how ecologically important traits evolve in natural populations.

Regulatory polymorphisms in human DBH affect peripheral gene expression and sympathetic activity.

Dopamine β-hydroxylase (DBH) catalyzes the conversion of dopamine to norepinephrine in the central nervous system and peripherally. DBH variants are associated with large changes in circulating DBH and implicated in multiple disorders; yet causal relationships and tissue-specific effects remain unresolved. To characterize regulatory variants in DBH, effect on mRNA expression, and role in modulating sympathetic tone and disease risk. Analysis of DBH mRNA in human tissues confirmed high expression in the locus coeruleus and adrenal gland, but also in sympathetically innervated organs (liver>lung>heart). Allele-specific mRNA assays revealed pronounced allelic expression differences in the liver (2- to 11-fold) attributable to promoter rs1611115 and exon 2 rs1108580, but only small differences in locus coeruleus and adrenals. These alleles were also associated with significantly reduced mRNA expression in liver and lung. Although DBH protein is expressed in other sympathetically innervated organs, mRNA levels were too low for analysis. In mice, hepatic Dbh mRNA levels correlated with cardiovascular risk phenotypes. The minor alleles of rs1611115 and rs1108580 were associated with sympathetic phenotypes, including angina pectoris. Testing combined effects of these variants suggested protection against myocardial infarction in 3 separate clinical cohorts. We demonstrate profound effects of DBH variants on expression in 2 sympathetically innervated organs, liver and lung, but not in adrenals and brain. Preliminary results demonstrate an association of these variants with clinical phenotypes responsive to peripheral sympathetic tone. We hypothesize that in addition to endocrine effects via circulating DBH and norepinephrine, the variants act in sympathetically innervated target organs.

Allelic expression mapping across cellular lineages to establish impact of non-coding SNPs.

Most complex disease-associated genetic variants are located in non-coding regions and are therefore thought to be regulatory in nature. Association mapping of differential allelic expression (AE) is a powerful method to identify SNPs with direct cis-regulatory impact (cis-rSNPs). We used AE mapping to identify cis-rSNPs regulating gene expression in 55 and 63 HapMap lymphoblastoid cell lines from a Caucasian and an African population, respectively, 70 fibroblast cell lines, and 188 purified monocyte samples and found 40–60% of these cis-rSNPs to be shared across cell types. We uncover a new class of cis-rSNPs, which disrupt footprint-derived de novo motifs that are predominantly bound by repressive factors and are implicated in disease susceptibility through overlaps with GWAS SNPs. Finally, we provide the proof-of-principle for a new approach for genome-wide functional validation of transcription factor–SNP interactions. By perturbing NFκB action in lymphoblasts, we identified 489 cis-regulated transcripts with altered AE after NFκB perturbation. Altogether, we perform a comprehensive analysis of cis-variation in four cell populations and provide new tools for the identification of functional variants associated to complex diseases.

Cellular adhesion gene SELP is associated with rheumatoid arthritis and displays differential allelic expression.

In rheumatoid arthritis (RA), a key event is infiltration of inflammatory immune cells into the synovial lining, possibly aggravated by dysregulation of cellular adhesion molecules. Therefore, single nucleotide polymorphisms of 14 genes involved in cellular adhesion processes (CAST, ITGA4, ITGB1, ITGB2, PECAM1, PTEN, PTPN11, PTPRC, PXN, SELE, SELP, SRC, TYK2, and VCAM1) were analyzed for association with RA. Association analysis was performed consecutively in three European RA family sample groups (Nfamilies = 407). Additionally, we investigated differential allelic expression, a possible functional consequence of genetic variants. SELP (selectin P, CD62P) SNP-allele rs6136-T was associated with risk for RA in two RA family sample groups as well as in global analysis of all three groups (ptotal = 0.003). This allele was also expressed preferentially (p<10−6) with a two- fold average increase in regulated samples. Differential expression is supported by data from Genevar MuTHER (p1 = 0.004; p2 = 0.0177). Evidence for influence of rs6136 on transcription factor binding was also found in silico and in public datasets reporting in vitro data. In summary, we found SELP rs6136-T to be associated with RA and with increased expression of SELP mRNA. SELP is located on the surface of endothelial cells and crucial for recruitment, adhesion, and migration of inflammatory cells into the joint. Genetically determined increased SELP expression levels might thus be a novel additional risk factor for RA.

RNA-Seq using two populations reveals genes and alleles controlling wood traits and growth in Eucalyptus nitens.

Eucalyptus nitens is a perennial forest tree species grown mainly for kraft pulp production in many parts of the world. Kraft pulp yield (KPY) is a key determinant of plantation profitability and increasing the KPY of trees grown in plantations is a major breeding objective. To speed up the breeding process, molecular markers that can predict KPY are desirable. To achieve this goal, we carried out RNA-Seq studies on trees at extremes of KPY in two different trials to identify genes and alleles whose expression correlated with KPY. KPY is positively correlated with growth measured as diameter at breast height (DBH) in both trials. In total, six RNA bulks from two treatments were sequenced on an Illumina HiSeq platform. At 5% false discovery rate level, 3953 transcripts showed differential expression in the same direction in both trials; 2551 (65%) were down-regulated and 1402 (35%) were up-regulated in low KPY samples. The genes up-regulated in low KPY trees were largely involved in biotic and abiotic stress response reflecting the low growth among low KPY trees. Genes down-regulated in low KPY trees mainly belonged to gene categories involved in wood formation and growth. Differential allelic expression was observed in 2103 SNPs (in 1068 genes) and of these 640 SNPs (30%) occurred in 313 unique genes that were also differentially expressed. These SNPs may represent the cis-acting regulatory variants that influence total gene expression. In addition we also identified 196 genes which had Ka/Ks ratios greater than 1.5, suggesting that these genes are under positive selection. Candidate genes and alleles identified in this study will provide a valuable resource for future association studies aimed at identifying molecular markers for KPY and growth.

Mouse pulmonary adenoma susceptibility 1 locus is an expression QTL modulating Kras-4A.

Pulmonary adenoma susceptibility 1 (Pas1) is the major locus responsible for lung tumor susceptibility in mice; among the six genes mapping in this locus, Kras is considered the best candidate for Pas1 function although how it determines tumor susceptibility remains unknown. In an (A/J×C57BL/6)F4 intercross population treated with urethane to induce lung tumors, Pas1 not only modulated tumor susceptibility (LOD score = 48, 69% of phenotypic variance explained) but also acted, in lung tumor tissue, as an expression quantitative trait locus (QTL) for Kras-4A, one of two alternatively spliced Kras transcripts, but not Kras-4B. Additionally, Kras-4A showed differential allelic expression in lung tumor tissue of (A/J×C57BL/6)F4 heterozygous mice, with significantly higher expression from the A/J-derived allele; these results suggest that cis-acting elements control Kras-4A expression. In normal lung tissue from untreated mice of the same cross, Kras-4A levels were also highly linked to the Pas1 locus (LOD score = 23.2, 62% of phenotypic variance explained) and preferentially generated from the A/J-derived allele, indicating that Pas1 is an expression QTL in normal lung tissue as well. Overall, the present findings shed new light on the genetic mechanism by which Pas1 modulates the susceptibility to lung tumorigenesis, through the fine control of Kras isoform levels.

Using next-generation RNA sequencing to identify imprinted genes.

Genomic imprinting is manifested as differential allelic expression (DAE) depending on the parent-of-origin. The most direct way to identify imprinted genes is to directly score the DAE in a context where one can identify which parent transmitted each allele. Because many genes display DAE, simply scoring DAE in an individual is not sufficient to identify imprinted genes. In this paper, we outline many technical aspects of a scheme for identification of imprinted genes that makes use of RNA sequencing (RNA-seq) from tissues isolated from F1 offspring derived from the pair of reciprocal crosses. Ideally, the parental lines are from two inbred strains that are not closely related to each other. Aspects of tissue purity, RNA extraction, library preparation and bioinformatic inference of imprinting are all covered. These methods have already been applied in a number of organisms, and one of the most striking results is the evolutionary fluidity with which novel imprinted genes are gained and lost within genomes. The general methodology is also applicable to a wide range of other biological problems that require quantification of allele-specific expression using RNA-seq, such as cis-regulation of gene expression, X chromosome inactivation and random monoallelic expression.

Seed evolution: parental conflicts in a multi-generational household.

Seeds are multi-generational structures containing a small embryonic plant enclosed in layers of diverse parental origins. The evolution of seeds was a pinnacle in an evolutionary trend towards a progressive retention of embryos and gametes within parental tissue. This strategy, which dates back to the first land plants, allowed an increased protection and nourishing of the developing embryo. Flowering plants took parental control one step further with the evolution of a biparental endosperm that derives from a second parallel fertilization event. The endosperm directly nourishes the developing embryo and allows not only the maternal genes, but also paternal genes, to play an active role during seed development. The appearance of an endosperm set the conditions for the manifestation of conflicts of interest between maternal and paternal genomes over the allocation of resources to the developing embryos. As a consequence, a dynamic balance was established between maternal and paternal gene dosage in the endosperm, and maintaining a correct balance became essential to ensure a correct seed development. This balance was achieved in part by changes in the genetic constitution of the endosperm and through epigenetic mechanisms that allow a differential expression of alleles depending on their parental origin. This review discusses the evolutionary steps that resulted in the appearance of seeds and endosperm, and the epigenetic and genetic mechanisms that allow a harmonious coinhabitance of multiple generations within a single seed.

Different timing and spatial separation of parental chromosomes in intergeneric somatic hybrids between Brassica napus and Orychophragmus violaceus

Experimental and newly formed hybrids and polyploids generated by wide crosses usually show varying degrees of cytological instability. The spatial separation of parental genomes and uniparental chromosome elimination in hybrid cells has been reported in many hybrids from plants and animals. Herein, the behavior of parental genomes in intergeneric somatic hybrids between Brassica napus and Orychophragmus violaceus was analyzed using genomic in situ hybridization (GISH). In mitotic and meiotic cells, the chromosomes from O. violaceus were distinguished from B. napus by their larger size and staining patterns. In interphase nuclei of the hybrid, O. violaceus-labeled chromatin appeared as large heterochromatic blocks that were nonrandomly distributed at prophase, typically distributed toward one side of the nucleus. In pollen mother cells at prophase I of meiosis, O. violaceus chromosomes appeared as one or two deeply stained chromatin blocks that resolved into bivalents at a late stage, after bivalents from B. napus were visible. Thereafter, bivalents of O. violaceus congressed to the equatorial plate and segregated at anaphase I after those from B. napus. The different behavior of O. violaceus chromosomes in the hybrids indicates that they have differential condensation states at interphase and progress later through the cell cycle and meiosis than B. napus chromosomes. This difference in behavior may restrict or prevent the formation of bivalents of mixed genome origin. Differential gene expression of parental alleles including rDNA loci may contribute to their distinct cytological behavior and to the phenotype of hybrids.