Rapid genome wide mapping of phosphine resistance loci by a simple regional averaging analysis in the red flour beetle, Tribolium castaneum

Background Next-generation sequencing technology is an important tool for the rapid, genome-wide identification of genetic variations. However, it is difficult to resolve the ‘signal’ of variations of interest and the ‘noise’ of stochastic sequencing and bioinformatic errors in the large datasets that are generated. We report a simple approach to identify regional linkage to a trait that requires only two pools of DNA to be sequenced from progeny of a defined genetic cross (i.e. bulk segregant analysis) at low coverage (<10×) and without parentage assignment of individual SNPs. The analysis relies on regional averaging of pooled SNP frequencies to rapidly scan polymorphisms across the genome for differential regional homozygosity, which is then displayed graphically. Results Progeny from defined genetic crosses of Tribolium castaneum (F4 and F19) segregating for the phosphine resistance trait were exposed to phosphine to select for the resistance trait while the remainders were left unexposed. Next generation sequencing was then carried out on the genomic DNA from each pool of selected and unselected insects from each generation. The reads were mapped against the annotated T. castaneum genome from NCBI (v3.0) and analysed for SNP variations. Since it is difficult to accurately call individual SNP frequencies when the depth of sequence coverage is low, variant frequencies were averaged across larger regions. Results from regional SNP frequency averaging identified two loci, tc_rph1 on chromosome 8 and tc_rph2 on chromosome 9, which together are responsible for high level resistance. Identification of the two loci was possible with only 5-7× average coverage of the genome per dataset. These loci were subsequently confirmed by direct SNP marker analysis and fine-scale mapping. Individually, homozygosity of tc_rph1 or tc_rph2 results in only weak resistance to phosphine (estimated at up to 1.5-2.5× and 3-5× respectively), whereas in combination they interact synergistically to provide a high-level resistance >200×. The tc_rph2 resistance allele resulted in a significant fitness cost relative to the wild type allele in unselected beetles over eighteen generations. Conclusion We have validated the technique of linkage mapping by low-coverage sequencing of progeny from a simple genetic cross. The approach relied on regional averaging of SNP frequencies and was used to successfully identify candidate gene loci for phosphine resistance in T. castaneum. This is a relatively simple and rapid approach to identifying genomic regions associated with traits in defined genetic crosses that does not require any specialised statistical analysis.

Genome-wide association analysis and differential expression analysis of resistance to Sclerotinia stem rot in Brassica napus

Brassica napus is one of the most important oil crops in the world, and stem rot caused by the fungus Sclerotinia sclerotiorum results in major losses in yield and quality. To elucidate resistance genes and pathogenesis-related genes, genome-wide association analysis of 347 accessions was performed using the Illumina 60K Brassica SNP (single nucleotide polymorphism) array. In addition, the detached stem inoculation assay was used to select five highly resistant (R) and susceptible (S) B. napus lines, 48 h postinoculation with S. sclerotiorum for transcriptome sequencing. We identified 17 significant associations for stem resistance on chromosomes A8 and C6, five of which were on A8 and 12 on C6. The SNPs identified on A8 were located in a 409-kb haplotype block, and those on C6 were consistent with previous QTL mapping efforts. Transcriptome analysis suggested that S. sclerotiorum infection activates the immune system, sulphur metabolism, especially glutathione (GSH) and glucosinolates in both R and S genotypes. Genes found to be specific to the R genotype related to the jasmonic acid pathway, lignin biosynthesis, defence response, signal transduction and encoding transcription factors. Twenty-four genes were identified in both the SNP-trait association and transcriptome sequencing analyses, including a tau class glutathione S-transferase (GSTU) gene cluster. This study provides useful insight into the molecular mechanisms underlying the plant's response to S. sclerotiorum.

Transcriptome analysis of Neotyphodium and Epichloë grass endophytes

Large-scale gene discovery has been performed for the grass fungal endophytes Neotyphodium coenophialum, Neotyphodium lolii, and Epichloë festucae. The resulting sequences have been annotated by comparison with public DNA and protein sequence databases and using intermediate gene ontology annotation tools. Endophyte sequences have also been analysed for the presence of simple sequence repeat and single nucleotide polymorphism molecular genetic markers. Sequences and annotation are maintained within a MySQL database that may be queried using a custom web interface. Two cDNA-based microarrays have been generated from this genome resource. They permit the interrogation of 3806 Neotyphodium genes (NchipTM microarray), and 4195 Neotyphodium and 920 Epichloë genes (EndoChipTM microarray), respectively. These microarrays provide tools for high-throughput transcriptome analysis, including genome-specific gene expression studies, profiling of novel endophyte genes, and investigation of the host grass–symbiont interaction. Comparative transcriptome analysis in Neotyphodium and Epichloë was performed

