The utility of PacBio circular consensus sequencing for characterizing complex gene families in non-model organisms

Background

Molecular characterization of highly diverse gene families can be time consuming, expensive, and difficult, especially when considering the potential for relatively large numbers of paralogs and/or pseudogenes. Here we investigate the utility of Pacific Biosciences single molecule real-time (SMRT) circular consensus sequencing (CCS) as an alternative to traditional cloning and Sanger sequencing PCR amplicons for gene family characterization. We target vomeronasal gene receptors, one of the most diverse gene families in mammals, with the goal of better understanding intra-specific V1R diversity of the gray mouse lemur (Microcebus murinus). Our study compares intragenomic variation for two V1R subfamilies found in the mouse lemur. Specifically, we compare gene copy variation within and between two individuals of M. murinus as characterized by different methods for nucleotide sequencing. By including the same individual animal from which the M. murinus draft genome was derived, we are able to cross-validate gene copy estimates from Sanger sequencing versus CCS methods.

Results

We generated 34,088 high quality circular consensus sequences of two diverse V1R subfamilies (here referred to as V1RI and V1RIX) from two individuals of Microcebus murinus. Using a minimum threshold of 7× coverage, we recovered approximately 90% of V1RI sequences previously identified in the draft M. murinus genome (59% being identical at all nucleotide positions). When low coverage sequences were considered (i.e. < 7× coverage) 100% of V1RI sequences identified in the draft genome were recovered. At least 13 putatively novel V1R loci were also identified using CCS technology.

Conclusions

Recent upgrades to the Pacific Biosciences RS instrument have improved the CCS technology and offer an alternative to traditional sequencing approaches. Our results suggest that the Microcebus murinus V1R repertoire has been underestimated in the draft genome. In addition to providing an improved understanding of V1R diversity in the mouse lemur, this study demonstrates the utility of CCS technology for characterizing complex regions of the genome. We anticipate that long-read sequencing technologies such as PacBio SMRT will allow for the assembly of multigene family clusters and serve to more accurately characterize patterns of gene copy variation in large gene families, thus revealing novel micro-evolutionary patterns within non-model organisms.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-720) contains supplementary material, which is available to authorized users.     

Construction of Pseudomolecule Sequences of the aus Rice Cultivar Kasalath for Comparative Genomics of Asian Cultivated Rice

Having a deep genetic structure evolved during its domestication and adaptation, the Asian cultivated rice (Oryza sativa) displays considerable physiological and morphological variations. Here, we describe deep whole-genome sequencing of the aus rice cultivar Kasalath by using the advanced next-generation sequencing (NGS) technologies to gain a better understanding of the sequence and structural changes among highly differentiated cultivars. The de novo assembled Kasalath sequences represented 91.1% (330.55 Mb) of the genome and contained 35 139 expressed loci annotated by RNA-Seq analysis. We detected 2 787 250 single-nucleotide polymorphisms (SNPs) and 7393 large insertion/deletion (indel) sites (>100 bp) between Kasalath and Nipponbare, and 2 216 251 SNPs and 3780 large indels between Kasalath and 93-11. Extensive comparison of the gene contents among these cultivars revealed similar rates of gene gain and loss. We detected at least 7.39 Mb of inserted sequences and 40.75 Mb of unmapped sequences in the Kasalath genome in comparison with the Nipponbare reference genome. Mapping of the publicly available NGS short reads from 50 rice accessions proved the necessity and the value of using the Kasalath whole-genome sequence as an additional reference to capture the sequence polymorphisms that cannot be discovered by using the Nipponbare sequence alone.     

Effect of cleavage enzyme, search algorithm and decoy database on mass spectrometric identification of wheat gluten proteins

While tandem mass spectrometry (MS/MS) is routinely used to identify proteins from complex mixtures, certain types of proteins present unique challenges for MS/MS analyses. The major wheat gluten proteins, gliadins and glutenins, are particularly difficult to distinguish by MS/MS. Each of these groups contains many individual proteins with similar sequences that include repetitive motifs rich in proline and glutamine. These proteins have few cleavable tryptic sites, often resulting in only one or two tryptic peptides that may not provide sufficient information for identification. Additionally, there are less than 14,000 complete protein sequences from wheat in the current NCBInr release. In this paper, MS/MS methods were optimized for the identification of the wheat gluten proteins. Chymotrypsin and thermolysin as well as trypsin were used to digest the proteins and the collision energy was adjusted to improve fragmentation of chymotryptic and thermolytic peptides. Specialized databases were constructed that included protein sequences derived from contigs from several assemblies of wheat expressed sequence tags (ESTs), including contigs assembled from ESTs of the cultivar under study. Two different search algorithms were used to interrogate the database and the results were analyzed and displayed using a commercially available software package (Scaffold). We examined the effect of protein database content and size on the false discovery rate. We found that as database size increased above 30,000 sequences there was a decrease in the number of proteins identified. Also, the type of decoy database influenced the number of proteins identified. Using three enzymes, two search algorithms and a specialized database allowed us to greatly increase the number of detected peptides and distinguish proteins within each gluten protein group.     

