Library preparation and data analysis packages for rapid genome sequencing

The arginine methyltransferase Prmt1 coordinates the germline arginine methylome essential for spermatogonial homeostasis and male fertility

Arginine methylation, catalyzed by the protein arginine methyltransferases (PRMTs), is a common post-translational protein modification (PTM) that is engaged in a plethora of biological events. However, little is known about how the methylarginine-directed signaling functions in germline development. In this study, we discover that Prmt1 is predominantly distributed in the nuclei of spermatogonia but weakly in the spermatocytes throughout mouse spermatogenesis. By exploiting a combination of three Cre-mediated Prmt1 knockout mouse lines, we unravel that Prmt1 is essential for spermatogonial establishment and maintenance, and that Prmt1-catalyzed asymmetric methylarginine coordinates inherent transcriptional homeostasis within spermatogonial cells. In conjunction with high-throughput CUT&Tag profiling and modified mini-bulk Smart-seq2 analyses, we unveil that the Prmt1-deposited H4R3me2a mark is permissively enriched at promoter and exon/intron regions, and sculpts a distinctive transcriptomic landscape as well as the alternative splicing pattern, in the mouse spermatogonia. Collectively, our study provides the genetic and mechanistic evidence that connects the Prmt1-deposited methylarginine signaling to the establishment and maintenance of a high-fidelity transcriptomic identity in orchestrating spermatogonial development in the mammalian germline.

The Neurospora Transcription Factor ADV-1 Transduces Light Signals and Temporal Information to Control Rhythmic Expression of Genes Involved in Cell Fusion

Light and the circadian clock have a profound effect on the biology of organisms through the regulation of large sets of genes. Toward understanding how light and the circadian clock regulate gene expression, we used genome-wide approaches to identify the direct and indirect targets of the light-responsive and clock-controlled transcription factor ADV-1 in Neurospora crassa A large proportion of ADV-1 targets were found to be light- and/or clock-controlled, and enriched for genes involved in development, metabolism, cell growth, and cell fusion. We show that ADV-1 is necessary for transducing light and/or temporal information to its immediate downstream targets, including controlling rhythms in genes critical to somatic cell fusion. However, while ADV-1 targets are altered in predictable ways in Δadv-1 cells in response to light, this is not always the case for rhythmic target gene expression. These data suggest that a complex regulatory network downstream of ADV-1 functions to generate distinct temporal dynamics of target gene expression relative to the central clock mechanism.

Histone modifications rather than the novel regional centromeres of Zymoseptoria tritici distinguish core and accessory chromosomes

BACKGROUND

Supernumerary chromosomes have been found in many organisms. In fungi, these "accessory" or "dispensable" chromosomes are present at different frequencies in populations and are usually characterized by higher repetitive DNA content and lower gene density when compared to the core chromosomes. In the reference strain of the wheat pathogen, Zymoseptoria tritici, eight discrete accessory chromosomes have been found. So far, no functional role has been assigned to these chromosomes; however, they have existed as separate entities in the karyotypes of Zymoseptoria species over evolutionary time. In this study, we addressed what-if anything-distinguishes the chromatin of accessory chromosomes from core chromosomes. We used chromatin immunoprecipitation combined with high-throughput sequencing ("ChIP-seq") of DNA associated with the centromere-specific histone H3, CENP-A (CenH3), to identify centromeric DNA, and ChIP-seq with antibodies against dimethylated H3K4, trimethylated H3K9 and trimethylated H3K27 to determine the relative distribution and proportion of euchromatin, obligate and facultative heterochromatin, respectively.

RESULTS

Centromeres of the eight accessory chromosomes have the same sequence composition and structure as centromeres of the 13 core chromosomes and they are of similar length. Unlike those of most other fungi, Z. tritici centromeres are not composed entirely of repetitive DNA; some centromeres contain only unique DNA sequences, and bona fide expressed genes are located in regions enriched with CenH3. By fluorescence microscopy, we showed that centromeres of Z. tritici do not cluster into a single chromocenter during interphase. We found dramatically higher enrichment of H3K9me3 and H3K27me3 on the accessory chromosomes, consistent with the twofold higher proportion of repetitive DNA and poorly transcribed genes. In contrast, no single histone modification tested here correlated with the distribution of centromeric nucleosomes.