Applications of pedigree-based genome mapping in wheat and barley breeding programs.

The aim of the pedigree-based genome mapping project is to investigate and develop systems for implementing marker assisted selection to improve the efficiency of selection and increase the rate of genetic gain in breeding programs. Pedigree-based whole genome marker application provides a vehicle for incorporating marker technologies into applied breeding programs by bridging the gap between marker-trait association and marker implementation. We report on the development of protocols for implementation of pedigree-based whole genome marker analysis in breeding programs within the Australian northern winter cereals region. Examples of applications from the Queensland DPI&F wheat and barley breeding programs are provided, commenting on the use of microsatellites and other types of molecular markers for routine genomic analysis, the integration of genotypic, phenotypic and pedigree information for targeted wheat and barley lines, the genomic impacts of strong selection pressure in case study pedigrees, and directions for future pedigree-based marker development and analysis.

Completion of the genome sequence of Lettuce necrotic yellows virus, type species of the genus Cytorhabdovirus

We completed the genome sequence of Lettuce necrotic yellows virus (LNYV) by determining the nucleotide sequences of the 4a (putative phosphoprotein), 4b, M (matrix protein), G (glycoprotein) and L (polymerase) genes. The genome consists of 12,807 nucleotides and encodes six genes in the order 3′ leader-N-4a(P)-4b-M-G-L-5′ trailer. Sequences were derived from clones of a cDNA library from LNYV genomic RNA and from fragments amplified using reverse transcription-polymerase chain reaction. The 4a protein has a low isoelectric point characteristic for rhabdovirus phosphoproteins. The 4b protein has significant sequence similarities with the movement proteins of capillo- and trichoviruses and may be involved in cell-to-cell movement. The putative G protein sequence contains a predicted 25 amino acids signal peptide and endopeptidase cleavage site, three predicted glycosylation sites and a putative transmembrane domain. The deduced L protein sequence shows similarities with the L proteins of other plant rhabdoviruses and contains polymerase module motifs characteristic for RNA-dependent RNA polymerases of negative-strand RNA viruses. Phylogenetic analysis of this motif among rhabdoviruses placed LNYV in a group with other sequenced cytorhabdoviruses, most closely related to Strawberry crinkle virus.

Integrating sorghum whole genome sequence information with a compendium of sorghum QTL studies reveals uneven distribution of QTL and of gene-rich regions with significant implications for crop improvement.

A comprehensive analysis was conducted using 48 sorghum QTL studies published from 1995 to 2010 to make information from historical sorghum QTL experiments available in a form that could be more readily used by sorghum researchers and plant breeders. In total, 771 QTL relating to 161 unique traits from 44 studies were projected onto a sorghum consensus map. Confidence intervals (CI) of QTL were estimated so that valid comparisons could be made between studies. The method accounted for the number of lines used and the phenotypic variation explained by individual QTL from each study. In addition, estimated centimorgan (cM) locations were calculated for the predicted sorghum gene models identified in Phytozome (JGI GeneModels SBI v1.4) and compared with QTL distribution genome-wide, both on genetic linkage (cM) and physical (base-pair/bp) map scales. QTL and genes were distributed unevenly across the genome. Heterochromatic enrichment for QTL was observed, with approximately 22% of QTL either entirely or partially located in the heterochromatic regions. Heterochromatic gene enrichment was also observed based on their predicted cM locations on the sorghum consensus map, due to suppressed recombination in heterochromatic regions, in contrast to the euchromatic gene enrichment observed on the physical, sequence-based map. The finding of high gene density in recombination-poor regions, coupled with the association with increased QTL density, has implications for the development of more efficient breeding systems in sorghum to better exploit heterosis. The projected QTL information described, combined with the physical locations of sorghum sequence-based markers and predicted gene models, provides sorghum researchers with a useful resource for more detailed analysis of traits and development of efficient marker-assisted breeding strategies.

Recent emergence of the wheat Lr34 multi-pathogen resistance: insights from haplotype analysis in wheat, rice, sorghum and Aegilops tauschii

Spontaneous sequence changes and the selection of beneficial mutations are driving forces of gene diversification and key factors of evolution. In highly dynamic co-evolutionary processes such as plant-pathogen interactions, the plant's ability to rapidly adapt to newly emerging pathogens is paramount. The hexaploid wheat gene Lr34, which encodes an ATP-binding cassette (ABC) transporter, confers durable field resistance against four fungal diseases. Despite its extensive use in breeding and agriculture, no increase in virulence towards Lr34 has been described over the last century. The wheat genepool contains two predominant Lr34 alleles of which only one confers disease resistance. The two alleles, located on chromosome 7DS, differ by only two exon-polymorphisms. Putatively functional homoeologs and orthologs of Lr34 are found on the B-genome of wheat and in rice and sorghum, but not in maize, barley and Brachypodium. In this study we present a detailed haplotype analysis of homoeologous and orthologous Lr34 genes in genetically and geographically diverse selections of wheat, rice and sorghum accessions. We found that the resistant Lr34 haplotype is unique to the wheat D-genome and is not found in the B-genome of wheat or in rice and sorghum. Furthermore, we only found the susceptible Lr34 allele in a set of 252 Ae. tauschii genotypes, the progenitor of the wheat D-genome. These data provide compelling evidence that the Lr34 multi-pathogen resistance is the result of recent gene diversification occurring after the formation of hexaploid wheat about 8,000 years ago.

