Validation of in silico-predicted genic SNPs in white clover (Trifolium repens L.), an outbreeding allopolyploid species

Emerging epigenetic mechanisms of long non-coding RNAs

Long non-coding RNAs (lncRNAs) have been increasingly appreciated as an integral component of gene regulatory networks. Genome-wide features of their origin and expression patterns ascribed a prominent role for lncRNAs to the regulation of protein-coding genes, and also suggest a potential link to many human diseases. Recent studies have begun to unravel the intricate regulatory mechanism of lncRNAs occurring at multiple levels. The brain is one of the richest sources of lncRNAs, many of which have already shown a close relationship with genes or genetic loci implicated in a wide range of neurological disorders. This review describes recently emerging mechanistic principles of lncRNA functions to provide neuroscientists with molecular insights that will help future research on lncRNAs in the brain.

An integrated genetic linkage map for white clover (Trifolium repens L.) with alignment to Medicago

BACKGROUND

White clover (Trifolium repens L.) is a temperate forage legume with an allotetraploid genome (2n=4×=32) estimated at 1093 Mb. Several linkage maps of various sizes, marker sources and completeness are available, however, no integrated map and marker set has explored consistency of linkage analysis among unrelated mapping populations. Such integrative analysis requires tools for homoeologue matching among populations. Development of these tools provides for a consistent framework map of the white clover genome, and facilitates in silico alignment with the model forage legume, Medicago truncatula.

RESULTS

This is the first report of integration of independent linkage maps in white clover, and adds to the literature on methyl filtered GeneThresher®-derived microsatellite (simple sequence repeat; SSR) markers for linkage mapping. Gene-targeted SSR markers were discovered in a GeneThresher® (TrGT) methyl-filtered database of 364,539 sequences, which yielded 15,647 SSR arrays. Primers were designed for 4,038 arrays and of these, 465 TrGT-SSR markers were used for parental consensus genetic linkage analysis in an F1 mapping population (MP2). This was merged with an EST-SSR consensus genetic map of an independent population (MP1), using markers to match homoeologues and develop a multi-population integrated map of the white clover genome. This integrated map (IM) includes 1109 loci based on 804 SSRs over 1274 cM, covering 97% of the genome at a moderate density of one locus per 1.2 cM. Eighteen candidate genes and one morphological marker were also placed on the IM. Despite being derived from disparate populations and marker sources, the component maps and the derived IM had consistent representations of the white clover genome for marker order and genetic length. In silico analysis at an E-value threshold of 1e-20 revealed substantial co-linearity with the Medicago truncatula genome, and indicates a translocation between T. repens groups 2 and 6 relative to M. truncatula.

CONCLUSIONS

This integrated genetic linkage analysis provides a consistent and comprehensive linkage analysis of the white clover genome, with alignment to a model forage legume. Associated marker locus information, particularly the homoeologue-specific markers, offers a new resource for forage legume research to enable genetic analysis and improvement of this forage and grassland species.

A hybrid next generation transcript sequencing-based approach to identify allelic and homeolog-specific single nucleotide polymorphisms in allotetraploid white clover

Background

White clover (Trifolium repens L.) is an allotetraploid species possessing two highly collinear ancestral sub-genomes. The apparent existence of highly similar homeolog copies for the majority of genes in white clover is problematic for the development of genome-based resources in the species. This is especially true for the development of genetic markers based on single nucleotide polymorphisms (SNPs), since it is difficult to distinguish between homeolog-specific and allelic variants. Robust methods for categorising single nucleotide variants as allelic or homeolog-specific in large transcript datasets are required. We illustrate one potential approach in this study.

Results

We used 454-pyrosequencing sequencing to generate ~760,000 transcript sequences from an 8th generation white clover inbred line. These were assembled and partially annotated to yield a reference transcript set comprising 71,545 sequences. We subsequently performed Illumina sequencing on three further white clover samples, generating 14 million transcript reads from a mixed sample comprising 24 divergent white clover genotypes, and 50 million reads on two further eighth generation white clover inbred lines. Mapping these reads to the reference transcript set allowed us to develop a significant SNP resource for white clover, and to partition the SNPs from the inbred lines into categories reflecting allelic or homeolog-specific variation. The potential for using haplotype reconstruction and progenitor genome comparison to assign haplotypes to specific ancestral sub-genomes of white clover is demonstrated for sequences corresponding to genes encoding dehydration responsive element binding protein and acyl-coA oxidase.

Conclusions

In total, 208,854 independent SNPs in 31,715 reference sequences were discovered, approximately three quarters of which were categorised as representing allelic or homeolog-specific variation using two inbred lines. This represents a significant resource for white clover genomics and genetics studies. We discuss the potential to extend the analysis to identify a “core set” of ancestrally derived homeolog specific variants in white clover.

Identification, characterization and interpretation of single-nucleotide sequence variation in allopolyploid crop species

