KHSRP loss increases neuronal growth and synaptic transmission and alters memory consolidation through RNA stabilization.

The KH-type splicing regulatory protein (KHSRP) is an RNA-binding protein linked to decay of mRNAs with AU-rich elements. KHSRP was previously shown to destabilize Gap43 mRNA and decrease neurite growth in cultured embryonic neurons. Here, we have tested functions of KHSRP in vivo. We find upregulation of 1460 mRNAs in neocortex of adult Khsrp-/- mice, of which 527 bind to KHSRP with high specificity. These KHSRP targets are involved in pathways for neuronal morphology, axon guidance, neurotransmission and long-term memory. Khsrp-/- mice show increased axon growth and dendritic spine density in vivo. Neuronal cultures from Khsrp-/- mice show increased axon and dendrite growth and elevated KHSRP-target mRNAs, including subcellularly localized mRNAs. Furthermore, neuron-specific knockout of Khsrp confirms these are from neuron-intrinsic roles of KHSRP. Consistent with this, neurons in the hippocampus and infralimbic cortex of Khsrp-/- mice show elevations in frequency of miniature excitatory postsynaptic currents. The Khsrp-/- mice have deficits in trace conditioning and attention set-shifting tasks compared Khsrp+/+ mice, indicating impaired prefrontal- and hippocampal-dependent memory consolidation with loss of KHSRP. Overall, these results indicate that deletion of KHSRP impairs neuronal development resulting in alterations in neuronal morphology and function by changing post-transcriptional control of neuronal gene expression.

High-Quality Reference Genome Sequence for the Oomycete Vegetable Pathogen Phytophthora capsici Strain LT1534.

The oomycete Phytophthora capsici is a destructive pathogen of a wide range of vegetable hosts, especially peppers and cucurbits. A 94.17-Mb genome assembly was constructed using PacBio and Illumina data and annotated with support from transcriptome sequencing (RNA-Seq) reads.

Corrigendum: Novel Meiotic miRNAs and Indications for a Role of PhasiRNAs in Meiosis.

[This corrects the article DOI: 10.3389/fpls.2016.00762.].

Dynamic Extreme Aneuploidy (DEA) in the vegetable pathogen Phytophthora capsici and the potential for rapid asexual evolution.

Oomycete plant pathogens are difficult to control and routine genetic research is challenging. A major problem is instability of isolates. Here we characterize >600 field and single zoospore isolates of Phytophthora capsici for inheritance of mating type, sensitivity to mefenoxam, chromosome copy number and heterozygous allele frequencies. The A2 mating type was highly unstable with 26% of 241 A2 isolates remaining A2. The A1 mating type was stable. Isolates intermediately resistant to mefenoxam produced fully resistant single-spore progeny. Sensitive isolates remained fully sensitive. Genome re-sequencing of single zoospore isolates revealed extreme aneuploidy; a phenomenon dubbed Dynamic Extreme Aneuploidy (DEA). DEA is characterized by the asexual inheritance of diverse intra-genomic combinations of chromosomal ploidy ranging from 2N to 3N and heterozygous allele frequencies that do not strictly correspond to ploidy. Isolates sectoring on agar media showed dramatically altered heterozygous allele frequencies. DEA can explain the rapid increase of advantageous alleles (e.g. drug resistance), mating type switches and copy neutral loss of heterozygosity (LOH). Although the mechanisms driving DEA are unknown, it can play an important role in adaptation and evolution and seriously hinders all aspects of P. capsici research.

Experimental drought reduces genetic diversity in the grassland foundation species Bouteloua eriopoda.

