Using multiple reference genomes to identify and resolve annotation inconsistencies

BACKGROUND

Advances in sequencing technologies have led to the release of reference genomes and annotations for multiple individuals within more well-studied systems. While each of these new genome assemblies shares significant portions of synteny between each other, the annotated structure of gene models within these regions can differ. Of particular concern are split-gene misannotations, in which a single gene is incorrectly annotated as two distinct genes or two genes are incorrectly annotated as a single gene. These misannotations can have major impacts on functional prediction, estimates of expression, and many downstream analyses.

RESULTS

We developed a high-throughput method based on pairwise comparisons of annotations that detect potential split-gene misannotations and quantifies support for whether the genes should be merged into a single gene model. We demonstrated the utility of our method using gene annotations of three reference genomes from maize (B73, PH207, and W22), a difficult system from an annotation perspective due to the size and complexity of the genome. On average, we found several hundred of these potential split-gene misannotations in each pairwise comparison, corresponding to 3-5% of gene models across annotations. To determine which state (i.e. one gene or multiple genes) is biologically supported, we utilized RNAseq data from 10 tissues throughout development along with a novel metric and simulation framework. The methods we have developed require minimal human interaction and can be applied to future assemblies to aid in annotation efforts.

CONCLUSIONS

Split-gene misannotations occur at appreciable frequency in maize annotations. We have developed a method to easily identify and correct these misannotations. Importantly, this method is generic in that it can utilize any type of short-read expression data. Failure to account for split-gene misannotations has serious consequences for biological inference, particularly for expression-based analyses.

Molecular characterization of the insecticidal activity of double-stranded RNA targeting the smooth septate junction of western corn rootworm (Diabrotica virgifera virgifera)

The western corn rootworm (WCR, Diabrotica virgifera virgifera) gene, dvssj1, is a putative homolog of the Drosophila melanogaster gene, snakeskin (ssk). This gene encodes a membrane protein associated with the smooth septate junction (SSJ) which is required for the proper barrier function of the epithelial lining of insect intestines. Disruption of DVSSJ integrity by RNAi technique has been shown previously to be an effective approach for corn rootworm control, by apparent suppression of production of DVSSJ1 protein leading to growth inhibition and mortality. To understand the mechanism that leads to the death of WCR larvae by dvssj1 double-stranded RNA, we examined the molecular characteristics associated with SSJ functions during larval development. Dvssj1 dsRNA diet feeding results in dose-dependent suppression of mRNA and protein; this impairs SSJ formation and barrier function of the midgut and results in larval mortality. These findings suggest that the malfunctioning of the SSJ complex in midgut triggered by dvssj1 silencing is the principal cause of WCR death. This study also illustrates that dvssj1 is a midgut-specific gene in WCR and its functions are consistent with biological functions described for ssk.

Shared and genetically distinct Zea mays transcriptome responses to ongoing and past low temperature exposure

BACKGROUND

Cold temperatures and their alleviation affect many plant traits including the abundance of protein coding gene transcripts. Transcript level changes that occur in response to cold temperatures and their alleviation are shared or vary across genotypes. In this study we identify individual transcripts and groups of functionally related transcripts that consistently respond to cold and its alleviation. Genes that respond differently to temperature changes across genotypes may have limited functional importance. We investigate if these genes share functions, and if their genotype-specific gene expression levels change in magnitude or rank across temperatures.

RESULTS

We estimate transcript abundances from over 22,000 genes in two unrelated Zea mays inbred lines during and after cold temperature exposure. Genotype and temperature contribute to many genes' abundances. Past cold exposure affects many fewer genes. Genes up-regulated in cold encode many cytokinin glucoside biosynthesis enzymes, transcription factors, signalling molecules, and proteins involved in diverse environmental responses. After cold exposure, protease inhibitors and cuticular wax genes are newly up-regulated, and environmentally responsive genes continue to be up-regulated. Genes down-regulated in response to cold include many photosynthesis, translation, and DNA replication associated genes. After cold exposure, DNA replication and translation genes are still preferentially downregulated. Lignin and suberin biosynthesis are newly down-regulated. DNA replication, reactive oxygen species response, and anthocyanin biosynthesis genes have strong, genotype-specific temperature responses. The ranks of genotypes' transcript abundances often change across temperatures.

CONCLUSIONS

We report a large, core transcriptome response to cold and the alleviation of cold. In cold, many of the core suite of genes are up or downregulated to control plant growth and photosynthesis and limit cellular damage. In recovery, core responses are in part to prepare for future stress. Functionally related genes are consistently and greatly up-regulated in a single genotype in response to cold or its alleviation, suggesting positive selection has driven genotype-specific temperature responses in maize.