CONCLUSIONS

All centromeres are similar in length and composed of a mixture of unique and repeat DNA, and most contain actively transcribed genes. Centromeres, subtelomeric regions or telomere repeat length cannot account for the differences in transfer fidelity between core and accessory chromosomes, but accessory chromosomes are greatly enriched in nucleosomes with H3K27 trimethylation. Genes on accessory chromosomes appear to be silenced by trimethylation of H3K9 and H3K27.

Sequence-based ultra-dense genetic and physical maps reveal structural variations of allopolyploid cotton genomes

Background

SNPs are the most abundant polymorphism type, and have been explored in many crop genomic studies, including rice and maize. SNP discovery in allotetraploid cotton genomes has lagged behind that of other crops due to their complexity and polyploidy. In this study, genome-wide SNPs are detected systematically using next-generation sequencing and efficient SNP genotyping methods, and used to construct a linkage map and characterize the structural variations in polyploid cotton genomes.

Results

We construct an ultra-dense inter-specific genetic map comprising 4,999,048 SNP loci distributed unevenly in 26 allotetraploid cotton linkage groups and covering 4,042 cM. The map is used to order tetraploid cotton genome scaffolds for accurate assembly of G. hirsutum acc. TM-1. Recombination rates and hotspots are identified across the cotton genome by comparing the assembled draft sequence and the genetic map. Using this map, genome rearrangements and centromeric regions are identified in tetraploid cotton by combining information from the publicly-available G. raimondii genome with fluorescent in situ hybridization analysis.

Conclusions

We report the genotype-by-sequencing method used to identify millions of SNPs between G. hirsutum and G. barbadense. We construct and use an ultra-dense SNP map to correct sequence mis-assemblies, merge scaffolds into pseudomolecules corresponding to chromosomes, detect genome rearrangements, and identify centromeric regions in allotetraploid cottons. We find that the centromeric retro-element sequence of tetraploid cotton derived from the D subgenome progenitor might have invaded the A subgenome centromeres after allotetrapolyploid formation. This study serves as a valuable genomic resource for genetic research and breeding of cotton.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0678-1) contains supplementary material, which is available to authorized users.

Chromatin analyses of Zymoseptoria tritici: Methods for chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq)

The presence or absence of specific transcription factors, chromatin remodeling machineries, chromatin modification enzymes, post-translational histone modifications and histone variants all play crucial roles in the regulation of pathogenicity genes. Chromatin immunoprecipitation (ChIP) followed by high-throughput sequencing (ChIP-seq) provides an important tool to study genome-wide protein-DNA interactions to help understand gene regulation in the context of native chromatin. ChIP-seq is a convenient in vivo technique to identify, map and characterize occupancy of specific DNA fragments with proteins against which specific antibodies exist or which can be epitope-tagged in vivo. We optimized existing ChIP protocols for use in the wheat pathogen Zymoseptoria tritici and closely related sister species. Here, we provide a detailed method, underscoring which aspects of the technique are organism-specific. Library preparation for Illumina sequencing is described, as this is currently the most widely used ChIP-seq method. One approach for the analysis and visualization of representative sequence is described; improved tools for these analyses are constantly being developed. Using ChIP-seq with antibodies against H3K4me2, which is considered a mark for euchromatin or H3K9me3 and H3K27me3, which are considered marks for heterochromatin, the overall distribution of euchromatin and heterochromatin in the genome of Z. tritici can be determined. Our ChIP-seq protocol was also successfully applied to Z. tritici strains with high levels of melanization or aberrant colony morphology, and to different species of the genus (Z. ardabiliae and Z. pseudotritici), suggesting that our technique is robust. The methods described here provide a powerful framework to study new aspects of chromatin biology and gene regulation in this prominent wheat pathogen.

Neurospora importin α is required for normal heterochromatic formation and DNA methylation