5S ribosomal gene clusters in wheat: pulsed field gel electrophoresis reveals a high degree of polymorphism

Summary The long-range structure of 5S rRNA gene clusters has been investigated in wheat (Triticum aestivum L.) by means of pulsed field gel electrophoresis. Using aneuploid stocks, 5S rRNA gene clusters were assigned to sites on chromosomes 1B, 1D, 513 and 5D. Cluster sizes were evaluated and the copy number of 5S DNA repeats was estimated at 4700-5200 copies for the short repeating unit (410 bp) and about 3100 copies for the long repeat (500 bp) per haploid genome. A comparison of wheat cultivars revealed extremely high levels of polymorphism in the 5S rRNA gene clusters. With one restriction enzyme digest all varieties tested gave unique banding patterns and, on a per fragment basis, 21-fold more polymorphism was detected among cultivars for 5S DNA compared to standard restriction fragment length polymorphisms (RFLPs) detected with single copy clones. Experiments with aneuploid stocks suggest that the 5S rRNA gene clusters at several chromosomal sites contribute to this polymorphism. A number of previous reports have shown that wheat cultivars are not easily distinguished by isozymes or RFLPs. The high level of variation detected in 5S rRNA gene clusters therefore offers the possibility of a sensitive fingerprinting method for wheat. 5S DNA and other macro-satellite sequences may also serve as hypervariable Mendelian markers for genetic and breeding experiments in wheat.     

小麦高分子量谷蛋白亚基及其基因的研究进展

主要介绍了小麦高分子量谷蛋白亚基（HMW－GS）及其基因的研究进展情况,目前,转基因小麦的技术已经逐渐成熟,由于分子生物学领域分子标记技术的迅速发展,尤其是PCR技术的广泛应用,为实现外源优良储藏蛋白基因导入改良品种提供了可能,利用已知小麦品种的基因序列设计引物,从众多的未知小麦品种中扩增出新基因加以研究并做外源优质HMW－GS基因的转入已成为一种趋势。     

Cell organisation,sulphur metabolism and ion transport-related genes are differentially expressed in Paracoccidioides brasiliensis mycelium and yeast cells

Background

Bread wheat (Triticum aestivum) is an important staple food. However, wheat gluten proteins cause celiac disease (CD) in 0.5 to 1% of the general population. Among these proteins, the α-gliadins contain several peptides that are associated to the disease.

Results

We obtained 230 distinct α-gliadin gene sequences from severaldiploid wheat species representing the ancestral A, B, and D genomes of the hexaploid bread wheat. The large majority of these sequences (87%) contained an internal stop codon. All α-gliadin sequences could be distinguished according to the genome of origin on the basis of sequence similarity, of the average length of the polyglutamine repeats, and of the differences in the presence of four peptides that have been identified as T cell stimulatory epitopes in CD patients through binding to HLA-DQ2/8. By sequence similarity, α-gliadins from the public database of hexaploid T. aestivum could be assigned directly to chromosome 6A, 6B, or 6D. T. monococcum (A genome) sequences, as well as those from chromosome 6A of bread wheat, almost invariably contained epitope glia-α9 and glia-α20, but never the intact epitopes glia-α and glia-α2. A number of sequences from T. speltoides, as well as a number of sequences fromchromosome 6B of bread wheat, did not contain any of the four T cell epitopes screened for. The sequences from T. tauschii (D genome), as well as those from chromosome 6D of bread wheat, were found to contain all of these T cell epitopes in variable combinations per gene. The differences in epitope composition resulted mainly from point mutations. These substitutions appeared to be genome specific.

Conclusion

Our analysis shows that α-gliadin sequences from the three genomes of bread wheat form distinct groups. The four known T cell stimulatory epitopes are distributed non-randomly across the sequences, indicating that the three genomes contribute differently to epitope content. A systematic analysis of all known epitopes in gliadins and glutenins will lead to better understanding of the differences in toxiCity among wheat varieties. On the basis of such insight, breeding strategies can be designed to generate less toxic varieties of wheat which may be tolerated by at least part of the CD patient population.     

