Hongmei Jiang

Genomewide Nonadditive Gene Regulation in Arabidopsis Allotetraploids

Polyploidy has occurred throughout the evolutionary history of all eukaryotes and is extremely common in plants. Reunification of the evolutionarily divergent genomes in allopolyploids creates regulatory incompatibilities that must be reconciled. Here we report genomewide gene expression analysis of Arabidopsis synthetic allotetraploids, using spotted 70-mer oligo-gene microarrays. We detected >15% transcriptome divergence between the progenitors, and 2105 and 1818 genes were highly expressed in Arabidopsis thaliana and A. arenosa, respectively. Approximately 5.2% (1362) and 5.6% (1469) genes displayed expression divergence from the midparent value (MPV) in two independently derived synthetic allotetraploids, suggesting nonadditive gene regulation following interspecific hybridization. Remarkably, the majority of nonadditively expressed genes in the allotetraploids also display expression changes between the parents, indicating that transcriptome divergence is reconciled during allopolyploid formation. Moreover, >65% of the nonadditively expressed genes in the allotetraploids are repressed, and >94% of the repressed genes in the allotetraploids match the genes that are expressed at higher levels in A. thaliana than in A. arenosa, consistent with the silencing of A. thaliana rRNA genes subjected to nucleolar dominance and with overall suppression of the A. thaliana phenotype in the synthetic allotetraploids and natural A. suecica. The nonadditive gene regulation is involved in various biological pathways, and the changes in gene expression are developmentally regulated. In contrast to the small effects of genome doubling on gene regulation in autotetraploids, the combination of two divergent genomes in allotetraploids by interspecific hybridization induces genomewide nonadditive gene regulation, providing a molecular basis for de novo variation and allopolyploid evolution.

Epigenetic natural variation in Arabidopsis thaliana

Cytosine methylation of repetitive sequences is widespread in plant genomes, occurring in both symmetric (CpG and CpNpG) as well as asymmetric sequence contexts. We used the methylation-dependent restriction enzyme McrBC to profile methylated DNA using tiling microarrays of Arabidopsis Chromosome 4 in two distinct ecotypes, Columbia and Landsberg erecta. We also used comparative genome hybridization to profile copy number polymorphisms. Repeated sequences and transposable elements (TEs), especially long terminal repeat retrotransposons, are densely methylated, but one third of genes also have low but detectable methylation in their transcribed regions. While TEs are almost always methylated, genic methylation is highly polymorphic, with half of all methylated genes being methylated in only one of the two ecotypes. A survey of loci in 96 Arabidopsis accessions revealed a similar degree of methylation polymorphism. Within-gene methylation is heritable, but is lost at a high frequency in segregating F 2 families. Promoter methylation is rare, and gene expression is not generally affected by differences in DNA methylation. Small interfering RNA are preferentially associated with methylated TEs, but not with methylated genes, indicating that most genic methylation is not guided by small interfering RNA. This may account for the instability of gene methylation, if occasional failure of maintenance methylation cannot be restored by other means.

Epigenomic consequences of immortalized plant cell suspension culture.

Plant cells grown in culture exhibit genetic and epigenetic instability. Using a combination of chromatin immunoprecipitation and DNA methylation profiling on tiling microarrays, we have mapped the location and abundance of histone and DNA modifications in a continuously proliferating, dedifferentiated cell suspension culture of Arabidopsis. We have found that euchromatin becomes hypermethylated in culture and that a small percentage of the hypermethylated genes become associated with heterochromatic marks. In contrast, the heterochromatin undergoes dramatic and very precise DNA hypomethylation with transcriptional activation of specific transposable elements (TEs) in culture. High throughput sequencing of small interfering RNA (siRNA) revealed that TEs activated in culture have increased levels of 21-nucleotide (nt) siRNA, sometimes at the expense of the 24-nt siRNA class. In contrast, TEs that remain silent, which match the predominant 24-nt siRNA class, do not change significantly in their siRNA profiles. These results implicate RNA interference and chromatin modification in epigenetic restructuring of the genome following the activation of TEs in immortalized cell culture.

Gene expression analysis at the intersection of ploidy and hybridity in maize