Parallel altitudinal clines reveal trends in adaptive evolution of genome size in Zea mays

While the vast majority of genome size variation in plants is due to differences in repetitive sequence, we know little about how selection acts on repeat content in natural populations. Here we investigate parallel changes in intraspecific genome size and repeat content of domesticated maize (Zea mays) landraces and their wild relative teosinte across altitudinal gradients in Mesoamerica and South America. We combine genotyping, low coverage whole-genome sequence data, and flow cytometry to test for evidence of selection on genome size and individual repeat abundance. We find that population structure alone cannot explain the observed variation, implying that clinal patterns of genome size are maintained by natural selection. Our modeling additionally provides evidence of selection on individual heterochromatic knob repeats, likely due to their large individual contribution to genome size. To better understand the phenotypes driving selection on genome size, we conducted a growth chamber experiment using a population of highland teosinte exhibiting extensive variation in genome size. We find weak support for a positive correlation between genome size and cell size, but stronger support for a negative correlation between genome size and the rate of cell production. Reanalyzing published data of cell counts in maize shoot apical meristems, we then identify a negative correlation between cell production rate and flowering time. Together, our data suggest a model in which variation in genome size is driven by natural selection on flowering time across altitudinal clines, connecting intraspecific variation in repetitive sequence to important differences in adaptive phenotypes.

Genome size in plants can vary by orders of magnitude, but this variation has long been considered to be of little functional consequence. Studying three independent adaptations to high altitude in Zea mays, we find that genome size experiences parallel pressures from natural selection, causing a reduction in genome size with increasing altitude. Though reductions in overall repetitive content are responsible for the genome size change, we find that only those individual loci contributing most to the variation in genome size are individually targeted by selection. To identify the phenotype influenced by genome size, we study how variation in genome size within a single wild population impacts leaf growth and cell division. We find that genome size variation correlates negatively with the rate of cell division, suggesting that individuals with larger genomes require longer to complete a mitotic cycle. Finally, we reanalyze data from maize inbreds to show that faster cell division is correlated with earlier flowering, connecting observed variation in genome size to an important adaptive phenotype.

Genomic abundance is not predictive of tandem repeat localization in grass genomes

Highly repetitive regions have historically posed a challenge when investigating sequence variation and content. High-throughput sequencing has enabled researchers to use whole-genome shotgun sequencing to estimate the abundance of repetitive sequence, and these methodologies have been recently applied to centromeres. Previous research has investigated variation in centromere repeats across eukaryotes, positing that the highest abundance tandem repeat in a genome is often the centromeric repeat. To test this assumption, we used shotgun sequencing and a bioinformatic pipeline to identify common tandem repeats across a number of grass species. We find that de novo assembly and subsequent abundance ranking of repeats can successfully identify tandem repeats with homology to known tandem repeats. Fluorescent in-situ hybridization shows that de novo assembly and ranking of repeats from non-model taxa identifies chromosome domains rich in tandem repeats both near pericentromeres and elsewhere in the genome.

Genetic Vulnerability and the Relationship of Commercial Germplasms of Maize in Brazil with the Nested Association Mapping Parents

A few breeding companies dominate the maize (Zea mays L.) hybrid market in Brazil: Monsanto^® (35%), DuPont Pioneer^® (30%), Dow Agrosciences^® (15%), Syngenta^® (10%) and Helix Sementes (4%). Therefore, it is important to monitor the genetic diversity in commercial germplasms as breeding practices, registration and marketing of new cultivars can lead to a significant reduction of the genetic diversity. Reduced genetic variation may lead to crop vulnerabilities, food insecurity and limited genetic gains following selection. The aim of this study was to evaluate the genetic vulnerability risk by examining the relationship between the commercial Brazilian maize germplasms and the Nested Association Mapping (NAM) Parents. For this purpose, we used the commercial hybrids with the largest market share in Brazil and the NAM parents. The hybrids were genotyped for 768 single nucleotide polymorphisms (SNPs), using the Illumina Goldengate^® platform. The NAM parent genomic data, comprising 1,536 SNPs for each line, were obtained from the Panzea data bank. The population structure, genetic diversity and the correlation between allele frequencies were analyzed. Based on the estimated effective population size and genetic variability, it was found that there is a low risk of genetic vulnerability in the commercial Brazilian maize germplasms. However, the genetic diversity is lower than those found in the NAM parents. Furthermore, the Brazilian germplasms presented no close relations with most NAM parents, except B73. This indicates that B73, or its heterotic group (Iowa Stiff Stalk Synthetic), contributed to the development of the commercial Brazilian germplasms.