Complex rearrangements and oncogene amplifications revealed by long-read DNA and RNA sequencing of a breast cancer cell line.

The SK-BR-3 cell line is one of the most important models for HER2+ breast cancers, which affect one in five breast cancer patients. SK-BR-3 is known to be highly rearranged, although much of the variation is in complex and repetitive regions that may be underreported. Addressing this, we sequenced SK-BR-3 using long-read single molecule sequencing from Pacific Biosciences and develop one of the most detailed maps of structural variations (SVs) in a cancer genome available, with nearly 20,000 variants present, most of which were missed by short-read sequencing. Surrounding the important ERBB2 oncogene (also known as HER2), we discover a complex sequence of nested duplications and translocations, suggesting a punctuated progression. Full-length transcriptome sequencing further revealed several novel gene fusions within the nested genomic variants. Combining long-read genome and transcriptome sequencing enables an in-depth analysis of how SVs disrupt the genome and sheds new light on the complex mechanisms involved in cancer genome evolution.

Improved maize reference genome with single-molecule technologies.

Complete and accurate reference genomes and annotations provide fundamental tools for characterization of genetic and functional variation. These resources facilitate the determination of biological processes and support translation of research findings into improved and sustainable agricultural technologies. Many reference genomes for crop plants have been generated over the past decade, but these genomes are often fragmented and missing complex repeat regions. Here we report the assembly and annotation of a reference genome of maize, a genetic and agricultural model species, using single-molecule real-time sequencing and high-resolution optical mapping. Relative to the previous reference genome, our assembly features a 52-fold increase in contig length and notable improvements in the assembly of intergenic spaces and centromeres. Characterization of the repetitive portion of the genome revealed more than 130,000 intact transposable elements, allowing us to identify transposable element lineage expansions that are unique to maize. Gene annotations were updated using 111,000 full-length transcripts obtained by single-molecule real-time sequencing. In addition, comparative optical mapping of two other inbred maize lines revealed a prevalence of deletions in regions of low gene density and maize lineage-specific genes.

A trimming-and-retrieving alignment scheme for reduced representation bisulfite sequencing.

UNLABELLED: Currently available bisulfite sequencing tools frequently suffer from low mapping rates and low methylation calls, especially for data generated from the Illumina sequencer, NextSeq. Here, we introduce a sequential trimming-and-retrieving alignment approach for investigating DNA methylation patterns, which significantly improves the number of mapped reads and covered CpG sites. The method is implemented in an automated analysis toolkit for processing bisulfite sequencing reads. AVAILABILITY AND IMPLEMENTATION: http://mysbfiles.stonybrook.edu/~xuefenwang/software.html and https://github.com/xfwang/BStools.

Whole genome de novo assemblies of three divergent strains of rice, Oryza sativa, document novel gene space of aus and indica.

BACKGROUND: The use of high throughput genome-sequencing technologies has uncovered a large extent of structural variation in eukaryotic genomes that makes important contributions to genomic diversity and phenotypic variation. When the genomes of different strains of a given organism are compared, whole genome resequencing data are typically aligned to an established reference sequence. However, when the reference differs in significant structural ways from the individuals under study, the analysis is often incomplete or inaccurate. RESULTS: Here, we use rice as a model to demonstrate how improvements in sequencing and assembly technology allow rapid and inexpensive de novo assembly of next generation sequence data into high-quality assemblies that can be directly compared using whole genome alignment to provide an unbiased assessment. Using this approach, we are able to accurately assess the "pan-genome" of three divergent rice varieties and document several megabases of each genome absent in the other two. CONCLUSIONS: Many of the genome-specific loci are annotated to contain genes, reflecting the potential for new biological properties that would be missed by standard reference-mapping approaches. We further provide a detailed analysis of several loci associated with agriculturally important traits, including the S5 hybrid sterility locus, the Sub1 submergence tolerance locus, the LRK gene cluster associated with improved yield, and the Pup1 cluster associated with phosphorus deficiency, illustrating the utility of our approach for biological discovery. All of the data and software are openly available to support further breeding and functional studies of rice and other species.

Establishment of highly tumorigenic human colorectal cancer cell line (CR4) with properties of putative cancer stem cells.

