Core circadian clock transcription factor BMAL1 regulates mammary epithelial cell growth, differentiation, and milk component synthesis

The role the mammary epithelial circadian clock plays in gland development and lactation is unknown. We hypothesized that mammary epithelial clocks function to regulate mammogenesis and lactogenesis, and propose the core clock transcription factor BMAL1:CLOCK regulates genes that control mammary epithelial development and milk synthesis. Our objective was to identify transcriptional targets of BMAL1 in undifferentiated (UNDIFF) and lactogen differentiated (DIFF) mammary epithelial cells (HC11) using ChIP-seq. Ensembl gene IDs with the nearest transcriptional start site to ChIP-seq peaks were explored as potential targets, and represented 846 protein coding genes common to UNDIFF and DIFF cells and 2773 unique to DIFF samples. Genes with overlapping peaks between samples (1343) enriched cell-cell adhesion, membrane transporters and lipid metabolism categories. To functionally verify targets, an HC11 line with Bmal1 gene knocked out (BMAL1-KO) using CRISPR-CAS was created. BMAL1-KO cultures had lower cell densities over an eight-day growth curve, which was associated with increased (p<0.05) levels of reactive oxygen species and lower expression of superoxide dismutase 3 (Sod3). RT-qPCR analysis also found lower expression of the putative targets, prolactin receptor (Prlr), Ppara, and beta-casein (Csn2). Findings support our hypothesis and highlight potential importance of clock in mammary development and substrate transport.

Characterization of the Gray Whale Eschrichtius robustus Genome and a Genotyping Array Based on Single-Nucleotide Polymorphisms in Candidate Genes

Genetic and genomic approaches have much to offer in terms of ecology, evolution, and conservation. To better understand the biology of the gray whale Eschrichtius robustus (Lilljeborg, 1861), we sequenced the genome and produced an assembly that contains ∼95% of the genes known to be highly conserved among eukaryotes. From this assembly, we annotated 22,711 genes and identified 2,057,254 single-nucleotide polymorphisms (SNPs). Using this assembly, we generated a curated list of candidate genes potentially subject to strong natural selection, including genes associated with osmoregulation, oxygen binding and delivery, and other aspects of marine life. From these candidate genes, we queried 92 autosomal protein-coding markers with a panel of 96 SNPs that also included 2 sexing and 2 mitochondrial markers. Genotyping error rates, calculated across loci and across 69 intentional replicate samples, were low (0.021%), and observed heterozygosity was 0.33 averaged over all autosomal markers. This level of variability provides substantial discriminatory power across loci (mean probability of identity of 1.6 × 10^-25 and mean probability of exclusion >0.999 with neither parent known), indicating that these markers provide a powerful means to assess parentage and relatedness in gray whales. We found 29 unique multilocus genotypes represented among our 36 biopsies (indicating that we inadvertently sampled 7 whales twice). In total, we compiled an individual data set of 28 western gray whales (WGSs) and 1 presumptive eastern gray whale (EGW). The lone EGW we sampled was no more or less related to the WGWs than expected by chance alone. The gray whale genomes reported here will enable comparative studies of natural selection in cetaceans, and the SNP markers should be highly informative for future studies of gray whale evolution, population structure, demography, and relatedness.

An epigenetic switch regulates de novo DNA methylation at a subset of pluripotency gene enhancers during embryonic stem cell differentiation

Coordinated regulation of gene expression that involves activation of lineage specific genes and repression of pluripotency genes drives differentiation of embryonic stem cells (ESC). For complete repression of pluripotency genes during ESC differentiation, chromatin at their enhancers is silenced by the activity of the Lsd1-Mi2/NuRD complex. The mechanism/s that regulate DNA methylation at these enhancers are largely unknown. Here, we investigated the affect of the Lsd1-Mi2/NuRD complex on the dynamic regulatory switch that induces the local interaction of histone tails with the Dnmt3 ATRX-DNMT3-DNMT3L (ADD) domain, thus promoting DNA methylation at the enhancers of a subset of pluripotency genes. This is supported by previous structural studies showing a specific interaction between Dnmt3-ADD domain with H3K4 unmethylated histone tails that is disrupted by histone H3K4 methylation and histone acetylation. Our data suggest that Dnmt3a activity is triggered by Lsd1-Mi2/NuRD-mediated histone deacetylation and demethylation at these pluripotency gene enhancers when they are inactivated during mouse ESC differentiation. Using Dnmt3 knockout ESCs and the inhibitors of Lsd1 and p300 histone modifying enzymes during differentiation of E14Tg2A and ZHBTc4 ESCs, our study systematically reveals this mechanism and establishes that Dnmt3a is both reader and effector of the epigenetic state at these target sites.

RNA-seq analysis identifies potential modulators of gravity response in spores of Ceratopteris (Parkeriaceae): evidence for modulation by calcium pumps and apyrase activity

PREMISE OF THE STUDY