Heterochromatin and associated gene silencing processes play roles in development, genome defense, and chromosome function. In many species, constitutive heterochromatin is decorated with histone H3 tri-methylated at lysine 9 (H3K9me3) and cytosine methylation. In Neurospora crassa, a five-protein complex, DCDC, catalyzes H3K9 methylation, which then directs DNA methylation. Here, we identify and characterize a gene important for DCDC function, dim-3 (defective in methylation-3), which encodes the nuclear import chaperone NUP-6 (Importin α). The critical mutation in dim-3 results in a substitution in an ARM repeat of NUP-6 and causes a substantial loss of H3K9me3 and DNA methylation. Surprisingly, nuclear transport of all known proteins involved in histone and DNA methylation, as well as a canonical transport substrate, appear normal in dim-3 strains. Interactions between DCDC members also appear normal, but the nup-6dim-3 allele causes the DCDC members DIM-5 and DIM-7 to mislocalize from heterochromatin and NUP-6^dim-3 itself is mislocalized from the nuclear envelope, at least in conidia. GCN-5, a member of the SAGA histone acetyltransferase complex, also shows altered localization in dim-3, raising the possibility that NUP-6 is necessary to localize multiple chromatin complexes following nucleocytoplasmic transport.

The epigenetic information contained in chromatin is essential for development of higher organisms, and if misregulated, can lead to the unregulated growth associated with human cancers. Chromatin is typically classified into two basic types: gene-rich 'euchromatin', and gene-poor heterochromatin, which is also rich in repeated DNA and 'repressive chromatin marks'. As in humans and eukaryotes generally, heterochromatin in Neurospora crassa is decorated with DNA methylation and histone H3 lysine 9 (H3K9) methylation, but unlike the case in mammals, loss of these epigenetic marks does not compromise viability. In Neurospora, the DCDC, a five-member Cul4-based protein complex, trimethylates H3K9. Little information is available on the regulation of DCDC or similar complexes in other organisms. Using forward genetics, we identified a novel role for Importin α (NUP-6) for the function of DCDC. Although NUP-6 typically functions in nucleocytoplasmic transport, the dim-3 strain, which contains an altered nup-6 gene that reduces DNA methylation and H3K9me3, shows normal nuclear transport of the heterochromatin machinery and a canonical transport substrate. Two DCDC members are mislocalized from heterochromatin in the dim-3 mutant, signifying that NUP-6 may be important for targeting key proteins to incipient heterochromatic DNA. The euchromatic complex SAGA has increased euchromatin localization in dim-3, suggesting that NUP-6 may localize multiple chromatin complexes to sub-nuclear genomic targets.

Two histone deacetylases, FfHda1 and FfHda2, are important for Fusarium fujikuroi secondary metabolism and virulence

Histone modifications are crucial for the regulation of secondary metabolism in various filamentous fungi. Here we studied the involvement of histone deacetylases (HDACs) in secondary metabolism in the phytopathogenic fungus Fusarium fujikuroi, a known producer of several secondary metabolites, including phytohormones, pigments, and mycotoxins. Deletion of three Zn(2+)-dependent HDAC-encoding genes, ffhda1, ffhda2, and ffhda4, indicated that FfHda1 and FfHda2 regulate secondary metabolism, whereas FfHda4 is involved in developmental processes but is dispensable for secondary-metabolite production in F. fujikuroi. Single deletions of ffhda1 and ffhda2 resulted not only in an increase or decrease but also in derepression of metabolite biosynthesis under normally repressing conditions. Moreover, double deletion of both the ffhda1 and ffhda2 genes showed additive but also distinct phenotypes with regard to secondary-metabolite biosynthesis, and both genes are required for gibberellic acid (GA)-induced bakanae disease on the preferred host plant rice, as Δffhda1 Δffhda2 mutants resemble the uninfected control plant. Microarray analysis with a Δffhda1 mutant that has lost the major HDAC revealed differential expression of secondary-metabolite gene clusters, which was subsequently verified by a combination of chemical and biological approaches. These results indicate that HDACs are involved not only in gene silencing but also in the activation of some genes. Chromatin immunoprecipitation with the Δffhda1 mutant revealed significant alterations in the acetylation state of secondary-metabolite gene clusters compared to the wild type, thereby providing insights into the regulatory mechanism at the chromatin level. Altogether, manipulation of HDAC-encoding genes constitutes a powerful tool to control secondary metabolism in filamentous fungi.

The Fusarium graminearum histone H3 K27 methyltransferase KMT6 regulates development and expression of secondary metabolite gene clusters