Understanding the resistance and resilience of foundation plant species to climate change is a critical issue because the loss of these species would fundamentally reshape communities and ecosystem processes. High levels of population genetic diversity may buffer foundation species against climate disruptions, but the strong selective pressures associated with climatic shifts may also rapidly reduce such diversity. We characterized genetic diversity and its responsiveness to experimental drought in the foundation plant, black grama grass (Bouteloua eriopoda), which dominates many western North American grasslands and shrublands. Previous studies suggested that in arid ecosystems, black grama reproduces largely asexually via stolons, and thus is likely to have low genetic variability, which might limit its potential to respond to climate disruptions. Using genotyping-by-sequencing, we demonstrated unexpectedly high genetic variability among black grama plants in a 1 ha site within the Sevilleta National Wildlife Refuge in central New Mexico, suggesting some level of sexual reproduction. Three years of experimental, growing season drought reduced black grama survival and biomass (the latter by 96%), with clear genetic differentiation (higher FST) between plants succumbing to drought and those remaining alive. Reduced genetic variability in the surviving plants in drought plots indicated that the experimental drought had forced black grama populations through selection bottlenecks. These results suggest that foundation grass species, such as black grama, may experience rapid evolutionary change if future climates include more severe droughts.

Exploring Frequented Regions in Pan-Genomic Graphs.

We consider the problem of identifying regions within a pan-genome De Bruijn graph that are traversed by many sequence paths. We define such regions and the subpaths that traverse them as frequented regions (FRs). In this work, we formalize the FR problem and describe an efficient algorithm for finding FRs. Subsequently, we propose some applications of FRs based on machine-learning and pan-genome graph simplification. We demonstrate the effectiveness of these applications using data sets for the organisms Staphylococcus aureus (bacterium) and Saccharomyces cerevisiae (yeast). We corroborate the biological relevance of FRs such as identifying introgressions in yeast that aid in alcohol tolerance, and show that FRs are useful for classification of yeast strains by industrial use and visualizing pan-genomic space.

Genome and evolution of the arbuscular mycorrhizal fungus Diversispora epigaea (formerly Glomus versiforme) and its bacterial endosymbionts.

Arbuscular mycorrhizal (AM) fungi form endosymbioses with most plants, and they themselves are hosts for Mollicutes/Mycoplasma-related endobacteria (MRE). Despite their significance, genomic information for AM fungi and their MRE are relatively sparse, which hinders our understanding of their biology and evolution. We assembled the genomes of the AM fungus Diversispora epigaea (formerly Glomus versiforme) and its MRE and performed comparative genomics and evolutionary analyses. The D. epigaea genome showed a pattern of substantial gene duplication and differential evolution of gene families, including glycosyltransferase family 25, whose activities are exclusively lipopolysaccharide biosynthesis. Genes acquired by horizontal transfer from bacteria possibly function in defense against foreign DNA or viruses. The MRE population was diverse, with multiple genomes displaying characteristics of differential evolution and encoding many MRE-specific genes as well as genes of AM fungal origin. Gene family expansion in D. epigaea may enhance adaptation to both external and internal environments, such as expansion of kinases for signal transduction upon external stimuli and expansion of nucleoside salvage pathway genes potentially for competition with MRE, whose genomes lack purine and pyrimidine biosynthetic pathways. Collectively, this metagenome provides high-quality references and begins to reveal the diversity within AM fungi and their MRE.

Histone Citrullination Represses MicroRNA Expression, Resulting in Increased Oncogene mRNAs in Somatolactotrope Cells.

Peptidylarginine deiminase (PAD) enzymes convert histone arginine residues into citrulline to modulate chromatin organization and gene expression. Although PADs are expressed in anterior pituitary gland cells, their functional role and expression in pituitary adenomas are unknown. To begin to address these issues, we first examined normal human pituitaries and pituitary adenomas and found that PAD2, PAD4, and citrullinated histones are highest in prolactinomas and somatoprolactinomas. In the somatoprolactinoma-derived GH3 cell line, PADs citrullinate histone H3, which is attenuated by a pan-PAD inhibitor. RNA sequencing and chromatin immunoprecipitation (ChIP) studies show that the expression of microRNAs (miRNAs) let-7c-2, 23b, and 29c is suppressed by histone citrullination. Our studies demonstrate that these miRNAs directly target the mRNA of the oncogenes encoding HMGA, insulin-like growth factor 1 (IGF-1), and N-MYC, which are highly implicated in human prolactinoma/somatoprolactinoma pathogenesis. Our results are the first to define a direct role for PAD-catalyzed histone citrullination in miRNA expression, which may underlie the etiology of prolactinoma and somatoprolactinoma tumors through regulation of oncogene expression.