Gravity regulates the magnitude and direction of a trans-cell calcium current in germinating spores of Ceratopteris richardii. Blocking this current with nifedipine blocks the spore's downward polarity alignment, a polarization that is fixed by gravity ∼10 h after light induces the spores to germinate. RNA-seq analysis at 10 h was used to identify genes potentially important for the gravity response. The data set will be valuable for other developmental and phylogenetic studies.

METHODS

De novo Newbler assembly of 958 527 reads from Roche 454 sequencing was executed. The sequences were identified and analyzed using in silico methods. The roles of endomembrane Ca(2+)-ATPase pumps and apyrases in the gravity response were further tested using pharmacological agents.

KEY RESULTS

Transcripts related to calcium signaling and ethylene biosynthesis were identified as notable constituents of the transcriptome. Inhibiting the activity of endomembrane Ca(2+)-ATPase pumps with 2,5-di-(t-butyl)-1,4-hydroquinone diminished the trans-cell current, but increased the orientation of the polar axis to gravity. The effects of applied nucleotides and purinoceptor antagonists gave novel evidence implicating extracellular nucleotides as regulators of the gravity response in these fern spores.

CONCLUSIONS

In addition to revealing general features of the transcriptome of germinating spores, the results highlight a number of calcium-responsive and light-receptive transcripts. Pharmacologic assays indicate endomembrane Ca(2+)-ATPases and extracellular nucleotides may play regulatory roles in the gravity response of Ceratopteris spores.

Orthologous comparisons of the Hd1 region across genera reveal Hd1 gene lability within diploid Oryza species and disruptions to microsynteny in Sorghum

Heading date is one of the most important quantitative traits responsible for the domestication of rice. We compared a 155-kb reference segment of the Oryza sativa ssp. japonica cv. Nipponbare genome surrounding Hd1, a major heading date gene in rice, with orthologous regions from nine diploid Oryza species that diverged over a relatively short time frame (∼16 My) to study sequence evolution around a domestication locus. The orthologous Hd1 region from Sorghum bicolor was included to compare and contrast the evolution in a more distant relative of rice. Consistent with other observations at the adh1/adh2, monoculm1, and sh2/a1 loci in grass species, we found high gene colinearity in the Hd1 region amidst size differences that were lineage specific and long terminal repeat retrotransposon driven. Unexpectedly, the Hd1 gene was deleted in O. glaberrima, whereas the O. rufipogon and O. punctata copies had degenerative mutations, suggesting that other heading date loci might compensate for the loss or nonfunctionality of Hd1 in these species. Compared with the japonica Hd1 region, the orthologous region in sorghum exhibited micro-rearrangements including gene translocations, seven additional genes, and a gene triplication and truncation event predating the divergence from Oryza.

Detailed analysis of a contiguous 22-Mb region of the maize genome

Most of our understanding of plant genome structure and evolution has come from the careful annotation of small (e.g., 100 kb) sequenced genomic regions or from automated annotation of complete genome sequences. Here, we sequenced and carefully annotated a contiguous 22 Mb region of maize chromosome 4 using an improved pseudomolecule for annotation. The sequence segment was comprehensively ordered, oriented, and confirmed using the maize optical map. Nearly 84% of the sequence is composed of transposable elements (TEs) that are mostly nested within each other, of which most families are low-copy. We identified 544 gene models using multiple levels of evidence, as well as five miRNA genes. Gene fragments, many captured by TEs, are prevalent within this region. Elimination of gene redundancy from a tetraploid maize ancestor that originated a few million years ago is responsible in this region for most disruptions of synteny with sorghum and rice. Consistent with other sub-genomic analyses in maize, small RNA mapping showed that many small RNAs match TEs and that most TEs match small RNAs. These results, performed on approximately 1% of the maize genome, demonstrate the feasibility of refining the B73 RefGen_v1 genome assembly by incorporating optical map, high-resolution genetic map, and comparative genomic data sets. Such improvements, along with those of gene and repeat annotation, will serve to promote future functional genomic and phylogenomic research in maize and other grasses.

Comparative physical mapping between Oryza sativa (AA genome type) and O. punctata (BB genome type)

A comparative physical map of the AA genome (Oryza sativa) and the BB genome (O. punctata) was constructed by aligning a physical map of O. punctata, deduced from 63,942 BAC end sequences (BESs) and 34,224 fingerprints, onto the O. sativa genome sequence. The level of conservation of each chromosome between the two species was determined by calculating a ratio of BES alignments. The alignment result suggests more divergence of intergenic and repeat regions in comparison to gene-rich regions. Further, this characteristic enabled localization of heterochromatic and euchromatic regions for each chromosome of both species. The alignment identified 16 locations containing expansions, contractions, inversions, and transpositions. By aligning 40% of the punctata BES on the map, 87% of the punctata FPC map covered 98% of the O. sativa genome sequence. The genome size of O. punctata was estimated to be 8% larger than that of O. sativa with individual chromosome differences of 1.5-16.5%. The sum of expansions and contractions observed in regions >500 kb were similar, suggesting that most of the contractions/expansions contributing to the genome size difference between the two species are small, thus preserving the macro-collinearity between these species, which diverged approximately 2 million years ago.

Integration of hybridization-based markers (overgos) into physical maps for comparative and evolutionary explorations in the genus Oryza and in Sorghum