The cereal pathogen Fusarium graminearum produces secondary metabolites toxic to humans and animals, yet coordinated transcriptional regulation of gene clusters remains largely a mystery. By chromatin immunoprecipitation and high-throughput DNA sequencing (ChIP-seq) we found that regions with secondary metabolite clusters are enriched for trimethylated histone H3 lysine 27 (H3K27me3), a histone modification associated with gene silencing. H3K27me3 was found predominantly in regions that lack synteny with other Fusarium species, generally subtelomeric regions. Di- or trimethylated H3K4 (H3K4me2/3), two modifications associated with gene activity, and H3K27me3 are predominantly found in mutually exclusive regions of the genome. To find functions for H3K27me3, we deleted the gene for the putative H3K27 methyltransferase, KMT6, a homolog of Drosophila Enhancer of zeste, E(z). The kmt6 mutant lacks H3K27me3, as shown by western blot and ChIP-seq, displays growth defects, is sterile, and constitutively expresses genes for mycotoxins, pigments and other secondary metabolites. Transcriptome analyses showed that 75% of 4,449 silent genes are enriched for H3K27me3. A subset of genes that were enriched for H3K27me3 in WT gained H3K4me2/3 in kmt6. A largely overlapping set of genes showed increased expression in kmt6. Almost 95% of the remaining 2,720 annotated silent genes showed no enrichment for either H3K27me3 or H3K4me2/3 in kmt6. In these cases mere absence of H3K27me3 was insufficient for expression, which suggests that additional changes are required to activate genes. Taken together, we show that absence of H3K27me3 allowed expression of an additional 14% of the genome, resulting in derepression of genes predominantly involved in secondary metabolite pathways and other species-specific functions, including putative secreted pathogenicity factors. Results from this study provide the framework for novel targeted strategies to control the “cryptic genome”, specifically secondary metabolite expression.

Changes in chromatin structure are required for time- and tissue-specific gene regulation. How exactly these changes are mediated is under intense scrutiny. The interplay between activating histone modifications, e.g. H3K4me, and the silencing H3K27me3 mark has been recognized as critical to orchestrate differentiation and development in plants and animals. Here we show that filamentous fungi, exemplified by the cereal pathogen Fusarium graminearum, can use H3K27 methylation to generate silenced, facultative heterochromatin, covering more than a third of the genome, much more than the 5–8% of Neurospora or metazoan genomes. Removal of the silencing mark by mutation of the methyltransferase subunit of the PRC2 silencing complex resulted in activation of more than 1,500 genes, 14% of the genome. We show that generation of facultative heterochromatin by H3K27 methylation is an ancestral process that has been lost in certain lineages (e.g. at least some hemiascomycetes, the genus Aspergillus and some basidiomycetes). Our studies will open the door to future precise “epigenetic engineering” of gene clusters that generate bioactive compounds, e.g. putative mycotoxins, antibiotics and industrial feedstocks. Availability of tractable fungal model systems for studies of the control and function of H3K27 methylation may accelerate mechanistic research.

Deciphering the cryptic genome: genome-wide analyses of the rice pathogen Fusarium fujikuroi reveal complex regulation of secondary metabolism and novel metabolites

The fungus Fusarium fujikuroi causes "bakanae" disease of rice due to its ability to produce gibberellins (GAs), but it is also known for producing harmful mycotoxins. However, the genetic capacity for the whole arsenal of natural compounds and their role in the fungus' interaction with rice remained unknown. Here, we present a high-quality genome sequence of F. fujikuroi that was assembled into 12 scaffolds corresponding to the 12 chromosomes described for the fungus. We used the genome sequence along with ChIP-seq, transcriptome, proteome, and HPLC-FTMS-based metabolome analyses to identify the potential secondary metabolite biosynthetic gene clusters and to examine their regulation in response to nitrogen availability and plant signals. The results indicate that expression of most but not all gene clusters correlate with proteome and ChIP-seq data. Comparison of the F. fujikuroi genome to those of six other fusaria revealed that only a small number of gene clusters are conserved among these species, thus providing new insights into the divergence of secondary metabolism in the genus Fusarium. Noteworthy, GA biosynthetic genes are present in some related species, but GA biosynthesis is limited to F. fujikuroi, suggesting that this provides a selective advantage during infection of the preferred host plant rice. Among the genome sequences analyzed, one cluster that includes a polyketide synthase gene (PKS19) and another that includes a non-ribosomal peptide synthetase gene (NRPS31) are unique to F. fujikuroi. The metabolites derived from these clusters were identified by HPLC-FTMS-based analyses of engineered F. fujikuroi strains overexpressing cluster genes. In planta expression studies suggest a specific role for the PKS19-derived product during rice infection. Thus, our results indicate that combined comparative genomics and genome-wide experimental analyses identified novel genes and secondary metabolites that contribute to the evolutionary success of F. fujikuroi as a rice pathogen.