BACKGROUND: Colorectal cancer (CRC) has the third highest mortality rates among the US population. According to the most recent concept of carcinogenesis, human tumors are organized hierarchically, and the top of it is occupied by malignant stem cells (cancer stem cells, CSCs, or cancer-initiating cells, CICs), which possess unlimited self-renewal and tumor-initiating capacities and high resistance to conventional therapies. To reflect the complexity and diversity of human tumors and to provide clinically and physiologically relevant cancer models, large banks of characterized patient-derived low-passage cell lines, and especially CIC-enriched cell lines, are urgently needed. PRINCIPAL FINDINGS: Here we report the establishment of a novel CIC-enriched, highly tumorigenic and clonogenic colon cancer cell line, CR4, derived from liver metastasis. This stable cell line was established by combining 3D culturing and 2D culturing in stem cell media, subcloning of cells with particular morphology, co-culture with carcinoma associated fibroblasts (CAFs) and serial transplantation to NOD/SCID mice. Using RNA-Seq complete transcriptome profiling of the tumorigenic fraction of the CR4 cells in comparison to the bulk tumor cells, we have identified about 360 differentially expressed transcripts, many of which represent stemness, pluripotency and resistance to treatment. Majority of the established CR4 cells express common markers of stemness, including CD133, CD44, CD166, EpCAM, CD24 and Lgr5. Using immunocytochemical, FACS and western blot analyses, we have shown that a significant ratio of the CR4 cells express key markers of pluripotency markers, including Sox-2, Oct3/4 and c-Myc. Constitutive overactivation of ABC transporters and NF-kB and absence of tumor suppressors p53 and p21 may partially explain exceptional drug resistance of the CR4 cells. CONCLUSIONS: The highly tumorigenic and clonogenic CIC-enriched CR4 cell line may provide an important new tool to support the discovery of novel diagnostic and/or prognostic biomarkers as well as the development of more effective therapeutic strategies.

Parallel comparison of Illumina RNA-Seq and Affymetrix microarray platforms on transcriptomic profiles generated from 5-aza-deoxy-cytidine treated HT-29 colon cancer cells and simulated datasets.

BACKGROUND: High throughput parallel sequencing, RNA-Seq, has recently emerged as an appealing alternative to microarray in identifying differentially expressed genes (DEG) between biological groups. However, there still exists considerable discrepancy on gene expression measurements and DEG results between the two platforms. The objective of this study was to compare parallel paired-end RNA-Seq and microarray data generated on 5-azadeoxy-cytidine (5-Aza) treated HT-29 colon cancer cells with an additional simulation study. METHODS: We first performed general correlation analysis comparing gene expression profiles on both platforms. An Errors-In-Variables (EIV) regression model was subsequently applied to assess proportional and fixed biases between the two technologies. Then several existing algorithms, designed for DEG identification in RNA-Seq and microarray data, were applied to compare the cross-platform overlaps with respect to DEG lists, which were further validated using qRT-PCR assays on selected genes. Functional analyses were subsequently conducted using Ingenuity Pathway Analysis (IPA). RESULTS: Pearson and Spearman correlation coefficients between the RNA-Seq and microarray data each exceeded 0.80, with 66%~68% overlap of genes on both platforms. The EIV regression model indicated the existence of both fixed and proportional biases between the two platforms. The DESeq and baySeq algorithms (RNA-Seq) and the SAM and eBayes algorithms (microarray) achieved the highest cross-platform overlap rate in DEG results from both experimental and simulated datasets. DESeq method exhibited a better control on the false discovery rate than baySeq on the simulated dataset although it performed slightly inferior to baySeq in the sensitivity test. RNA-Seq and qRT-PCR, but not microarray data, confirmed the expected reversal of SPARC gene suppression after treating HT-29 cells with 5-Aza. Thirty-three IPA canonical pathways were identified by both microarray and RNA-Seq data, 152 pathways by RNA-Seq data only, and none by microarray data only. CONCLUSIONS: These results suggest that RNA-Seq has advantages over microarray in identification of DEGs with the most consistent results generated from DESeq and SAM methods. The EIV regression model reveals both fixed and proportional biases between RNA-Seq and microarray. This may explain in part the lower cross-platform overlap in DEG lists compared to those in detectable genes.

Gene expression in mouse oocytes by RNA-Seq.

Second-generation sequencers such as the Illumina GAIIx make possible the study of all transcribed loci in a genome across an almost endless dynamic range. Although typical protocols call for starting from at least 1 µg of total RNA, this is not possible when studying small tissues or rare cell types. This chapter explains how to prepare Illumina sequencing libraries from mouse oocytes. The protocol is also suitable for mural and cumulus cells, flow sorted or laser captured cells.