BACKGROUND

With the completion of the genome sequence for rice (Oryza sativa L.), the focus of rice genomics research has shifted to the comparison of the rice genome with genomes of other species for gene cloning, breeding, and evolutionary studies. The genus Oryza includes 23 species that shared a common ancestor 8-10 million years ago making this an ideal model for investigations into the processes underlying domestication, as many of the Oryza species are still undergoing domestication. This study integrates high-throughput, hybridization-based markers with BAC end sequence and fingerprint data to construct physical maps of rice chromosome 1 orthologues in two wild Oryza species. Similar studies were undertaken in Sorghum bicolor, a species which diverged from cultivated rice 40-50 million years ago.

RESULTS

Overgo markers, in conjunction with fingerprint and BAC end sequence data, were used to build sequence-ready BAC contigs for two wild Oryza species. The markers drove contig merges to construct physical maps syntenic to rice chromosome 1 in the wild species and provided evidence for at least one rearrangement on chromosome 1 of the O. sativa versus Oryza officinalis comparative map. When rice overgos were aligned to available S. bicolor sequence, 29% of the overgos aligned with three or fewer mismatches; of these, 41% gave positive hybridization signals. Overgo hybridization patterns supported colinearity of loci in regions of sorghum chromosome 3 and rice chromosome 1 and suggested that a possible genomic inversion occurred in this syntenic region in one of the two genomes after the divergence of S. bicolor and O. sativa.

CONCLUSION

The results of this study emphasize the importance of identifying conserved sequences in the reference sequence when designing overgo probes in order for those probes to hybridize successfully in distantly related species. As interspecific markers, overgos can be used successfully to construct physical maps in species which diverged less than 8 million years ago, and can be used in a more limited fashion to examine colinearity among species which diverged as much as 40 million years ago. Additionally, overgos are able to provide evidence of genomic rearrangements in comparative physical mapping studies.

Analysis of expressed sequence tags from Gibberella zeae (anamorph Fusarium graminearum)

Gibberella zeae is a broad host range pathogen that infects many crop plants, including wheat and barley, and causes head blight and rot diseases throughout the world. To better understand fungal development and pathogenicity, we have generated 7996 ESTs from three cDNA libraries. Two libraries were generated from carbon-(C-) and nitrogen- (N-) starved mycelia and one library was generated from cultures of maturing perithecia (P). In other fungal pathogens, starvation conditions have been shown to act as cues to induce infection-related gene expression. To assign putative function to cDNAs, sequences were initially assembled using StackPack. The estimated total number of genes identified from the three EST databases was 2110: 1088 contigs and 1022 singleton sequences. These 2110 sequences were compared to a yeast protein sequence reference set and to the GenBank nonredundant database using BLASTX. Based on presumptive gene function identified by this process, we found that the two starved cultures had similar, but not identical, patterns of gene expression, whereas the developmental cultures were distinct in their pattern of expression. Of the three libraries, the perithecium library had the greatest percentage (46%) of ESTS falling into the "unclassified" category. Homologues of some known fungal virulence or pathogenicity factors were found primarily in the N- and C-libraries. Comparisons also were made with ESTs from the related fungi, Neurospora crassa and Magnaporthe grisea and the genomic sequence of N. crassa.

Lymphoepithelioma-like carcinoma of the uterine cervix: a case report studied by in situ hybridization and polymerase chain reaction for Epstein-Barr virus

Lymphoepithelioma-like carcinomas have been reported outside the nasopharynx in many sites, including the uterine cervix. The association with the Epstein-Barr virus in the latter site is still controversial. To date, Epstein-Barr virus genome has only been demonstrated in Asian patients. We report a case of lymphoepithelioma-like carcinoma of the uterine cervix in a white woman in whom the Epstein-Barr virus infection was tested for by in situ hybridization and polymerase chain reaction. The results of both techniques were negative. Our case and a review of the literature support the contention that cervical lymphoepithelioma-like carcinoma is not associated with Epstein-Barr virus infection in non-Asian patients.

The y1 gene of maize codes for phytoene synthase

The cloned y1 locus of maize was sequenced and found to encode phytoene synthase. Different "wild-type" alleles of the locus were found to differ by the insertion of transposable elements in their promoter and polyA addition regions, and by the length of a CCA tandem repeat series, without any obvious effect on function of the gene. A dominant Y1 ("wild-type") allele was observed to be expressed at highest levels in the seedling but also in the embryo and endosperm. The Mu3 transposable element insertion responsible for a pastel allele of y1, which gives lowered levels of carotenoids in the endosperm of kernels and seedlings grown at high temperatures, was located in the 5' end of the gene. Although the size of the transcript from this y1 mutation suggests that the Mu3 element provides the promoter for this allele, leaf tissue in this mutant line contained approximately normal amounts of y1 mRNA. A recessive allele of y1, which conditions normal levels of carotenoids in the embryo and seedling, but almost no carotenoids in the endosperm, was found to accumulate normal amounts of y1 mRNA in the seedling and embryo, while y1 transcripts were not detected in the endosperm.