Genome-wide transposon tagging reveals location-dependent effects on transcription and chromatin organization in Arabidopsis

The interphase nucleus exists as a highly dynamic system, the physical properties of which have functional importance in gene regulation. Not only can gene expression be influenced by the local sequence context, but also by the architecture of the nucleus in three-dimensions (3D), and by the interactions between these levels via chromatin modifications. A challenging task is to resolve the complex interplay between sequence- and genome structure-based control mechanisms. Here, we created a collection of 277 Arabidopsis lines that allow the visual tracking of individual loci in living plants while comparing gene expression potential at these locations, via an identical reporter cassette. Our studies revealed regional gene silencing near a heterochromatin island, via DNA methylation, that is correlated with mobility constraint and nucleolar association. We also found an example of nucleolar association that does not correlate with gene suppression, suggesting that distinct mechanisms exist that can mediate interactions between chromatin and the nucleolus. These studies demonstrate the utility of this novel resource in unifying structural and functional studies towards a more comprehensive model of how global chromatin organization may coordinate gene expression over large scales.     

Sugarcane genome sequencing by methylation filtration provides tools for genomic research in the genus Saccharum

Many economically important crops have large and complex genomes that hamper their sequencing by standard methods such as whole genome shotgun (WGS). Large tracts of methylated repeats occur in plant genomes that are interspersed by hypomethylated gene‐rich regions. Gene‐enrichment strategies based on methylation profiles offer an alternative to sequencing repetitive genomes. Here, we have applied methyl filtration with McrBC endonuclease digestion to enrich for euchromatic regions in the sugarcane genome. To verify the efficiency of methylation filtration and the assembly quality of sequences submitted to gene‐enrichment strategy, we have compared assemblies using methyl‐filtered (MF) and unfiltered (UF) libraries. The use of methy filtration allowed a better assembly by filtering out 35% of the sugarcane genome and by producing 1.5× more scaffolds and 1.7× more assembled Mb in length compared with unfiltered dataset. The coverage of sorghum coding sequences (CDS) by MF scaffolds was at least 36% higher than by the use of UF scaffolds. Using MF technology, we increased by 134× the coverage of gene regions of the monoploid sugarcane genome. The MF reads assembled into scaffolds that covered all genes of the sugarcane bacterial artificial chromosomes (BACs), 97.2% of sugarcane expressed sequence tags (ESTs), 92.7% of sugarcane RNA‐seq reads and 98.4% of sorghum protein sequences. Analysis of MF scaffolds from encoded enzymes of the sucrose/starch pathway discovered 291 single‐nucleotide polymorphisms (SNPs) in the wild sugarcane species, S. spontaneum and S. officinarum. A large number of microRNA genes was also identified in the MF scaffolds. The information achieved by the MF dataset provides a valuable tool for genomic research in the genus Saccharum and for improvement of sugarcane as a biofuel crop.     

Validation and assessment of variant calling pipelines for next-generation sequencing

Background

The processing and analysis of the large scale data generated by next-generation sequencing (NGS) experiments is challenging and is a burgeoning area of new methods development. Several new bioinformatics tools have been developed for calling sequence variants from NGS data. Here, we validate the variant calling of these tools and compare their relative accuracy to determine which data processing pipeline is optimal.

Results

We developed a unified pipeline for processing NGS data that encompasses four modules: mapping, filtering, realignment and recalibration, and variant calling. We processed 130 subjects from an ongoing whole exome sequencing study through this pipeline. To evaluate the accuracy of each module, we conducted a series of comparisons between the single nucleotide variant (SNV) calls from the NGS data and either gold-standard Sanger sequencing on a total of 700 variants or array genotyping data on a total of 9,935 single-nucleotide polymorphisms. A head to head comparison showed that Genome Analysis Toolkit (GATK) provided more accurate calls than SAMtools (positive predictive value of 92.55% vs. 80.35%, respectively). Realignment of mapped reads and recalibration of base quality scores before SNV calling proved to be crucial to accurate variant calling. GATK HaplotypeCaller algorithm for variant calling outperformed the UnifiedGenotype algorithm. We also showed a relationship between mapping quality, read depth and allele balance, and SNV call accuracy. However, if best practices are used in data processing, then additional filtering based on these metrics provides little gains and accuracies of >99% are achievable.

Conclusions

Our findings will help to determine the best approach for processing NGS data to confidently call variants for downstream analyses. To enable others to implement and replicate our results, all of our codes are freely available at http://metamoodics.org/wes.

     

Transpositional landscape of the rice genome revealed by paired-end mapping of high-throughput re-sequencing data

Transposable elements (TEs) are mobile entities that densely populate most eukaryotic genomes and contribute to both their structural and functional dynamics. However, most TE-related sequences in both plant and animal genomes correspond to inactive, degenerated elements, due to the combined effect of silencing pathways and elimination through deletions. One of the major difficulties in fully characterizing the molecular basis of genetic diversity of a given species lies in establishing its genome-wide transpositional activity. Here, we provide an extensive survey of the transpositional landscape of a plant genome using a deep sequencing strategy. This was achieved through paired-end mapping of a fourfold coverage of the genome of rice mutant line derived from an in vitro callus culture using Illumina technology. Our study shows that at least 13 TE families are active in this genotype, causing 34 new insertions. This next-generation sequencing-based strategy provides new opportunities to quantify the impact of TEs on the genome dynamics of the species.     

Draft genome sequence of pigeonpea (Cajanus cajan), an orphan legume crop of resource-poor farmers