High-resolution crossover mapping reveals similarities and differences of male and female recombination in maize.

Meiotic crossovers (COs) are not uniformly distributed across the genome. Factors affecting this phenomenon are not well understood. Although many species exhibit large differences in CO numbers between sexes, sex-specific aspects of CO landscape are particularly poorly elucidated. Here, we conduct high-resolution CO mapping in maize. Our results show that CO numbers as well as their overall distribution are similar in male and female meioses. There are, nevertheless, dissimilarities at local scale. Male and female COs differ in their locations relative to transcription start sites in gene promoters and chromatin marks, including nucleosome occupancy and tri-methylation of lysine 4 of histone H3 (H3K4me3). Our data suggest that sex-specific factors not only affect male-female CO number disparities but also cause fine differences in CO positions. Differences between male and female CO landscapes indicate that recombination has distinct implications for population structure and gene evolution in male and in female meioses.

Correction to: The transcriptome landscape of early maize meiosis.

CORRECTION: Following publication of the original article [1], the authors reported that the number of genes overlaying the bar graph in Fig. 3A were incorrectly counted and inserted (i.e. including a title tile, and in inverse order). The corrected numbers are below and match with the listed genes supplied in Additional File: Table S2.

Comparative Genomics of an Unusual Biogeographic Disjunction in the Cotton Tribe (Gossypieae) Yields Insights into Genome Downsizing.

Long-distance insular dispersal is associated with divergence and speciation because of founder effects and strong genetic drift. The cotton tribe (Gossypieae) has experienced multiple transoceanic dispersals, generating an aggregate geographic range that encompasses much of the tropics and subtropics worldwide. Two genera in the Gossypieae, Kokia and Gossypioides, exhibit a remarkable geographic disjunction, being restricted to the Hawaiian Islands and Madagascar/East Africa, respectively. We assembled and use de novo genome sequences to address questions regarding the divergence of these two genera from each other and from their sister-group, Gossypium. In addition, we explore processes underlying the genome downsizing that characterizes Kokia and Gossypioides relative to other genera in the tribe. Using 13,000 gene orthologs and synonymous substitution rates, we show that the two disjuncts last shared a common ancestor ∼5 Ma, or half as long ago as their divergence from Gossypium. We report relative stasis in the transposable element fraction. In comparison to Gossypium, there is loss of ∼30% of the gene content in the two disjunct genera and a history of genome-wide accumulation of deletions. In both genera, there is a genome-wide bias toward deletions over insertions, and the number of gene losses exceeds the number of gains by ∼2- to 4-fold. The genomic analyses presented here elucidate genomic consequences of the demographic and biogeographic history of these closest relatives of Gossypium, and enhance their value as phylogenetic outgroups.

Strategies for optimizing BioNano and Dovetail explored through a second reference quality assembly for the legume model, Medicago truncatula.