Distribution,functional impact,and origin mechanisms of copy number variation in the barley genome

Background

There is growing evidence for the prevalence of copy number variation (CNV) and its role in phenotypic variation in many eukaryotic species. Here we use array comparative genomic hybridization to explore the extent of this type of structural variation in domesticated barley cultivars and wild barleys.

Results

A collection of 14 barley genotypes including eight cultivars and six wild barleys were used for comparative genomic hybridization. CNV affects 14.9% of all the sequences that were assessed. Higher levels of CNV diversity are present in the wild accessions relative to cultivated barley. CNVs are enriched near the ends of all chromosomes except 4H, which exhibits the lowest frequency of CNVs. CNV affects 9.5% of the coding sequences represented on the array and the genes affected by CNV are enriched for sequences annotated as disease-resistance proteins and protein kinases. Sequence-based comparisons of CNV between cultivars Barke and Morex provided evidence that DNA repair mechanisms of double-strand breaks via single-stranded annealing and synthesis-dependent strand annealing play an important role in the origin of CNV in barley.

Conclusions

We present the first catalog of CNVs in a diploid Triticeae species, which opens the door for future genome diversity research in a tribe that comprises the economically important cereal species wheat, barley, and rye. Our findings constitute a valuable resource for the identification of CNV affecting genes of agronomic importance. We also identify potential mechanisms that can generate variation in copy number in plant genomes.     

Genomic regions influencing gene expression of the HMW glutenins in wheat

Bread wheat (Triticum aestivum L.) produces glutenin storage proteins in the endosperm. The HMW glutenins confer distinct viscoelastic properties to bread dough. The genetics of HMW glutenin proteins have been extensively studied, and information has accumulated about individual subunits, chromosomal locations and DNA sequences, but little is known about the regulators of the HMW glutenins. This investigation addressed the question of glutenin regulators. Expression of the glutenins was analyzed using QRT-PCR in ditelosomic (dt) Chinese Spring (CS) lines. Primers were designed for each of 4 CS glutenin genes and a control, non-storage protein endosperm-specific gene Agp-L (ADP-glucose pyrophosphorylase). Each line represents CS wheat, lacking one chromosome arm. The effect of a missing arm could feasibly cause an increase, decrease or no change in expression. For each HMW glutenin, results indicated there were, on average, 8 chromosome arms with an up-regulatory effect and only one instance of a down-regulatory effect. There were significant correlations between orthologous and paralogous HMW glutenins for effects of chromosome groups B and D. Some or all the glutenin alleles shared regulatory loci on chromosome arms 2BS, 7BS, 4DS, 5DS and 6DS, and Agp-L shared regulatory loci with glutenins on arms 7AS, 7BS, 2DS, 3DS, 4DS and 5DS. These results suggest a few chromosome arms contain putative regulatory genes affecting the expression of conserved cis elements of 4 HMW glutenin and Agp-L genes in CS. Regulation by common genes implies the regulators have diverged little from the common wheat ancestor, and furthermore, some regulation may be shared by endosperm-specific-genes. Significant common regulators have practical implications.     

Conserved globulin gene across eight grass genomes identify fundamental units of the loci encoding seed storage proteins

The wheat high molecular weight (HMW) glutenins are important seed storage proteins that determine bread-making quality in hexaploid wheat (Triticum aestivum). In this study, detailed comparative sequence analyses of large orthologous HMW glutenin genomic regions from eight grass species, representing a wide evolutionary history of grass genomes, reveal a number of lineage-specific sequence changes. These lineage-specific changes, which resulted in duplications, insertions, and deletions of genes, are the major forces disrupting gene colinearity among grass genomes. Our results indicate that the presence of the HMW glutenin gene in Triticeae genomes was caused by lineage-specific duplication of a globulin gene. This tandem duplication event is shared by Brachypodium and Triticeae genomes, but is absent in rice, maize, and sorghum, suggesting the duplication occurred after Brachypodium and Triticeae genomes diverged from the other grasses ~35 Ma ago. Aside from their physical location in tandem, the sequence similarity, expression pattern, and conserved cis-acting elements responsible for endosperm-specific expression further support the paralogous relationship between the HMW glutenin and globulin genes. While the duplicated copy in Brachypodium has apparently become nonfunctional, the duplicated copy in wheat has evolved to become the HMW glutenin gene by gaining a central prolamin repetitive domain.     