Complete genome sequence and integrated protein localization and interaction map for alfalfa dwarf virus, which combines properties of both cytoplasmic and nuclear plant rhabdoviruses

Summary We have determined the full-length 14,491-nucleotide genome sequence of a new plant rhabdovirus, alfalfa dwarf virus (ADV). Seven open reading frames (ORFs) were identified in the antigenomic orientation of the negative-sense, single-stranded viral RNA, in the order 3′-N-P-P3-M-G-P6-L-5′. The ORFs are separated by conserved intergenic regions and the genome coding region is flanked by complementary 3′ leader and 5′ trailer sequences. Phylogenetic analysis of the nucleoprotein amino acid sequence indicated that this alfalfa-infecting rhabdovirus is related to viruses in the genus Cytorhabdovirus. When transiently expressed as GFP fusions in Nicotiana benthamiana leaves, most ADV proteins accumulated in the cell periphery, but unexpectedly P protein was localized exclusively in the nucleus. ADV P protein was shown to have a homotypic, and heterotypic nuclear interactions with N, P3 and M proteins by bimolecular fluorescence complementation. ADV appears unique in that it combines properties of both cytoplasmic and nuclear plant rhabdoviruses.

Genetic differentiation analysis for the identification of complementary parental pools for sorghum hybrid breeding in Ethiopia

Key message The potential for exploiting heterosis for sorghum hybrid production in Ethiopia with improved local adaptation and farmers preferences has been investigated and populations suitable for initial hybrid development have been identified. Abstract Hybrids in sorghum have demonstrated increased productivity and stability of performance in the developed world. In Ethiopia, the uptake of hybrid sorghum has been limited to date, primarily due to poor adaptation and absence of farmer’s preferred traits in existing hybrids. This study aimed to identify complementary parental pools to develop locally adapted hybrids, through an analysis of whole genome variability of 184 locally adapted genotypes and introduced hybrid parents (R and B). Genetic variability was assessed using genetic distance, model-based STRUCTURE analysis and pair-wise comparison of groups. We observed a high degree of genetic similarity between the Ethiopian improved inbred genotypes and a subset of landraces adapted to lowland agro-ecology with the introduced R lines. This coupled with the genetic differentiation from existing B lines, indicated that these locally adapted genotype groups are expected to have similar patterns of heterotic expression as observed between introduced R and B line pools. Additionally, the hybrids derived from these locally adapted genotypes will have the benefit of containing farmers preferred traits. The groups most divergent from introduced B lines were the Ethiopian landraces adapted to highland and intermediate agro-ecologies and a subset of lowland-adapted genotypes, indicating the potential for increased heterotic response of their hybrids. However, these groups were also differentiated from the R lines, and hence are different from the existing complementary heterotic pools. This suggests that although these groups could provide highly divergent parental pools, further research is required to investigate the extent of heterosis and their hybrid performance.

Complete genome sequence and intracellular protein localization of Datura yellow vein nucleorhabdovirus

A limited number of plant rhabdovirus genomes have been fully sequenced, making taxonomic classification, evolutionary analysis and molecular characterization of this virus group difficult. We have for the first time determined the complete genome sequence of 13,188 nucleotides of Datura yellow vein nucleorhabdovirus (DYVV). DYVV genome organization resembles that of its closest relative, Sonchus yellow net virus (SYNV), with six ORFs in antigenomic orientation, separated by highly conserved intergenic regions and flanked by complementary 3′ leader and 5′ trailer sequences. As is typical for nucleorhabdoviruses, all viral proteins, except the glycoprotein, which is targeted to the endoplasmic reticulum, are localized to the nucleus. Nucleocapsid (N) protein, matrix (M) protein and polymerase, as components of nuclear viroplasms during replication, have predicted strong canonical nuclear localization signals, and N and M proteins exclusively localize to the nucleus when transiently expressed as GFP fusions. As in all nucleorhabdoviruses studied so far, N and phosphoprotein P interact when co-expressed, significantly increasing P nuclear localization in the presence of N protein. This research adds to the list of complete genomes of plant-infecting rhabdoviruses, provides molecular tools for further characterization and supports classification of DYVV as a nucleorhabdovirus closely related to but with some distinct differences from SYNV.