An understanding of nature and extent of nucleotide sequence variation is required for programmes of discovery and characterization of single nucleotide polymorphisms (SNPs), which provide the most versatile class of molecular genetic marker. A majority of higher plant species are polyploids, and allopolyploidy, because of hybrid formation between closely related taxa, is very common. Mutational variation may arise both between allelic (homologous) sequences within individual subgenomes and between homoeologous sequences among subgenomes, in addition to paralogous variation between duplicated gene copies. Successful SNP validation in allopolyploids depends on differentiation of the sequence variation classes. A number of biological factors influence the feasibility of discrimination, including degree of gene family complexity, inbreeding or outbreeding reproductive habit, and the level of knowledge concerning progenitor diploid species. In addition, developments in high-throughput DNA sequencing and associated computational analysis provide general solutions for the genetic analysis of allopolyploids. These issues are explored in the context of experience from a range of allopolyploid species, representing grain (wheat and canola), forage (pasture legumes and grasses), and horticultural (strawberry) crop. Following SNP discovery, detection in routine genotyping applications also presents challenges for allopolyploids. Strategies based on either design of subgenome-specific SNP assays through homoeolocus-targeted polymerase chain reaction (PCR) amplification, or detection of incremental changes in nucleotide variant dosage, are described.

Identification of genetic factors influencing salt stress tolerance in white clover (Trifolium repens L.) by QTL analysis

Allotetraploid (2n = 4x = 32) white clover (Trifolium repens L.) is the most commonly cultivated legume component of temperate pastures, sown in swards with a companion grass species. Genetic control of growth performance of white clover on saline land is highly important for dairy industries, due to increasing soil salinity problems. The objective of this study was to identify quantitative trait loci (QTLs) for salinity tolerance in terms of vegetative growth under stress. Two parental genetic maps consisting of 213 and 159 marker loci and spanning 1,973.0 and 1,837.6 cM, respectively, were constructed using simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers from a two-way pseudo-test cross F(1) population derived from pair-crossing of the Haifa(2) and LCL(2) genotypes. A total of 8 unique genomic regions on 8 linkage groups (LGs) of the Haifa(2) parental map and 6 unique regions on 5 LGs in the LCL(2) parental map were associated with plant growth under salt stress and relative growth under stress, as compared to control conditions. The results of this study indicate that salt tolerance in white clover is controlled by multiple QTLs, some at common locations, but each of limited magnitude. Location of these QTLs provides the genetic basis and potential for pyramiding of salt tolerance genes in breeding improvement.

Sampling strategies for whole genome association studies in aquaculture and outcrossing plant species

A number of farmed species are characterized by breeding populations of large full-sib families, including aquaculture species and outcrossing plant species. Whole genome association studies in such species must account for stratification arising from the full-sib family structure to avoid high rates of false discovery. Here, we demonstrate the value of selective genotyping strategies which balance the contribution of families across high and low phenotypes to greatly reduce rates of false discovery with a minimal effect on power.

Genetic map construction and QTL mapping of resistance to blackleg (Leptosphaeria maculans) disease in Australian canola (Brassica napus L.) cultivars

Genetic map construction and identification of quantitative trait loci (QTLs) for blackleg resistance were performed for four mapping populations derived from five different canola source cultivars. Three of the populations were generated from crosses between single genotypes from the blackleg-resistant cultivars Caiman, Camberra and (AV)Sapphire and the blackleg-susceptible cultivar Westar(10). The fourth population was derived from a cross between genotypes from two blackleg resistant varieties (Rainbow and (AV)Sapphire). Different types of DNA-based markers were designed and characterised from a collection of 20,000 EST sequences generated from multiple Brassica species, including a new set of 445 EST-SSR markers of high value to the international community. Multiple molecular genetic marker systems were used to construct linkage maps with locus numbers varying between 219 and 468, and coverage ranging from 1173 to 1800 cM. The proportion of polymorphic markers assigned to map locations varied from 70 to 89% across the four populations. Publicly available simple sequence repeat markers were used to assign linkage groups to reference nomenclature, and a sub-set of mapped markers were also screened on the Tapidor x Ningyou (T x N) reference population to assist this process. QTL analysis was performed based on percentage survival at low and high disease pressure sites. Multiple QTLs were identified across the four mapping populations, accounting for 13-33% of phenotypic variance (V (p)). QTL-linked marker data are suitable for implementation in breeding for disease resistance in Australian canola cultivars. However, the likelihood of shifts in pathogen race structure across different geographical locations may have implications for the long-term durability of such associations.

EST, COSII, and arbitrary gene markers give similar estimates of nucleotide diversity in cultivated tomato (Solanum lycopersicum L.)

