Detection of alternative splicing in western corn rootworm (Diabrotica virgifera virgifera LeConte) in association with eCry3.1Ab resistance using RNA-seq and PacBio Iso-Seq.

Alternative splicing is a common feature in eukaryotes that not only increases the transcript diversity, but also has functional consequences. In insects, alternative splicing has been found associated with resistance to pesticides and Bt toxins. Up to date, the alternative splicing in western corn rootworm (Diabrotica virgifera virgifera LeConte) has not been studied. To investigate its alternative splicing pattern and relation to Bt resistance, we carried out single-molecule real-time (SMRT) transcript sequencing and Iso-seq analysis on resistant, eCry3.1Ab-selected and susceptible, unselected, western corn rootworm neonate midguts which fed on seedling maize with and without eCry3.1Ab for 12 and 24 h. We present transcriptome-wide alternative splicing patterns of western corn rootworm midgut in response to feeding on eCry3.1Ab-expressing corn using a comprehensive approach that combines both RNA-seq and SMRT transcript sequencing techniques. The results showed genes in western corn rootworm are highly alternatively spliced, which happens on 67.73% of multi-exon genes. One of the alternative splicing events we identified was a novel peritrophic matrix protein with two alternative splicing isoforms. Analysis of differential exon usage between resistant and susceptible colonies showed that in eCry3.1Ab-resistant western corn rootworm, expression of one isoform was significantly higher than in the susceptible colony, while no significant differences between colonies were observed with the other isoform. Our results provide the first survey of alternative splicing in western corn rootworm and suggest that the observed alternatively spliced isoforms of peritrophic matrix protein may be associated with eCry3.1Ab resistance in western corn rootworm.

Using online tools at the Bovine Genome Database to manually annotate genes in the new reference genome.

With the availability of a new highly contiguous Bos taurus reference genome assembly (ARS-UCD1.2), it is the opportune time to upgrade the bovine gene set by seeking input from researchers. Furthermore, advances in graphical genome annotation tools now make it possible for researchers to leverage sequence data generated with the latest technologies to collaboratively curate genes. For many years the Bovine Genome Database (BGD) has provided tools such as the Apollo genome annotation editor to support manual bovine gene curation. The goal of this paper is to explain the reasoning behind the decisions made in the manual gene curation process while providing examples using the existing BGD tools. We will describe the sources of gene annotation evidence provided at the BGD, including RNA-seq and Iso-Seq data. We will also explain how to interpret various data visualizations when curating gene models, and will demonstrate the value of manual gene annotation. The process described here can be applied to manual gene curation for other species with similar tools. With a better understanding of manual gene annotation, researchers will be encouraged to edit gene models and contribute to the enhancement of livestock gene sets.

Families of clustered microsatellites in a conifer genome.

Microsatellite clustering may account for genetic maps which do not coalesce into the expected number of linkage groups. Microsatellite organization within the large genome of Pinus taeda (1C = 20,000 Mb) was determined by (1) testing whether repeat motifs were sequestered within the low-copy DNA kinetic component and (2) testing for repeat motif clusters within DNA fragments regardless of copy number. Within the low-copy kinetic component, either (AC)n or (AG)n repeat units were present in 32% of sequences. No repeat motifs were found in the total genome control. Clustered repeat motifs were frequent; the (ATG)n triplet repeat motif was located upstream from a CG-rich trinucleotide microsatellite in 26 out of 44 microsatellite sequences. Fourteen of the clustered (ATG)n sequences could be assembled into four microsatellite sequence families based on similarities in the flanking regions. Consistent with the DNA turnover model, family members shared similar flanking regions but differed in repeat motif composition and length.

Retroelements contribute to the excess low-copy-number DNA in pine.

