Hiroshi Shinozuka

Fine-scale comparative genetic and physical mapping supports map-based cloning strategies for the self-incompatibility loci of perennial ryegrass (lolium perenne L.)

Perennial ryegrass is an obligate outbreeding pasture grass of the Poaceae family, with a two-locus (S and Z) gametophytic self-incompatibility (SI) mechanism. This system has provided a major obstacle to targeted varietal development, and enhanced knowledge is expected to support more efficient breeding strategies. Comparative genetics and physical mapping approaches have been developed to permit molecular cloning of the SI genes. SI gene-linked genetic markers based on heterologous cDNA restriction fragment length polymorphisms (RFLPs) and homologous genomic DNA-derived simple sequence repeats (SSRs) were converted to single nucleotide polymorphism (SNP) format for efficient genotyping. Genetic mapping identified the location of SI loci and demonstrated macrosynteny between related grass species. S- and Z-linked bacterial artificial chromosome (BAC) clones were sequenced using massively parallel pyrosequencing technology to provide the first physical mapping data for Poaceae SI loci. The sequence assembly process suggested a lower prevalence of middle repetitive sequences in the Z locus region and hence precedence for positional cloning strategy. In silico mapping using data from rice, Brachypodium distachyon and Sorghum revealed high sequence conservation in the vicinity of the Z locus region between SI and self-compatible (SC) grass species. Physical mapping identified a total of nine genes encoded in the Z locus region. Expression profiling and nucleotide diversity assessment identified two Z-linked genes, LpTC116908 and LpDUF247, as plausible candidates for the male and female determinants of the S-Z SI system.

Quantitative Trait Locus (QTL) meta-analysis and comparative genomics for candidate gene prediction in perennial ryegrass (Lolium perenne L.)

BackgroundIn crop species, QTL analysis is commonly used for identification of factors contributing to variation of agronomically important traits. As an important pasture species, a large number of QTLs have been reported for perennial ryegrass based on analysis of biparental mapping populations. Further characterisation of those QTLs is, however, essential for utilisation in varietal improvement programs.ResultsA bibliographic survey of perennial ryegrass trait-dissection studies identified a total of 560 QTLs from previously published papers, of which 189, 270 and 101 were classified as morphology-, physiology- and resistance/tolerance-related loci, respectively. The collected dataset permitted a subsequent meta-QTL study and implementation of a cross-species candidate gene identification approach. A meta-QTL analysis based on use of the BioMercator software was performed to identify two consensus regions for pathogen resistance traits. Genes that are candidates for causal polymorphism underpinning perennial ryegrass QTLs were identified through in silico comparative mapping using rice databases, and 7 genes were assigned to the p150/112 reference map. Markers linked to the Lp DGL1, Lp Ph1 and Lp PIPK1 genes were located close to plant size, leaf extension time and heading date-related QTLs, respectively, suggesting that these genes may be functionally associated with important agronomic traits in perennial ryegrass.ConclusionsFunctional markers are valuable for QTL meta-analysis and comparative genomics. Enrichment of such genetic markers may permit further detailed characterisation of QTLs. The outcomes of QTL meta-analysis and comparative genomics studies may be useful for accelerated development of novel perennial ryegrass cultivars with desirable traits.

Design of an F1 hybrid breeding strategy for ryegrasses based on selection of self-incompatibility locus-specific alleles

Relatively modest levels of genetic gain have been achieved in conventional ryegrass breeding when compared to cereal crops such as maize, current estimates indicating an annual improvement of 0.25–0.6% in dry matter production. This property is partially due to an inability to effectively exploit heterosis through the formation of F1 hybrids. Controlled crossing of ryegrass lines from geographically distant origins has demonstrated the occurrence of heterosis, which can result in increases of dry matter production in the order of 25%. Although capture of hybrid vigor offers obvious advantages for ryegrass cultivar production, to date there have been no effective and commercially suitable methods for obtaining high proportions of F1 hybrid seed. Continued advances in fine-scale genetic and physical mapping of the gametophytic self-incompatibility (SI) loci (S and Z) of ryegrasses are likely in the near future to permit the identification of closely linked genetic markers that define locus-specific haplotypes, allowing prediction of allelic variants and hence compatibility between different plant genotypes. Given the availability of such information, a strategy for efficient generation of ryegrass cultivars with a high proportion of F1 hybrid individuals has been simulated, which is suitable for commercial implementation. Through development of two parental pools with restricted diversity at the SI loci, relative crossing compatibility between pools is increased. Based on simulation of various levels of SI allele diversity restriction, the most effective scheme will generate 83.33% F1 hybrids. Results from the study, including the impact of varying flowering time, are discussed along with a proposed breeding design for commercial application.

A simple method for semi-random DNA amplicon fragmentation using the methylation-dependent restriction enzyme MspJI