Characterisation of two gene subunits on the 1R chromosome of rye as orthologs of each of the Glu-1 genes of hexaploid wheat

The visco-elastic properties of bread flour are firmly associated with the presence or absence of certain HMW subunits coded by the Glu-1 genes. Identifying allelic specific molecular markers (AS-PCR) associated with the presence of Glu-1 genes can serve as a valuable tool for the selection of useful genotypes. This paper reports the use of primers designed from nucleotide sequences of the Glu-D1 gene of wheat (AS-PCR for Glu-D1y10) that recognise and amplify homologous sequences of the Glu-R1 gene subunits of rye. The primers amplify the complete coding regions and provided two products of different size in rye, in wheats carrying the substitution 1R(1D) and in rye-wheat aneuploid lines carrying the long arm of chromosome 1R. The location, the molecular characterisation of these sequences and their expression during grain ripening seem to demonstrate that the amplification products correspond to structural genes encoding the high-molecular-weight (HMW) glutenins of rye. The homology of the rye gene to subunits encoding HMW glutenins in wheat was confirmed by Southern blots and sequencing. The amplification-products were cloned, sequenced and characterised, and the sequences compared with the main glutenin subunits of wheat and related species. Further, an RT-PCR experiment was performed using primers designed from the sequence of both amplified products. This assay demonstrated that both sequences are expressed in endosperm during grain ripening. The results of these analyses suggest that both gene subunits correspond to x- and y-type genes of the Glu-R1 locus of rye. Received: 11 December 2000 / Accepted: 17 April 2001     

Towards coeliac‐safe bread

Gluten‐free foods cannot substitute for products made from wheat flour. When wheat products are digested, the remaining peptides can trigger an autoimmune disease in 1% of the North American and European population, called coeliac disease. Because wheat proteins are encoded by a large gene family, it has been impossible to use conventional breeding to select wheat varieties that are coeliac‐safe. However, one can test the properties of protein variants by expressing single genes in coeliac‐safe cereals like maize. One source of protein that can be considered as coeliac‐safe and has bread‐making properties is teff (Eragrostis tef), a grain consumed in Ethiopia. Here, we show that teff α‐globulin3 (Etglo3) forms storage vacuoles in maize that are morphologically similar to those of wheat. Using transmission electron microscopy, immunogold labelling shows that Etglo3 is almost exclusively deposited in the storage vacuole as electron‐dense aggregates. Of maize seed storage proteins, 27‐kDa γ‐zein is co‐deposited with Etglo3. Etglo3 polymerizes via intermolecular disulphide bonds in maize, similar to wheat HMW glutenins under non‐reducing conditions. Crossing maize Etglo3 transgenic lines with α‐, β‐ and γ‐zein RNA interference (RNAi) lines reveals that Etglo3 accumulation is only dramatically reduced in γ‐zein RNAi background. This suggests that Etglo3 and 27‐kDa γ‐zein together cause storage vacuole formation and behave similar to the interactions of glutenins and gliadins in wheat. Therefore, expression of teff α‐globulins in maize presents a major step in the development of a coeliac‐safe grain with bread‐making properties.     

Evidence for the presence of two different types of protein bodies in wheat endosperm

Storage proteins of wheat grains (Triticum L. em Thell) are deposited in protein bodies inside vacuoles. However, the subcellular sites and mechanisms of their aggregation into protein bodies are not clear. In the present report, we provide evidence for two different types of protein bodies, low- and high-density types that accumulate concurrently and independently in developing wheat endosperm cells. Gliadins were present in both types of protein bodies, whereas the high molecular weight glutenins were localized mainly in the dense ones. Pulse-chase experiments verified that the dense protein bodies were not formed by a gradual increase in density but, presumably, by a distinct, quick process of storage protein aggregation. Subcellular fractionation and electron microscopy studies revealed that the wheat homolog of immunoglobulin heavy-chain-binding protein, an endoplasmic reticulum-resident protein, was present within the dense protein bodies, implying that these were formed by aggregation of storage proteins within the endoplasmic reticulum. The present results suggest that a large part of wheat storage proteins aggregate into protein bodies within the rough endoplasmic reticulum. Because these protein bodies are too large to enter the Golgi, they are likely to be transported directly to vacuoles. This route may operate in concert with the known Golgi-mediated transport to vacuoles in which the storage proteins apparently condense into protein bodies at a postendoplasmic reticulum location. Our results further suggest that although gliadins are transported by either one of these routes, the high molecular weight glutenins use only the Golgi bypass route.     