Heterosis and polyploidy are two important aspects of plant evolution. To examine these issues, we conducted a global gene expression study of a maize ploidy series as well as a set of tetraploid inbred and hybrid lines. This gene expression analysis complements an earlier phenotypic study of these same materials. We find that ploidy change affects a large fraction of the genome, albeit at low levels; gene expression changes rarely exceed 2-fold and are typically not statistically significant. The most common gene expression profile we detected is greater than linear increase from monoploid to diploid, and reductions from diploid to triploid and from triploid to tetraploid, a trend that mirrors plant stature. When examining heterosis in tetraploid maize lines, we found a large fraction of the genome impacted but the majority of changes were not statistically significant at 2-fold or less. Non-additive expression was common in the hybrids, and the extent of non-additivity increased both in number and magnitude from duplex to quadruplex hybrids. Overall, we find that gene expression trends mirror observations from the phenotypic studies; however, obvious mechanistic connections remain unknown.

Review of current methods, applications, and data management for the bioinformatics analysis of whole exome sequencing

The advent of next-generation sequencing technologies has greatly promoted advances in the study of human diseases at the genomic, transcriptomic, and epigenetic levels. Exome sequencing, where the coding region of the genome is captured and sequenced at a deep level, has proven to be a cost-effective method to detect disease-causing variants and discover gene targets. In this review, we outline the general framework of whole exome sequence data analysis. We focus on established bioinformatics tools and applications that support five analytical steps: raw data quality assessment, preprocessing, alignment, post-processing, and variant analysis (detection, annotation, and prioritization). We evaluate the performance of open-source alignment programs and variant calling tools using simulated and benchmark datasets, and highlight the challenges posed by the lack of concordance among variant detection tools. Based on these results, we recommend adopting multiple tools and resources to reduce false positives and increase the sensitivity of variant calling. In addition, we briefly discuss the current status and solutions for big data management, analysis, and summarization in the field of bioinformatics.

Multiple SET Methyltransferases Are Required to Maintain Normal Heterochromatin Domains in the Genome of Drosophila melanogaster

Methylation of histone H3 lysine 9 (H3K9) is a key feature of silent chromatin and plays an important role in stabilizing the interaction of heterochromatin protein 1 (HP1) with chromatin. Genomes of metazoans such as the fruit fly Drosophila melanogaster generally encode three types of H3K9-specific SET domain methyltransferases that contribute to chromatin homeostasis during the life cycle of the organism. SU(VAR)3-9, dG9a, and dSETDB1 all function in the generation of wild-type H3K9 methylation levels in the Drosophila genome. Two of these enzymes, dSETDB1 and SU(VAR)3-9, govern heterochromatin formation in distinct but overlapping patterns across the genome. H3K9 methylation in the small, heterochromatic fourth chromosome of D. melanogaster is governed mainly by dSETDB1, whereas dSETDB1 and SU(VAR)3-9 function in concert to methylate H3K9 in the pericentric heterochromatin of all chromosomes, with dG9a having little impact in these domains, as shown by monitoring position effect variegation. To understand how these distinct heterochromatin compartments may be differentiated, we examined the developmental timing of dSETDB1 function using a knockdown strategy. dSETDB1 acts to maintain heterochromatin during metamorphosis, at a later stage in development than the reported action of SU(VAR)3-9. Surprisingly, depletion of both of these enzymes has less deleterious effect than depletion of one. These results imply that dSETDB1 acts as a heterochromatin maintenance factor that may be required for the persistence of earlier developmental events normally governed by SU(VAR)3-9. In addition, the genetic interactions between dSETDB1 and Su(var)3-9 mutations emphasize the importance of maintaining the activities of these histone methyltransferases in balance for normal genome function.

miR-9a Minimizes the Phenotypic Impact of Genomic Diversity by Buffering a Transcription Factor