BACKGROUND: Third generation sequencing technologies, with sequencing reads in the tens- of kilo-bases, facilitate genome assembly by spanning ambiguous regions and improving continuity. This has been critical for plant genomes, which are difficult to assemble due to high repeat content, gene family expansions, segmental and tandem duplications, and polyploidy. Recently, high-throughput mapping and scaffolding strategies have further improved continuity. Together, these long-range technologies enable quality draft assemblies of complex genomes in a cost-effective and timely manner. RESULTS: Here, we present high quality genome assemblies of the model legume plant, Medicago truncatula (R108) using PacBio, Dovetail Chicago (hereafter, Dovetail) and BioNano technologies. To test these technologies for plant genome assembly, we generated five assemblies using all possible combinations and ordering of these three technologies in the R108 assembly. While the BioNano and Dovetail joins overlapped, they also showed complementary gains in continuity and join numbers. Both technologies spanned repetitive regions that PacBio alone was unable to bridge. Combining technologies, particularly Dovetail followed by BioNano, resulted in notable improvements compared to Dovetail or BioNano alone. A combination of PacBio, Dovetail, and BioNano was used to generate a high quality draft assembly of R108, a M. truncatula accession widely used in studies of functional genomics. As a test for the usefulness of the resulting genome sequence, the new R108 assembly was used to pinpoint breakpoints and characterize flanking sequence of a previously identified translocation between chromosomes 4 and 8, identifying more than 22.7 Mb of novel sequence not present in the earlier A17 reference assembly. CONCLUSIONS: Adding Dovetail followed by BioNano data yielded complementary improvements in continuity over the original PacBio assembly. This strategy proved efficient and cost-effective for developing a quality draft assembly compared to traditional reference assemblies.

Hybrid assembly with long and short reads improves discovery of gene family expansions.

BACKGROUND: Long-read and short-read sequencing technologies offer competing advantages for eukaryotic genome sequencing projects. Combinations of both may be appropriate for surveys of within-species genomic variation. METHODS: We developed a hybrid assembly pipeline called "Alpaca" that can operate on 20X long-read coverage plus about 50X short-insert and 50X long-insert short-read coverage. To preclude collapse of tandem repeats, Alpaca relies on base-call-corrected long reads for contig formation. RESULTS: Compared to two other assembly protocols, Alpaca demonstrated the most reference agreement and repeat capture on the rice genome. On three accessions of the model legume Medicago truncatula, Alpaca generated the most agreement to a conspecific reference and predicted tandemly repeated genes absent from the other assemblies. CONCLUSION: Our results suggest Alpaca is a useful tool for investigating structural and copy number variation within de novo assemblies of sampled populations.

Exploring structural variation and gene family architecture with De Novo assemblies of 15 Medicago genomes.

BACKGROUND: Previous studies exploring sequence variation in the model legume, Medicago truncatula, relied on mapping short reads to a single reference. However, read-mapping approaches are inadequate to examine large, diverse gene families or to probe variation in repeat-rich or highly divergent genome regions. De novo sequencing and assembly of M. truncatula genomes enables near-comprehensive discovery of structural variants (SVs), analysis of rapidly evolving gene families, and ultimately, construction of a pan-genome. RESULTS: Genome-wide synteny based on 15 de novo M. truncatula assemblies effectively detected different types of SVs indicating that as much as 22% of the genome is involved in large structural changes, altogether affecting 28% of gene models. A total of 63 million base pairs (Mbp) of novel sequence was discovered, expanding the reference genome space for Medicago by 16%. Pan-genome analysis revealed that 42% (180 Mbp) of genomic sequences is missing in one or more accession, while examination of de novo annotated genes identified 67% (50,700) of all ortholog groups as dispensable - estimates comparable to recent studies in rice, maize and soybean. Rapidly evolving gene families typically associated with biotic interactions and stress response were found to be enriched in the accession-specific gene pool. The nucleotide-binding site leucine-rich repeat (NBS-LRR) family, in particular, harbors the highest level of nucleotide diversity, large effect single nucleotide change, protein diversity, and presence/absence variation. However, the leucine-rich repeat (LRR) and heat shock gene families are disproportionately affected by large effect single nucleotide changes and even higher levels of copy number variation. CONCLUSIONS: Analysis of multiple M. truncatula genomes illustrates the value of de novo assemblies to discover and describe structural variation, something that is often under-estimated when using read-mapping approaches. Comparisons among the de novo assemblies also indicate that different large gene families differ in the architecture of their structural variation.

Gene Evolutionary Trajectories and GC Patterns Driven by Recombination in Zea mays.