Characterization of high molecular weight glutenin subunit genes from the Ns genome of Psathyrostachys juncea

The Ns genome of the genus Psathyrostachys possesses superior traits useful for wheat improvement. However, very little is known about the high molecular weight (HMW) subunits of glutenin encoded by the Ns genome. In this paper, we report the isolation of four alleles of HMW glutenin subunit gene from Psathyrostachys juncea. Sequence alignment data shows the four alleles have similar primary structure with those in wheat and other wheat-related grasses, with some unique modifications. All four sequences more closely resemble y-type, rather than x-type, glutenins. However, our results show three of the subunits (1Ns2-4) contain an extra glutamine residue in the N-terminal region not found on typical y-type subunits, as well as the x-type subunit specific sequence LAAQLPAMCRL. These three subunits likely represent an intermediate state in the divergence between x- and y-type subunits. Results also indicate that the Ns genome is more closely related to the St genome of Pseudoroegneria than any other Triticeae genomes.     

Copy number variation of transposable elements in Triticum–Aegilops genus suggests evolutionary and revolutionary dynamics following allopolyploidization

Key message

Here, we report on copy number variation of transposable elements and on the genome-specific proliferation in wheat. In addition, we report on revolutionary and evolutionary dynamics of transposons.

Abstract

Wheat is a valuable model for understanding the involvement of transposable elements (TEs) in speciation as wheat species (Triticum–Aegilops group) have diverged from a common ancestor, have undergone two events of speciation through allopolyploidy, and contain a very high fraction of TEs. However, an unbiased genome-wide examination of TE variation among these species has not been conducted. Our research utilized quantitative real time PCR to assess the relative copy numbers of 16 TE families in various Triticum and Aegilops species. We found (1) high variation and genome-specificity of TEs in wheat species, suggesting they were active throughout the evolution of wheat, (2) neither Ae. searsii nor Ae. speltoides by themselves can be the only contributors of the B genome to wheat, and (3) nonadditive changes in TE quantities in polyploid wheat. This study indicates the apparent involvement of large TEs in creating genetic variation in revolutionary and evolutionary scales following allopolyploidization events, presumably assisting in the diploidization of homeologous chromosomes.     

Nucleic acid sequence and chromosome assignment of a wheat storage protein gene.

A cloned gliadin gene was isolated from a wheat genomic library, and 2.4 kb of its primary sequence determined. The gene, alpha-1Y, was found by Southern analysis to be located on chromosome 6A, and its derived amino acid sequence identifies it as a member of the A-gliadin subgroup of alpha-gliadins located on the short arm of that chromosome. alpha-1Y is apparently functional, and contains consensus TATA and CAAT boxes, and polyadenylation signals. This gliadin gene has no introns, and its noncoding flanking regions contain several short repeats and inverted sequences. The gene is contained in a 6.2 kb EcoRI genomic fragment whose apparent copy number varies in different wheat cultivars.     

Wheat cultivar-specific proteins in grain revealed by 2-DE and their application to cultivar identification of flour

Wheat flour proteins were studied to identify the cultivar-specific proteins and use them to identify cultivars in flours. Proteins extracted from flours of Japanese wheat (cultivars Hokushin, Horoshirikomugi, Kitanokaori and Kachikei 33) and Canadian wheat (Canada Western Red Spring Wheat No. 1; 1CW) were analyzed by 2-DE with IEF gels over three pH ranges: pH 4-7, pH 5-8, and pH 6-11. This system enabled detection of more than 1600 protein spots. We recognized that among 50 protein spots showing cultivar-dependent qualitative changes, 25 proteins were wheat cultivar specific. These 50 protein spots were analyzed by N-terminal Edman degradation microsequencing and MALDI-TOF-MS; 21 protein spots were storage proteins, such as gliadin and low-molecular mass glutenin subunit. Five protein spots were identified as dehydroascorbate reductase (Triticum aestivum), triticin precursor (T. aestivum), alpha-amylase inhibitor (Oryza sativa), DNA-binding with one finger (Dof) zinc family protein (O. sativa), and nonphototropic hypocotyl 1 (NPH1) protein (Avena sativa). The other protein spots appeared to be hypothetical proteins (O. sativa or Arabidopsis thaliana) or functional unknown proteins. These specific proteins can be used as markers to identify wheat cultivars in blended flour composed of two or three flours.