Gene expression has to withstand stochastic, environmental, and genomic perturbations. For example, in the latter case, 0.5%-1% of the human genome is typically variable between any two unrelated individuals. Such diversity might create problematic variability in the activity of gene regulatory networks and, ultimately, in cell behaviors. Using multigenerational selection experiments, we find that for the Drosophila proneural network, the effect of genomic diversity is dampened by miR-9a-mediated regulation of senseless expression. Reducing miR-9a regulation of the Senseless transcription factor frees the genomic landscape to exert greater phenotypic influence. Whole-genome sequencing identified genomic loci that potentially exert such effects. A larger set of sequence variants, including variants within proneural network genes, exhibits these characteristics when miR-9a concentration is reduced. These findings reveal that microRNA-target interactions may be a key mechanism by which the impact of genomic diversity on cell behavior is dampened.Gene expression has to withstand stochastic, environmental, and genomic perturbations. For example, in the latter case, 0.5%-1% of the human genome is typically variable between any two unrelated individuals. Such diversity might create problematic variability in the activity of gene regulatory networks and, ultimately, in cell behaviors. Using multigenerational selection experiments, we find that for the Drosophila proneural network, the effect of genomic diversity is dampened by miR-9a-mediated regulation of senseless expression. Reducing miR-9a regulation of the Senseless transcription factor frees the genomic landscape to exert greater phenotypic influence. Whole-genome sequencing identified genomic loci that potentially exert such effects. A larger set of sequence variants, including variants within proneural network genes, exhibits these characteristics when miR-9a concentration is reduced. These findings reveal that microRNA-target interactions may be a key mechanism by which the impact of genomic diversity on cell behavior is dampened.

Genome-wide study of DNA methylation alterations in response to diazinon exposure in vitro

Pesticide exposure has repeatedly been associated with cancers. However, molecular mechanisms are largely undetermined. In this study, we examined whether exposure to diazinon, a common organophosphate that has been associated with cancers, could induce DNA methylation alterations. We conducted genome-wide DNA methylation analyses on DNA samples obtained from human hematopoietic K562 cell exposed to diazinon and ethanol using the Illumina Infinium HumanMethylation27 BeadChip. Bayesian-adjusted t-tests were used to identify differentially methylated gene promoter CpG sites. We identified 1069 CpG sites in 984 genes with significant methylation changes in diazinon-treated cells. Gene ontology analysis demonstrated that some genes are tumor suppressor genes, such as TP53INP1 (3.0-fold, q-value <0.001) and PTEN (2.6-fold, q-value <0.001), some genes are in cancer-related pathways, such as HDAC3 (2.2-fold, q-value=0.002), and some remain functionally unknown. Our results provided direct experimental evidence that diazinon may modify gene promoter DNA methylation levels, which may play a pathological role in cancer development.

Changes in gene expression within the nucleus accumbens and striatum following sexual experience

Sexual experience, like repeated drug use, produces long‐term changes including sensitization in the nucleus accumbens and dorsal striatum. To better understand the molecular mechanisms underlying the neuroadaptations following sexual experience, we employed a DNA microarray approach to identify genes differentially expressed between sexually experienced and sexually naïve female hamsters within the nucleus accumbens and dorsal striatum. For 6 weeks, a stimulus male was placed in the home cage of one‐half of the hormonally primed, ovariectomized female hamsters. On the seventh week, the two experimental groups were subdivided, with one half paired with a stimulus male. In comparison with sexually naïve animals, sexually experienced hamsters receiving a stimulus male on week 7 exhibited an increase in a large number of genes. Conversely, sexually experienced female hamsters not receiving a stimulus male on week 7 exhibited a reduction in the expression of many genes. For directional changes and the categories of genes regulated by the experimental conditions, data were consistent across the nucleus accumbens and dorsal striatum. However, the specific genes exhibiting changes in expression were disparate. These experiments, among the first to profile genes regulated by female sexual behavior, will provide insight into the mechanisms by which both motivated behaviors and drugs of abuse induce long‐term changes in the mesolimbic and nigrostriatal dopamine pathways.

Cost-effectiveness of population-based BRCA1/2 testing and ovarian cancer prevention for Ashkenazi Jews: A call for dialogue

Purpose: About half of unaffected BRCA1/2 carriers have a negative family history, confounding efforts toward presymptomatic carrier identification. Ovarian cancer is preventable for known carriers but is otherwise highly lethal. Cost-effectiveness and gains in life expectancy are important factors in evaluating the desirability of population-based genetic screening, currently the only viable strategy to identify carriers with unrevealing family histories.Methods: Cost-utility analysis for a population-based genetic screening program offered to American Ashkenazi Jewish women aged 35–55 years measuring cancer incidence, life expectancy, and cost.Results: Our model predicts that a genetic screening program would result in 2811 fewer cases of ovarian cancer, with a life expectancy gain of 1.83 quality-adjusted life years among carriers. At a cost of