Recombination occurring during meiosis is critical for creating genetic variation and plays an essential role in plant evolution. In addition to creating novel gene combinations, recombination can affect genome structure through altering GC patterns. In maize (Zea mays) and other grasses, another intriguing GC pattern exists. Maize genes show a bimodal GC content distribution that has been attributed to nucleotide bias in the third, or wobble, position of the codon. Recombination may be an underlying driving force given that recombination sites are often associated with high GC content. Here we explore the relationship between recombination and genomic GC patterns by comparing GC gene content at each of the three codon positions (GC1, GC2, and GC3, collectively termed GCx) to instances of a variable GC-rich motif that underlies double strand break (DSB) hotspots and to meiocyte-specific gene expression. Surprisingly, GCx bimodality in maize cannot be fully explained by the codon wobble hypothesis. High GCx genes show a strong overlap with the DSB hotspot motif, possibly providing a mechanism for the high evolutionary rates seen in these genes. On the other hand, genes that are turned on in meiosis (early prophase I) are biased against both high GCx genes and genes with the DSB hotspot motif, possibly allowing important meiotic genes to avoid DSBs. Our data suggests a strong link between the GC-rich motif underlying DSB hotspots and high GCx genes.

Transcriptome Sequencing of Diverse Peanut (Arachis) Wild Species and the Cultivated Species Reveals a Wealth of Untapped Genetic Variability.

To test the hypothesis that the cultivated peanut species possesses almost no molecular variability, we sequenced a diverse panel of 22 Arachis accessions representing Arachis hypogaea botanical classes, A-, B-, and K- genome diploids, a synthetic amphidiploid, and a tetraploid wild species. RNASeq was performed on pools of three tissues, and de novo assembly was performed. Realignment of individual accession reads to transcripts of the cultivar OLin identified 306,820 biallelic SNPs. Among 10 naturally occurring tetraploid accessions, 40,382 unique homozygous SNPs were identified in 14,719 contigs. In eight diploid accessions, 291,115 unique SNPs were identified in 26,320 contigs. The average SNP rate among the 10 cultivated tetraploids was 0.5, and among eight diploids was 9.2 per 1000 bp. Diversity analysis indicated grouping of diploids according to genome classification, and cultivated tetraploids by subspecies. Cluster analysis of variants indicated that sequences of B genome species were the most similar to the tetraploids, and the next closest diploid accession belonged to the A genome species. A subset of 66 SNPs selected from the dataset was validated; of 782 SNP calls, 636 (81.32%) were confirmed using an allele-specific discrimination assay. We conclude that substantial genetic variability exists among wild species. Additionally, significant but lesser variability at the molecular level occurs among accessions of the cultivated species. This survey is the first to report significant SNP level diversity among transcripts, and may explain some of the phenotypic differences observed in germplasm surveys. Understanding SNP variants in the Arachis accessions will benefit in developing markers for selection.

Augmenting Chinese hamster genome assembly by identifying regions of high confidence.

Chinese hamster Ovary (CHO) cell lines are the dominant industrial workhorses for therapeutic recombinant protein production. The availability of genome sequence of Chinese hamster and CHO cells will spur further genome and RNA sequencing of producing cell lines. However, the mammalian genomes assembled using shot-gun sequencing data still contain regions of uncertain quality due to assembly errors. Identifying high confidence regions in the assembled genome will facilitate its use for cell engineering and genome engineering. We assembled two independent drafts of Chinese hamster genome by de novo assembly from shotgun sequencing reads and by re-scaffolding and gap-filling the draft genome from NCBI for improved scaffold lengths and gap fractions. We then used the two independent assemblies to identify high confidence regions using two different approaches. First, the two independent assemblies were compared at the sequence level to identify their consensus regions as "high confidence regions" which accounts for at least 78 % of the assembled genome. Further, a genome wide comparison of the Chinese hamster scaffolds with mouse chromosomes revealed scaffolds with large blocks of collinearity, which were also compiled as high-quality scaffolds. Genome scale collinearity was complemented with EST based synteny which also revealed conserved gene order compared to mouse. As cell line sequencing becomes more commonly practiced, the approaches reported here are useful for assessing the quality of assembly and potentially facilitate the engineering of cell lines.