Excess DNA in the single-copy component is rarely recognized as a contributor to the C-value paradox yet the single-copy component of the pine genome is reported to comprise over 3000 Mb of DNA, in large excess over the estimated 100 Mb required for gene expression. Two hypotheses regarding the factors that might contribute to the excess low-copy-number DNA were tested. The first hypothesis proposes that the excess low-copy kinetic component is actually overestimated by reassociation data analysis. To test this, a previously published C0t curve for Pinus strobus was reanalyzed using a new estimate of genome size based on laser flow cytometry. Part of the excess low-copy-number DNA in the pine genome could be attributed to the choice of parameters used in the analysis of the reassociation data. The second hypothesis holds that diverged retrotransposons contribute to the excess low-copy DNA. Sequences randomly sampled from single-copy and low-repetitive kinetic components of the P. taeda genome were characterized. Twelve of 46 fragments cloned from these fractions were found to show sequence similarity to retroelements: hence diverged retroelements contribute to the excess low-repetitive kinetic component in the pine genome. Similarity search was shown to be a conservative method for identifying retroelements, and thus the number of retroelements in the low-copy component was actually underestimated. Most of the retroelements in this fraction were nonfunctional. divergent from known retroelement families and previously reported only for flowering plants. Divergent retrotransposons are thus a major factor contributing to the expansion of the low-repetitive DNA component in higher plants.

Low-copy microsatellite markers for Pinus taeda L.

Eighteen low-copy and genomic microsatellite markers were tested for Mendelian inheritance and then assayed in 41 Pinus taeda L. samples drawn from five regions in the southern United States. The PCR products had multiple alleles, high levels of polymorphism, and little non-specific priming. Fifteen of the 18 markers were informative for a P. taeda three-generation RFLP (restriction fragment length polymorphism) pedigree, and a P. taeda population survey revealed three to 28 alleles per locus. The highest allele numbers and polymorphic information content (PIC) values were associated with complex repeat sequences and (or) with sequences consisting of the longer strings of perfect repeats. The abundance of low- to rare-frequency alleles also accounted for high PIC values in both types of markers. Low-copy microsatellites are useful for the large, complex pine genome, especially in the absence of entire gene sequences in public databases and with the low levels of polymorphism in markers developed from expressed sequence tags (ESTs).

Microsatellite analysis of Pinus taeda L. in Zimbabwe.

Deducing the origin of early 20th century introductions of Pinus taeda into Zimbabwe is possible given microsatellite markers and clear population differentiation in ancestral U.S. populations. This study was designed to determine whether P. taeda introductions into Zimbabwe came from one U.S. region or whether the present-day population is an admixture of introductions from east and west of the Mississippi River Valley. Principal components analysis, Cavalli-Sforza and Edwards' chord distances and presence of diagnostic alleles each indicate that the Zimbabwe population is an admixture. There were five novel alleles in the Zimbabwe population not represented in the indigenous U.S. populations, possibly because of de novo mutation, introgression with other introduced North American pines or sampling error.

Leaf-level and whole-plant gas exchange characteristics of shortleaf pine exposed to ozone and simulated acid rain.

Field-grown shortleaf pine (Pinus echinata Mill.) seedlings were exposed to ozone (O(3)) and simulated acid rain (SAR) in open-top chambers over three growing seasons. Ranges of O(3) and SAR spanned ambient levels found in the southern USA. Effects of O(3) on leaf-level and whole-plant gas exchange were characterized for a single measurement period immediately before the third summer of exposure. Decreased photosynthesis rates were attributed to O(3), but not SAR. Stomatal conductance decreased in response to O(3) exposure, and either increased or was unaffected by SAR. Increased internal CO(2) concentration (c(i)) in response to O(3) treatment indicated a greater effect of O(3) on photosynthetic capacity than stomatal conductance. Whole-seedling gas exchange characteristics indicated that whole-seedling carbon assimilation was more severely affected by O(3) than was evident from leaf-level gas exchange characteristics. Seedlings exposed to O(3) retained fewer flushes than seedlings grown in charcoal-filtered air.