BackgroundFragmentation at random nucleotide locations is an essential process for preparation of DNA libraries to be used on massively parallel short-read DNA sequencing platforms. Although instruments for physical shearing, such as the Covaris S2 focused-ultrasonicator system, and products for enzymatic shearing, such as the Nextera technology and NEBNext dsDNA Fragmentase kit, are commercially available, a simple and inexpensive method is desirable for high-throughput sequencing library preparation. MspJI is a recently characterised restriction enzyme which recognises the sequence motif CNNR (where R = G or A) when the first base is modified to 5-methylcytosine or 5-hydroxymethylcytosine.ResultsA semi-random enzymatic DNA amplicon fragmentation method was developed based on the unique cleavage properties of MspJI. In this method, random incorporation of 5-methyl-2’-deoxycytidine-5’-triphosphate is achieved through DNA amplification with DNA polymerase, followed by DNA digestion with MspJI. Due to the recognition sequence of the enzyme, DNA amplicons are fragmented in a relatively sequence-independent manner. The size range of the resulting fragments was capable of control through optimisation of 5-methyl-2’-deoxycytidine-5’-triphosphate concentration in the reaction mixture. A library suitable for sequencing using the Illumina MiSeq platform was prepared and processed using the proposed method. Alignment of generated short reads to a reference sequence demonstrated a relatively high level of random fragmentation.ConclusionsThe proposed method may be performed with standard laboratory equipment. Although the uniformity of coverage was slightly inferior to the Covaris physical shearing procedure, due to efficiencies of cost and labour, the method may be more suitable than existing approaches for implementation in large-scale sequencing activities, such as bacterial artificial chromosome (BAC)-based genome sequence assembly, pan-genomic studies and locus-targeted genotyping-by-sequencing.

Genotyping-by-sequencing through transcriptomics: implementation in a range of crop species with varying reproductive habits and ploidy levels

Summary The application of genomics in crops has the ability to significantly improve genetic gain for agriculture. Many marker‐dense tools have been developed, but few have seen broad adoption in plant genomics due to issues of significant variations of genome size, levels of ploidy, single nucleotide polymorphism (SNP) frequency and reproductive habit. When combined with limited breeding activities, small research communities and scant sequence resources, the suitability of popular systems is often suboptimal and routinely fails to effectively balance cost‐effectiveness and sample throughput. Genotyping‐by‐sequencing (GBS) encompasses a range of protocols including resequencing of the transcriptome. This study describes a skim GBS‐transcriptomics (GBS‐t) approach developed to be broadly applicable, cost‐effective and high‐throughput while still assaying a significant number of SNP loci. A range of crop species with differing levels of ploidy and degree of inbreeding/outbreeding were chosen, including perennial ryegrass, a diploid outbreeding forage grass; phalaris, a putative segmental allotetraploid outbreeding forage grass; lentil, a diploid inbreeding grain legume; and canola, an allotetraploid partially outbreeding oilseed. GBS‐t was validated as a simple and largely automated, cost‐effective method which generates sufficient SNPs (from 89 738 to 231 977) with acceptable levels of missing data and even genome coverage from c. 3 million sequence reads per sample. GBS‐t is therefore a broadly applicable system suitable for many crops, offering advantages over other systems. The correct choice of subsequent sequence analysis software is important, and the bioinformatics process should be iterative and tailored to the specific challenges posed by ploidy variation and extent of heterozygosity.

Comparative Genomics in Perennial Ryegrass (Lolium perenne L.): Identification and Characterisation of an Orthologue for the Rice Plant Architecture-Controlling Gene OsABCG5

Perennial ryegrass is an important pasture grass in temperate regions. As a forage biomass-generating species, plant architecture-related characters provide key objectives for breeding improvement. In silico comparative genomics analysis predicted colocation between a previously identified QTL for plant type (erect versus prostrate growth) and the ortholocus of the rice OsABCG5 gene (LpABCG5), as well as related QTLs in other Poaceae species. Sequencing of an LpABCG5-containing BAC clone identified presence of a paralogue (LpABCG6) in the vicinity of the LpABCG5 locus, in addition to three other gene-like sequences. Comparative genomics involving five other 5 grass species (rice, Brachypodium, sorghum, maize, and foxtail millet) revealed conserved microsynteny in the ABCG5 ortholocus-flanking region. Gene expression profiling and phylogenetic analysis suggested that the two paralogues are functionally distinct. Fourteen additional ABCG5 gene family members, which may interact with the LpABCG5 gene, were identified through sequencing of transcriptomes from perennial ryegrass leaf, anther, and pistils. A larger-scale phylogenetic analysis of the ABCG gene family suggested conservation between major branches of the Poaceae family. This study identified the LpABCG5 gene as a candidate for the plant type determinant, suggesting that manipulation of gene expression may provide valuable phenotypes for perennial ryegrass breeding.