Novel Meiotic miRNAs and Indications for a Role of PhasiRNAs in Meiosis.

Small RNAs (sRNA) add additional layers to the regulation of gene expression, with siRNAs directing gene silencing at the DNA level by RdDM (RNA-directed DNA methylation), and micro RNAs (miRNAs) directing post-transcriptional regulation of specific target genes, mostly by mRNA cleavage. We used manually isolated male meiocytes from maize (Zea mays) to investigate sRNA and DNA methylation landscapes during zygotene, an early stage of meiosis during which steps of meiotic recombination and synapsis of paired homologous chromosomes take place. We discovered two novel miRNAs from meiocytes, zma-MIR11969 and zma-MIR11970, and identified putative target genes. Furthermore, we detected abundant phasiRNAs of 21 and 24 nt length. PhasiRNAs are phased small RNAs which occur in 21 or 24 nt intervals, at a few hundred loci, specifically in male reproductive tissues in grasses. So far, the function of phasiRNAs remained elusive. Data from isolated meiocytes now revealed elevated DNA methylation at phasiRNA loci, especially in the CHH context, suggesting a role for phasiRNAs in cis DNA methylation. In addition, we consider a role of these phasiRNAs in chromatin remodeling/dynamics during meiosis. However, this is not well supported yet and will need more additional data. Here, we only lay out the idea due to other relevant literature and our additional observation of a peculiar GC content pattern at phasiRNA loci. Chromatin remodeling is also indicated by the discovery that histone genes were enriched for sRNA of 22 nt length. Taken together, we gained clues that lead us to hypothesize sRNA-driven DNA methylation and possibly chromatin remodeling during male meiosis in the monocot maize which is in line with and extends previous knowledge.

Characterization, correction and de novo assembly of an Oxford Nanopore genomic dataset from Agrobacterium tumefaciens.

The MinION is a portable single-molecule DNA sequencing instrument that was released by Oxford Nanopore Technologies in 2014, producing long sequencing reads by measuring changes in ionic flow when single-stranded DNA molecules translocate through the pores. While MinION long reads have an error rate substantially higher than the ones produced by short-read sequencing technologies, they can generate de novo assemblies of microbial genomes, after an initial correction step that includes alignment of Illumina sequencing data or detection of overlaps between Oxford Nanopore reads to improve accuracy. In this study, MinION reads were generated from the multi-chromosome genome of Agrobacterium tumefaciens strain LBA4404. Errors in the consensus two-directional (sense and antisense) "2D" sequences were first characterized by way of comparison with an internal reference assembly. Both Illumina-based correction and self-correction were performed and the resulting corrected reads assembled into high-quality hybrid and non-hybrid assemblies. Corrected read datasets and assemblies were subsequently compared. The results shown here indicate that both hybrid and non-hybrid methods can be used to assemble Oxford Nanopore reads into informative multi-chromosome assemblies, each with slightly different outcomes in terms of contiguity and accuracy.

Ancient horizontal transfers of retrotransposons between birds and ancestors of human pathogenic nematodes.

Parasite host switches may trigger disease emergence, but prehistoric host ranges are often unknowable. Lymphatic filariasis and loiasis are major human diseases caused by the insect-borne filarial nematodes Brugia, Wuchereria and Loa. Here we show that the genomes of these nematodes and seven tropical bird lineages exclusively share a novel retrotransposon, AviRTE, resulting from horizontal transfer (HT). AviRTE subfamilies exhibit 83-99% nucleotide identity between genomes, and their phylogenetic distribution, paleobiogeography and invasion times suggest that HTs involved filarial nematodes. The HTs between bird and nematode genomes took place in two pantropical waves, >25-22 million years ago (Myr ago) involving the Brugia/Wuchereria lineage and >20-17 Myr ago involving the Loa lineage. Contrary to the expectation from the mammal-dominated host range of filarial nematodes, we hypothesize that these major human pathogens may have independently evolved from bird endoparasites that formerly infected the global breadth of avian biodiversity.