Pigeonpea is an important legume food crop grown primarily by smallholder farmers in many semi-arid tropical regions of the world. We used the Illumina next-generation sequencing platform to generate 237.2 Gb of sequence, which along with Sanger-based bacterial artificial chromosome end sequences and a genetic map, we assembled into scaffolds representing 72.7% (605.78 Mb) of the 833.07 Mb pigeonpea genome. Genome analysis predicted 48,680 genes for pigeonpea and also showed the potential role that certain gene families, for example, drought tolerance-related genes, have played throughout the domestication of pigeonpea and the evolution of its ancestors. Although we found a few segmental duplication events, we did not observe the recent genome-wide duplication events observed in soybean. This reference genome sequence will facilitate the identification of the genetic basis of agronomically important traits, and accelerate the development of improved pigeonpea varieties that could improve food security in many developing countries.     

Expressed sequence tag analysis in <Emphasis Type="Italic">Cycas</Emphasis>, the most primitive living seed plant

Background

Cycads are ancient seed plants (living fossils) with origins in the Paleozoic. Cycads are sometimes considered a 'missing link' as they exhibit characteristics intermediate between vascular non-seed plants and the more derived seed plants. Cycads have also been implicated as the source of 'Guam's dementia', possibly due to the production of S(+)-beta-methyl-alpha, beta-diaminopropionic acid (BMAA), which is an agonist of animal glutamate receptors.

Results

A total of 4,200 expressed sequence tags (ESTs) were created from Cycas rumphii and clustered into 2,458 contigs, of which 1,764 had low-stringency BLAST similarity to other plant genes. Among those cycad contigs with similarity to plant genes, 1,718 cycad 'hits' are to angiosperms, 1,310 match genes in gymnosperms and 734 match lower (non-seed) plants. Forty-six contigs were found that matched only genes in lower plants and gymnosperms. Upon obtaining the complete sequence from the clones of 37/46 contigs, 14 still matched only gymnosperms. Among those cycad contigs common to higher plants, ESTs were discovered that correspond to those involved in development and signaling in present-day flowering plants. We purified a cycad EST for a glutamate receptor (GLR)-like gene, as well as ESTs potentially involved in the synthesis of the GLR agonist BMAA.

Conclusions

Analysis of cycad ESTs has uncovered conserved and potentially novel genes. Furthermore, the presence of a glutamate receptor agonist, as well as a glutamate receptor-like gene in cycads, supports the hypothesis that such neuroactive plant products are not merely herbivore deterrents but may also serve a role in plant signaling.

     

Identification of mesoderm development (mesd) candidate genes by comparative mapping and genome sequence analysis

The proximal albino deletions identify several functional regions on mouse Chromosome 7 critical for differentiation of mesoderm (mesd), development of the hypothalamus neuroendocrine lineage (nelg), and function of the liver (hsdr1). Using comparative mapping and genomic sequence analysis, we have identified four novel genes and Il16 in the mesd deletion interval. Two of the novel genes, mesdc1 and mesdc2, are located within the mesd critical region defined by BAC transgenic rescue. We have investigated the fetal role of genes located outside the mesd critical region using BAC transgenic complementation of the mesd early embryonic lethality. Using human radiation hybrid mapping and BAC contig construction, we have identified a conserved region of human chromosome 15 homologous to the mesd, nelg, and hsdr1 functional regions. Three human diseases cosegregate with microsatellite markers used in construction of the human BAC/YAC physical map, including autosomal dominant nocturnal frontal lobe epilepsy (ENFL2; also known as ADNFLE), a syndrome of mental retardation, spasticity, and tapetoretinal degeneration (MRST); and a pyogenic arthritis, pyoderma gangrenosum, and acne syndrome (PAPA).     

A Complementary Method for Production of Tetraploid Crassostrea gigas Using Crosses Between Diploids and Tetraploids with Cytochalasin B Treatments

We present a new method to produce tetraploid Crassostrea gigas by cytochalasin B inhibition of polar body 2 expulsion in diploid females crossed with tetraploid males. This offers a means of direct introgression of genetic characters from selected diploid to tetraploid lines, avoiding a triploid step. Offspring larval ploidy shifted over time and depended on size, with tetraploids more frequent among the smaller larvae and triploids among the large. Viable tetraploids were found at 4 and 6 months, indicating the technique was successful. The possibility that gynogenesis occurred was tested by microsatellite analysis to confirm the presence of paternally inherited alleles. These were present in all animals of the 2n × 4n + CB (female first) cross. However, a 4n × 2n + CB cross produced triploids, including some gynogens. Our method illustrates for the first time that diploid C. gigas eggs, if selected for large size, can give viable tetraploid offspring.     

Secretion trap tagging of secreted and membrane-spanning proteins using Arabidopsis gene traps

Secreted and membrane-spanning proteins play fundamental roles in plant development but pose challenges for genetic identification and characterization. We describe a "secretion trap" screen for gene trap insertions in genes encoding proteins routed through the secretory pathway. The gene trap transposon encodes a beta-glucuronidase reporter enzyme that is inhibited by N-linked glycosylation specific to the secretory pathway. Treatment of seedlings with tunicamycin inhibits glycosylation, resulting in increased activity of secreted beta-glucuronidase fusions that result from gene trap integration downstream of exons encoding signal peptides. In the 2,059 gene trap lines that we screened, 32 secretion trap expression patterns were identified in a wide variety of tissues including embryos, meristems, and the developing vasculature. Genes disrupted by the secretion traps encode putative extracellular signaling proteins, membrane transport proteins, and novel secreted proteins of unknown function missed by conventional mutagenesis and gene prediction. Secretion traps provide a unique reagent for gene expression studies and can guide the genetic combination of loss of function alleles in related genes.