Reference transcriptome assembly and annotation for perennial ryegrass

RNA-Seq methodology has been used to generate a comprehensive transcriptome sequence resource for perennial ryegrass, an important temperate pasture grass species. A total of 931 547 255 reads were obtained from libraries corresponding to 19 distinct tissue samples, including both vegetative and reproductive stages of development. Assembly of data generated a final filtered reference set of 48 713 contigs and scaffolds. The transcriptome resource will support whole genome sequence assembly, comparative genomics, implementation of genotyping-by-sequencing (GBS) methods based on transcript sampling, and identification of candidate genes for multiple biological functions.

Use of the melting curve assay as a means for high-throughput quantification of Illumina sequencing libraries

Background. Multiplexed sequencing is commonly performed on massively parallel short-read sequencing platforms such as Illumina, and the efficiency of library normalisation can affect the quality of the output dataset. Although several library normalisation approaches have been established, none are ideal for highly multiplexed sequencing due to issues of cost and/or processing time. Methods. An inexpensive and high-throughput library quantification method has been developed, based on an adaptation of the melting curve assay. Sequencing libraries were subjected to the assay using the Bio-Rad Laboratories CFX ConnectTM Real-Time PCR Detection System. The library quantity was calculated through summation of reduction of relative fluorescence units between 86 and 95 °C. Results.PCR-enriched sequencing libraries are suitable for this quantification without pre-purification of DNA. Short DNA molecules, which ideally should be eliminated from the library for subsequent processing, were differentiated from the target DNA in a mixture on the basis of differences in melting temperature. Quantification results for long sequences targeted using the melting curve assay were correlated with those from existing methods (R2 > 0.77), and that observed from MiSeq sequencing (R2 = 0.82). Discussion.The results of multiplexed sequencing suggested that the normalisation performance of the described method is equivalent to that of another recently reported high-throughput bead-based method, BeNUS. However, costs for the melting curve assay are considerably lower and processing times shorter than those of other existing methods, suggesting greater suitability for highly multiplexed sequencing applications.

Nucleotide diversity of vernalization and flowering-time-related genes in a germplasm collection of meadow fescue (Festuca pratensis Huds. syn. Lolium pratense (Huds.) Darbysh.).

In plant species, control of flowering time is an important factor for adaptation to local natural environments. The Vrn1, CO, FT1 and CK2α genes are key components in the flowering-specific signaling pathway of grass species. Meadow fescue is an agronomically important forage grass species, which is naturally distributed across Europe and Western Asia. In this study, meadow fescue flowering-time-related genes were resequenced to assess nucleotide diversity in European and Western Asian subpopulations. Identified sequence polymorphisms were then converted into PCR-based molecular genetic markers, and a meadow fescue germplasm collection was genotyped to investigate global allelic variation. Lower nucleotide diversities were observed for the Vrn1 and CO orthologs, while relatively higher values were observed for the FT1 and casein kinase II α-subunit (CK2α) orthologs. The nucleotide diversity for FT1 orthologs in the Western Asian subpopulation was significantly higher than those of the European subpopulation. Similarly, significant differences in nucleotide diversity for the remaining genes were observed between several combinations of subpopulation. The global allele distribution pattern was consistent with observed level of nucleotide diversity. These results suggested that the degree of purifying selection acting on the genes differs according to geographical location. As previously shown for model plant species, functional specificities of flowering-time-related genes may also vary according to environmental conditions.

Homology-based enzymatic DNA fragment assembly-based illumina sequencing library preparation

Abstract The current Illumina HiSeq and MiSeq platforms can generate paired-end reads of up to 2 x 250 bp and 2 x 300 bp in length, respectively. These read lengths may be substantially longer than genomic regions of interest when a DNA sequencing library is prepared through a target enrichment-based approach. A sequencing library preparation method has been developed based on the homology-based enzymatic DNA fragment assembly scheme to allow processing of multiple PCR products within a single read. Target sequences were amplified using locus-specific PCR primers with 8 bp tags, and using the tags, homology-based enzymatic DNA assembly was performed with DNA polymerase, T7 exonuclease and T4 DNA ligase. Short PCR amplicons can hence be assembled into a single molecule, along with sequencing adapters specific to the Illumina platforms. As a proof-of-concept experiment, short PCR amplicons (57–66 bp in length) derived from genomic DNA templates of field pea and containing variable nucleotide locations were assembled and sequenced on the MiSeq platform. The results were validated with other genotyping methods. When 5 PCR amplicons were assembled, 4.3 targeted sequences (single-nucleotide polymorphisms) on average were successfully identified within each read. The utility of this for sequencing of short fragments has consequently been demonstrated.