Development and application of bovine and porcine oligonucleotide arrays with protein-based annotation.

The design of oligonucleotide sequences for the detection of gene expression in species with disparate volumes of genome and EST sequence information has been broadly studied. However, a congruous strategy has yet to emerge to allow the design of sensitive and specific gene expression detection probes. This study explores the use of a phylogenomic approach to align transcribed sequences to vertebrate protein sequences for the detection of gene families to design genomewide 70-mer oligonucleotide probe sequences for bovine and porcine. The bovine array contains 23,580 probes that target the transcripts of 16,341 genes, about 72% of the total number of bovine genes. The porcine array contains 19,980 probes targeting 15,204 genes, about 76% of the genes in the Ensembl annotation of the pig genome. An initial experiment using the bovine array demonstrates the specificity and sensitivity of the array.

Molecular mechanisms regulating bovine ovarian follicular selection.

Transcription profiling of ovarian follicles. Understanding the mechanisms by which a single follicle is selected for further ovulation is important to control fertility in mammals. However, development of new treatments is limited by our poor understanding of molecular mechanisms regulating follicular selection. Our hypothesis is that genes involved in the control of cell proliferation and apoptosis are differentially regulated during follicular selection. Our objective was to identify these new genes. Bovine follicles were collected and gene expression levels were measured using microarrays. First, follicles were allocated to three groups, according to the time spent from the initiation of follicular wave to surgery (24 H, 36 H, and 48-60 H). Fifty-seven genes are differentially expressed at a false discovery rate of 5%. These genes are involved in the control of lipid metabolism (P-value = 0.0005), cell proliferation (0.007), cell death (0.003), cell morphology (0.003), and immune response (0.003). Follicles were also grouped into four categories, according to the expected time of deviation (early deviation; 8 mm, mid-deviation; 8.5 mm, late deviation; 9 mm, dominant follicles; >or=10 mm). One hundred and twenty eight genes are differentially expressed between these four groups, including genes involved in cell proliferation (0.00002), cell death (0.0006), cell-to-cell signaling (0.003), cell morphology (0.003), lipid metabolism (0.0004), and immune response (0.00007). The expression levels of 10 genes were confirmed using quantitative real time PCR. As expected, we identified new differentially regulated genes involved in the control of cell growth and apoptosis. We also discovered a potential role for immune cells, and in particular macrophages, in follicular selection.

Aberrant DNA methylation in porcine in vitro-, parthenogenetic-, and somatic cell nuclear transfer-produced blastocysts.

Early embryonic development in the pig requires DNA methylation remodeling of the maternal and paternal genomes. Aberrant remodeling, which can be exasperated by in vitro technologies, is detrimental to development and can result in physiological and anatomic abnormalities in the developing fetus and offspring. Here, we developed and validated a microarray based approach to characterize on a global scale the CpG methylation profiles of porcine gametes and blastocyst stage embryos. The relative methylation in the gamete and blastocyst samples showed that 18.5% (921/4,992) of the DNA clones were found to be significantly different (P < 0.01) in at least one of the samples. Furthermore, for the different blastocyst groups, the methylation profile of the in vitro-produced blastocysts was less similar to the in vivo-produced blastocysts as compared to the parthenogenetic- and somatic cell nuclear transfer (SCNT)-produced blastocysts. The microarray results were validated by using bisulfite sequencing for 12 of the genomic regions in liver, sperm, and in vivo-produced blastocysts. These results suggest that a generalized change in global methylation is not responsible for the low developmental potential of blastocysts produced by using in vitro techniques. Instead, the appropriate methylation of a relatively small number of genomic regions in the early embryo may enable early development to occur.

Correlation of developmental differences of nuclear transfer embryos cells to the methylation profiles of nuclear transfer donor cells in Swine.