Because cultivated tomato (Solanum lycopersicum L.) is low in genetic diversity, public, verified single nucleotide polymorphism (SNP) markers within the species are in demand. To promote marker development we resequenced approximately 23 kb in a diverse set of 31 tomato lines including TA496. Three classes of markers were sampled: (1) 26 expressed-sequence tag (EST), all of which were predicted to be polymorphic based on TA496, (2) 14 conserved ortholog set II (COSII) or unigene, and (3) ten published sequences, composed of nine fruit quality genes and one anonymous RFLP marker. The latter two types contained mostly noncoding DNA. In total, 154 SNPs and 34 indels were observed. The distributions of nucleotide diversity estimates among marker types were not significantly different from each other. Ascertainment bias of SNPs was evaluated for the EST markers. Despite the fact that the EST markers were developed using SNP prediction within a sample consisting of only one TA496 allele and one additional allele, the majority of polymorphisms in the 26 EST markers were represented among the other 30 tomato lines. Fifteen EST markers with published SNPs were more closely examined for bias. Mean SNP diversity observations were not significantly different between the original discovery sample of two lines (53 SNPs) and the 31 line diversity panel (56 SNPs). Furthermore, TA496 shared its haplotype with at least one other line at 11 of the 15 markers. These data demonstrate that public EST databases and noncoding regions are a valuable source of unbiased SNP markers in tomato.

Comparison of genome structure between white clover and Medicago truncatula supports homoeologous group nomenclature based on conserved synteny

Computational analysis has been used to align the genetic map of white clover (Trifolium repens L.) with the draft genome sequence of the model legume species Medicago truncatula Gaertn. In silico comparison based on white clover expressed sequence tags that contain simple sequence repeat loci revealed substantial macrosynteny between the genomes of these two species, which are closely related within the Trifolieae tribe of the Fabaceae family. Six of the eight homoeologous chromosome groups (HGs) of allotetraploid white clover show predominant relationships with single M. truncatula (Mt) chromosomes, while the two remaining groups may have participated in an evolutionary reciprocal translocation event. On this basis, a new chromosome nomenclature system for allotetraploid white clover is proposed such that HG A = 3, HG B = 8, HG C = 7, HG D = 4, HG E = 1, HG F = 2, HG G = 5, and HG H = 6. A rationalized linkage map ordering system has also been demonstrated. Improved knowledge of the relationships between agricultural and model forage legume genomes will facilitate prediction of gene location for key agronomic traits for pasture production.

Identification of homologous, homoeologous and paralogous sequence variants in an outbreeding allopolyploid species based on comparison with progenitor taxa

The combination of homologous, homoeologous and paralogous classes of sequence variation presents major challenges for SNP discovery in outbreeding allopolyploid species. Previous in vitro gene-associated SNP discovery studies in the allotetraploid forage legume white clover (Trifolium repens L.) were vulnerable to such effects, leading to prohibitive levels of attrition during SNP validation. Identification of T. occidentale and T. pallescens as the putative diploid progenitors of white clover has permitted discrimination of the different sequence variant categories. Amplicons from selected abiotic stress tolerance-related genes were obtained using mapping family parents and individuals from each diploid species. Following cloning, progenitor comparison allowed tentative assignment of individual haplotypes to one or other sub-genome, as well as to gene copies within sub-genomes. A high degree of coincidence and identity between SNPs and HSVs was observed. Close similarity was observed between the genome of T. occidentale and one white clover sub-genome, but the affinity between T. pallescens and the other sub-genome was weaker, suggesting that a currently uncharacterised taxon may be the true second progenitor. Selected validated SNPs were attributed to individual sub-genomes by assignment to and naming of homoeologous linkage groups, providing the basis for improved genetic trait-dissection studies. The approach described in this study is broadly applicable to a range of allopolyploid taxa of equivocal ancestry.

Evidence on the molecular basis of the Ac/ac adaptive cyanogenesis polymorphism in white clover (Trifolium repens L)

White clover is polymorphic for cyanogenesis, with both cyanogenic and acyanogenic plants occurring in nature. This chemical defense polymorphism is one of the longest-studied and best-documented examples of an adaptive polymorphism in plants. It is controlled by two independently segregating genes: Ac/ac controls the presence/absence of cyanogenic glucosides; and Li/li controls the presence/absence of their hydrolyzing enzyme, linamarase. Whereas Li is well characterized at the molecular level, Ac has remained unidentified. Here we report evidence that Ac corresponds to a gene encoding a cytochrome P450 of the CYP79D protein subfamily (CYP79D15), and we describe the apparent molecular basis of the Ac/ac polymorphism. CYP79D orthologs catalyze the first step in cyanogenic glucoside biosynthesis in other cyanogenic plant species. In white clover, Southern hybridizations indicate that CYP79D15 occurs as a single-copy gene in cyanogenic plants but is absent from the genomes of ac plants. Gene-expression analyses by RT-PCR corroborate this finding. This apparent molecular basis of the Ac/ac polymorphism parallels our previous findings for the Li/li polymorphism, which also arises through the presence/absence of a single-copy gene. The nature of these polymorphisms may reflect white clover's evolutionary origin as an allotetraploid derived from cyanogenic and acyanogenic diploid progenitors.