Transcriptome characterization by RNA-Seq reveals the involvement of the complement components in noise-traumatized rat cochleae

Acoustic trauma, a leading cause of sensorineural hearing loss in adults, induces a complex degenerative process in the cochlea. Although previous investigations have identified multiple stress pathways, a comprehensive analysis of cochlear responses to acoustic injury is still lacking. In the current study, we used the next-generation RNA-sequencing (RNA-Seq) technique to sequence the whole transcriptome of the normal and noise-traumatized cochlear sensory epithelia (CSE). CSE tissues were collected from rat inner ears 1d after the rats were exposed to a 120-dB (sound pressure level) noise for 2 h. The RNA-Seq generated over 176 million sequence reads for the normal CSE and over 164 million reads for the noise-traumatized CSE. Alignment of these sequences with the rat Rn4 genome revealed the expression of over 17,000 gene transcripts in the CSE, over 2000 of which were exclusively expressed in either the normal or noise-traumatized CSE. Seventy-eight gene transcripts were differentially expressed (70 upregulated and 8 downregulated) after acoustic trauma. Many of the differentially expressed genes are related to the innate immune system. Further expression analyses using quantitative real time PCR confirmed the constitutive expression of multiple complement genes in the normal organ of Corti and the changes in the expression levels of the complement factor I (Cfi) and complement component 1, s subcomponent (C1s) after acoustic trauma. Moreover, protein expression analysis revealed strong expression of Cfi and C1s proteins in the organ of Corti. Importantly, these proteins exhibited expression changes following acoustic trauma. Collectively, the results of the current investigation suggest the involvement of the complement components in cochlear responses to acoustic trauma.

Functional analyses of Trichoderma reesei LAE1 reveal conserved and contrasting roles of this regulator

The putative methyltransferase LaeA is a global regulator that affects the expression of multiple secondary metabolite gene clusters in several fungi, and it can modify heterochromatin structure in Aspergillus nidulans. We have recently shown that the LaeA ortholog of Trichoderma reesei (LAE1), a fungus that is an industrial producer of cellulase and hemicellulase enzymes, regulates the expression of cellulases and polysaccharide hydrolases. To learn more about the function of LAE1 in T. reesei, we assessed the effect of deletion and overexpression of lae1 on genome-wide gene expression. We found that in addition to positively regulating 7 of 17 polyketide or nonribosomal peptide synthases, genes encoding ankyrin-proteins, iron uptake, heterokaryon incompatibility proteins, PTH11-receptors, and oxidases/monoxygenases are major gene categories also regulated by LAE1. chromatin immunoprecipitation sequencing with antibodies against histone modifications known to be associated with transcriptionally active (H3K4me2 and -me3) or silent (H3K9me3) chromatin detected 4089 genes bearing one or more of these methylation marks, of which 75 exhibited a correlation between either H3K4me2 or H3K4me3 and regulation by LAE1. Transformation of a laeA-null mutant of A. nidulans with the T. reesei lae1 gene did not rescue sterigmatocystin formation and further impaired sexual development. LAE1 did not interact with A. nidulans VeA in yeast two-hybrid assays, whereas it interacted with the T. reesei VeA ortholog, VEL1. LAE1 was shown to be required for the expression of vel1, whereas the orthologs of velB and VosA are unaffected by lae1 deletion. Our data show that the biological roles of A. nidulans LaeA and T. reesei LAE1 are much less conserved than hitherto thought. In T. reesei, LAE1 appears predominantly to regulate genes increasing relative fitness in its environment.