Methylation of DNA is the most commonly studied epigenetic mechanism of developmental competence and somatic cell nuclear transfer (SCNT). Previous studies of epigenetics and the SCNT procedures have examined the effects of different culture media on donor cells and reconstructed embryos, and the methylation status of specific genes in the fetus or live offspring. Here we used a microarray based approach to identify the methylation profiles of SCNT donor cells including three clonal porcine fetal fibroblast-like cell sublines and adult somatic cells selected from kidney and mammary tissues. The methylation profiles of the donor cells were then analyzed with respect to their ability to direct development to the blastocyst stage after nuclear transfer. Clonal cell lines A2, A7 and A8 had blastocyst rates of 11.7%(a), 16.7%(ab) and 20.0%(b), respectively ((ab) p < 0.05). Adult somatic cells included kidney, mammary (large), and mammary (small) also had different blastocyst rates (ab p < 0.05) of 4.2% (a), 10.7% (ab) and 18.3% (b), respectively. For clonal donor cells and for adult somatic cell groups the donor cells with the highest blastocyst rates also had methylation profiles with the lowest similarity to the methylation profiles of the in vivo-produced blastocysts. Conversely, the donor cells with the lowest blastocyst rates had methylation profiles with the highest similarity to the methylation profiles of the in vivo-produced blastocysts. Our findings show there is an inverse correlation to the similarity of the methylation profiles of the donor cells and the in vivo-produced embryos, and to the blastocyst rates following SCNT.

A microarray analysis for genes regulated by interferon-tau in ovine luminal epithelial cells.

Interferon-tau (IFNT) is released by preimplantation conceptuses of ruminant species and prepares the mother for pregnancy. Although one important function is to protect the corpus luteum from the luteolytic activity of prostaglandin-F 2alpha, IFNT most likely regulates a range of other physiological processes in endometrium. Here, an immortalized cell line from ovine uterine luminal epithelial cells was treated with IFNT for either 8 or 24 h. RNA was subjected to cDNA microarray analysis, with RNA from untreated cells as the reference standard. Of 15 634 genes, 1274 (8%) were IFNT responsive at P<0.01 and 585 at P<0.001 to at least one treatment. Of the latter, 356 were up-regulated and 229 down-regulated. Increasing IFNT concentrations from 10 ng/ml to 10 microg/ml had minor effects, and most genes up- or down-regulated at 8 h were regulated similarly at 24 h. Although IFNT influences many genes implicated in antiviral activity and apoptosis, its action also likely regulates prostaglandin metabolism, growth factors and their receptors, apoptosis and the nuclear factor (NF)-kappaB cascade, extracellular matrix accretion, angiogenesis, blood coagulation, and inflammation. In particular, it increased mRNA concentrations of genes related to the vascular endothelial growth factor R2 pathway of angiogenesis and down-regulated ones associated with hypoxia. Two genes implicated in the antiluteolytic actions of IFNT (encoding cyclooxygenase-2 and the oxytocin receptor respectively) were down-regulated in response to all treatments. IFNT targets a complex range of physiological processes during the establishment of pregnancy.

Expression profiling of heat stress effects on mice fed ergot alkaloids.

Fescue toxicosis affects wild and domestic animals consuming ergot alkaloids contained in tall fescue forage infected with the endophytic fungus, Neotyphodium coenophialum. When animals are consuming infected fescue (E+) forage during periods of elevated ambient temperatures (summer), a range of phenotypic disorders collectively called summer slump is observed. It is characterized by hyperthermia, with an accompanying decrease in feed intake, growth, milk yield, and reproductive fitness. Laboratory mice also exhibit symptoms of fescue toxicosis at thermoneutral (TN) temperature, as indicated by reduced growth rate and reproductive fitness. Our goal was to characterize the differences in gene expression in liver of mice exposed to summer-type heat stress (HS) and E+ when compared to mice fed E+ at TN temperature. Mice were fed E+ diet under HS (34 +/- 1 degrees C; n = 13; E+HS) or TN conditions (24 +/- 1 degrees C; n = 14; E+TN) for a period of 2 weeks between 47 and 60 days of age. Genes differentially expressed between E+HS versus E+TN were identified using DNA microarrays. Forty-one genes were differentially expressed between treatment groups. Expressions of eight genes were measured using quantitative real-time PCR. Genes coding for phase I detoxification enzymes were upregulated in E+HS mouse liver. This detoxification pathway is known to produce reactive oxidative species. We observed an upregulation of genes involved in the protection against reactive oxidative species. Key genes involved in de novo lipogenesis and lipid transport were also upregulated. Finally, genes involved in DNA damage control and unfolded protein responses were downregulated.

Luteinization of porcine preovulatory follicles leads to systematic changes in follicular gene expression.

The LH surge initiates the luteinization of preovulatory follicles and causes hormonal and structural changes that ultimately lead to ovulation and the formation of corpora lutea. The objective of the study was to examine gene expression in ovarian follicles (n = 11) collected from pigs (Sus scrofa domestica) approaching estrus (estrogenic preovulatory follicle; n = 6 follicles from two sows) and in ovarian follicles collected from pigs on the second day of estrus (preovulatory follicles that were luteinized but had not ovulated; n = 5 follicles from two sows). The follicular status within each follicle was confirmed by follicular fluid analyses of estradiol and progesterone ratios. Microarrays were made from expressed sequence tags that were isolated from cDNA libraries of porcine ovary. Gene expression was measured by hybridization of fluorescently labeled cDNA (preovulatory estrogenic or -luteinized) to the microarray. Microarray analyses detected 107 and 43 genes whose expression was decreased or increased (respectively) during the transition from preovulatory estrogenic to -luteinized (P<0.01). Cells within preovulatory estrogenic follicles had a gene-expression profile of proliferative and metabolically active cells that were responding to oxidative stress. Cells within preovulatory luteinized follicles had a gene-expression profile of nonproliferative and migratory cells with angiogenic properties. Approximately, 40% of the discovered genes had unknown function.

Effect of interferon-tau administration on endometrium of nonpregnant ewes: a comparison with pregnant ewes.

In ruminants, conceptus interferon-tau (IFNT) alters maternal physiology to accommodate a pregnancy. We hypothesized that the effectiveness of IFNT on extending corpus luteum (CL) life span in nonpregnant ewes would depend upon the dose and manner of administration and would be correlated with the response in gene expression in endometrium. We anticipated that IFNT, whether administered im or by uterine infusion, would mimic changes observed in pregnancy. Ewes were assigned to five treatments: 1) uterine infusion of saline; 2) uterine infusion of ovine IFNT4 (200 microg/d); 3) saline im injection; 4) im injection of IFNT4 at low dose (200 microg/d); and 5) high dose (2 mg/d). CL life span was increased in groups 2 and 5, but not in 1, 3, and 4. Endometrial RNA extracted from groups 1-5 on d 14 and from d 14 pregnant and nonbred (cyclic) ewes was used to assess expression of 70 genes on microarrays. When pregnant and cyclic ewes were compared, 30 genes were up-regulated and nine down-regulated during pregnancy. Responses were slightly less in groups 2 and 5 but were much lower in group 4. The majority of the highly up-regulated genes were associated with antiviral responses. Those down-regulated included ones for IGF-II, hypoxia-inducible factor 1alpha, oxytocin receptor, prostaglandin F synthase, and cyclooxygenase-2. Quantitative PCR for selected genes confirmed these data and revealed that similar gene expression changes occurred in the CL of pregnant and group 2 ewes. IFNT treatment mimics pregnancy, but relatively high doses of im-injected IFNT are required to elicit a full endometrial response.

Identification and characterization of two novel splicing isoforms of human estrogen-related receptor beta.

CONTEXT: Estrogen-related receptor beta (ERRbeta) was one of the first two orphan nuclear receptors reported and is believed to play important roles in estrogen-regulated pathways. Embryo lethality of ERRbeta-null mice indicated that ERRbeta is essential for embryo development. OBJECTIVE: Two novel splicing isoforms of human (h) ERRbeta, hERRbeta2-Delta10 and short-form hERRbeta, were identified during the cloning of previously reported hERRbeta-hERRbeta2. We aim to investigate the functional differences of these three human ERRbeta-splicing isoforms. RESULTS AND CONCLUSIONS: A genomic sequence comparison within and flanking the ERRbeta genes of eight species demonstrated that short-form hERRbeta lacks an F domain and is the matched homolog of mouse and rat ERRbeta proteins in humans. However, hERRbeta2-Delta10 and the previously reported hERRbeta2 isoforms are primate specific. RT-PCR analysis showed that short-form hERRbeta has a wide distribution in the 24 of 27 human tissues and cell lines tested, whereas hERRbeta2 and hERRbeta2-Delta10 were only expressed in testis and kidney. The three human ERRbeta-splicing isoforms have different transcriptional activities when measured on an estrogen response element-driven luciferase reporter in transfection assays. The localization of a nuclear localization signal of short-form hERRbeta was also determined. Interestingly, the F domain of hERRbeta2 alters the function of the nuclear localization signal. Therefore, the ERRbeta isoforms are likely to have diverse biological functions in vivo, and characterizing the three isoforms of ERRbeta will lead to an understanding of the multiple levels of gene regulation involved in steroid receptor-signaling pathways in humans and may provide novel therapeutic targets for human diseases.