Whole-genome sequencing reveals untapped genetic potential in Africa’s indigenous cereal crop sorghum

Sorghum is a food and feed cereal crop adapted to heat and drought and a staple for 500 million of the world’s poorest people. Its small diploid genome and phenotypic diversity make it an ideal C4 grass model as a complement to C3 rice. Here we present high coverage (16-45 × ) resequenced genomes of 44 sorghum lines representing the primary gene pool and spanning dimensions of geographic origin, end-use and taxonomic group. We also report the first resequenced genome of S. propinquum, identifying 8 M high-quality SNPs, 1.9 M indels and specific gene loss and gain events in S. bicolor. We observe strong racial structure and a complex domestication history involving at least two distinct domestication events. These assembled genomes enable the leveraging of existing cereal functional genomics data against the novel diversity available in sorghum, providing an unmatched resource for the genetic improvement of sorghum and other grass species.

Domestication and the storage starch biosynthesis pathway: Signatures of selection from a whole sorghum genome sequencing strategy

Next-generation sequencing of complete genomes has given researchers unprecedented levels of information to study the multifaceted evolutionary changes that have shaped elite plant germplasm. In conjunction with population genetic analytical techniques and detailed online databases, we can more accurately capture the effects of domestication on entire biological pathways of agronomic importance. In this study, we explore the genetic diversity and signatures of selection in all predicted gene models of the storage starch synthesis pathway of Sorghum bicolor, utilizing a diversity panel containing lines categorized as either ‘Landraces’ or ‘Wild and Weedy’ genotypes. Amongst a total of 114 genes involved in starch synthesis, 71 had at least a single signal of purifying selection and 62 a signal of balancing selection and others a mix of both. This included key genes such as STARCH PHOSPHORYLASE 2 (SbPHO2, under balancing selection), PULLULANASE (SbPUL, under balancing selection) and ADP-glucose pyrophosphorylases (SHRUNKEN2, SbSH2 under purifying selection). Effectively, many genes within the primary starch synthesis pathway had a clear reduction in nucleotide diversity between the Landraces and wild and weedy lines indicating that the ancestral effects of domestication are still clearly identifiable. There was evidence of the positional rate variation within the well-characterized primary starch synthesis pathway of sorghum, particularly in the Landraces, whereby low evolutionary rates upstream and high rates downstream in the metabolic pathway were expected. This observation did not extend to the wild and weedy lines or the minor starch synthesis pathways.

Domestication and the storage starch biosynthesis pathway: signatures of selection from a whole sorghum genome sequencing strategy

Next-generation sequencing of complete genomes has given researchers unprecedented levels of information to study the multifaceted evolutionary changes that have shaped elite plant germplasm. In conjunction with population genetic analytical techniques and detailed online databases, we can more accurately capture the effects of domestication on entire biological pathways of agronomic importance. In this study, we explore the genetic diversity and signatures of selection in all predicted gene models of the storage starch synthesis pathway of Sorghum bicolor, utilizing a diversity panel containing lines categorized as either ‘Landraces’ or ‘Wild and Weedy’ genotypes. Amongst a total of 114 genes involved in starch synthesis, 71 had at least a single signal of purifying selection and 62 a signal of balancing selection and others a mix of both. This included key genes such as STARCH PHOSPHORYLASE 2 (SbPHO2, under balancing selection), PULLULANASE (SbPUL, under balancing selection) and ADP-glucose pyrophosphorylases (SHRUNKEN2, SbSH2 under purifying selection). Effectively, many genes within the primary starch synthesis pathway had a clear reduction in nucleotide diversity between the Landraces and wild and weedy lines indicating that the ancestral effects of domestication are still clearly identifiable. There was evidence of the positional rate variation within the well-characterized primary starch synthesis pathway of sorghum, particularly in the Landraces, whereby low evolutionary rates upstream and high rates downstream in the metabolic pathway were expected. This observation did not extend to the wild and weedy lines or the minor starch synthesis pathways.

Optimizing Hill Seeding Density for High-Yielding Hybrid Rice in a Single Rice Cropping System in South China

Mechanical hill direct seeding of hybrid rice could be the way to solve the problems of high seeding rates and uneven plant establishment now faced in direct seeded rice; however, it is not clear what the optimum hill seeding density should be for high-yielding hybrid rice in the single-season rice production system. Experiments were conducted in 2010 and 2011 to determine the effects of hill seeding density (25 cm 615 cm, 25 cm 617 cm, 25 cm 619 cm, 25 cm 621 cm, and 25 cm 623 cm; three to five seeds per hill) on plant growth and grain yield of a hybrid variety, Nei2you6, in two fields with different fertility (soil fertility 1 and 2). In addition, in 2012 and 2013, comparisons among mechanical hill seeding, broadcasting, and transplanting were conducted with three hybrid varieties to evaluate the optimum seeding density. With increases in seeding spacing from 25 cm615 cm to 25 cm623 cm, productive tillers per hill increased by 34.2% and 50.0% in soil fertility 1 and 2. Panicles per m2 declined with increases in seeding spacing in soil fertility 1. In soil fertility 2, no difference in panicles per m2 was found at spacing ranging from 25 cm617 cm to 25 cm623 cm, while decreases in the area of the top three leaves and aboveground dry weight per shoot at flowering were observed. Grain yield was the maximum at 25 cm 617 cm spacing in both soil fertility fields. Our results suggest that a seeding density of 25 cm617 cm was suitable for high-yielding hybrid rice. These results were verified through on-farm demonstration experiments, in which mechanical hill-seeded rice at this density had equal or higher grain yield than transplanted rice

Construction of microsatellite linkage maps for Corymbia

The genus Corymbia is closely related to the genus Eucalyptus, and like Eucalyptus contains tree species that are important for sub-tropical forestry. Corymbia's close relationship with Eucalyptus suggests genetic studies in Corymbia should benefit from transfer of genetic information from its more intensively studied relatives. Here we report a genetic map for Corymbia spp. based on microsatellite markers identified de novo in Corymbia sp or transferred from Eucalyptus. A framework consensus map was generated from an outbred F 2 population (n = 90) created by crossing two unrelated Corymbia torelliana x C. citriodora subsp. variegata F1 trees. The map had a total length of 367 cM (Kosambi) and was composed of 46 microsatellite markers distributed across 13 linkage groups (LOD 3). A high proportion of Eucalyptus microsatellites (90%) transferred to Corymbia. Comparative analysis between the Corymbia map and a published Eucalyptus map identified eight homeologous linkage groups in Corymbia with 13 markers mapping on one or both maps. Further comparative analysis was limited by low power to detect linkage due to low genome coverage in Corymbia, however, there was no convincing evidence for chromosomal structural differences because instances of non-synteny were associated with large distances on the Eucalyptus map. Segregation distortion was primarily restricted to a single linkage group and due to a deficit of hybrid genotypes, suggesting that hybrid inviability was one factor shaping the genetic composition of the F2 population in this inter-subgeneric hybrid. The conservation of microsatellite loci and synteny between Corymbia and Eucalyptus suggests there will be substantial value in exchanging information between the two groups.

CattleTickBase: An integrated Internet-based bioinformatics resource for Rhipicephalus (Boophilus) microplus

The Rhipicephalus micro plus genome is large and complex in structure, making it difficult to assemble a genome sequence and costly to resource the required bioinformatics. In light of this, a consortium of international collaborators was formed to pool resources to begin sequencing this genome. We have acquired and assembled genomic DNA into contigs that represent over 1.8 Gigabase pairs of DNA from gene-enriched regions of the R. micro plus genome. We also have several datasets containing transcript sequences from a number of gene expression experiments conducted by the consortium. A web-based resource was developed to enable the scientific community to access our datasets and conduct analysis through a web-based bioinformatics environment called YABI. The collective bioinformatics resource is termed CattleTickBase. Our consortium has acquired genomic and transcriptomic sequence data at approximately 0.9X coverage of the gene-coding regions of the R. microplus genome. The YABI tool will facilitate access and manipulation of cattle tick genome sequence data as the genome sequencing of R. microplus proceeds. During this process the CattleTickBase resource will continue to be updated. Published by Elsevier Ltd. on behalf of Australian Society for Parasitology Inc.

Hybridisation, paternal leakage and mitochondrial DNA linearization in three anomalous fish (Scombridae)

Using mitochondrial DNA for species identification and population studies assumes that the genome is maternally inherited, circular, located in the cytoplasm and lacks recombination. This study explores the mitochondrial genomes of three anomalous mackerel. Complete mitochondrial genome sequencing plus nuclear microsatellite genotyping of these fish identified them as Scomberomorus munroi (spotted mackerel). Unlike normal S. munroi, these three fish also contained different linear, mitochondrial genomes of Scomberomorus semifasciatus (grey mackerel). The results are best explained by hybridisation, paternal leakage and mitochondrial DNA linearization. This unusual observation may provide an explanation for mtDNA outliers in animal population studies. © 2013.

Development of single nucleotide polymorphism (SNP) markers from the mango (Mangifera indica) transcriptome for mapping and estimation of genetic diversity

The development of molecular markers for genomic studies in Mangifera indica (mango) will allow marker-assisted selection and identification of genetically diverse germplasm, greatly aiding mango breeding programs. We report here our identification of thousands of unambiguous molecular markers that can be easily assayed across genotypes of the species. With origin centered in Southeast Asia, mangos are grown throughout the tropics and subtropics as a nutritious fruit that exhibits remarkable intraspecific phenotypic diversity. With the goal of building a high density genetic map, we have undertaken discovery of sequence variation in expressed genes across a broad range of mango cultivars. A transcriptome sequence reference was built de novo from extensive sequencing and assembly of RNA from cultivar 'Tommy Atkins'. Single nucleotide polymorphisms (SNPs) in protein coding transcripts were determined from alignment of RNA reads from 24 mango cultivars of diverse origins: 'Amin Abrahimpur' (India), 'Aroemanis' (Indonesia), 'Burma' (Burma), 'CAC' (Hawaii), 'Duncan' (Florida), 'Edward' (Florida), 'Everbearing' (Florida), 'Gary' (Florida), 'Hodson' (Florida), 'Itamaraca' (Brazil), 'Jakarata' (Florida), 'Long' (Jamaica), 'M. Casturi Purple' (Borneo), 'Malindi' (Kenya), 'Mulgoba' (India), 'Neelum' (India), 'Peach' (unknown), 'Prieto' (Cuba), 'Sandersha' (India), 'Tete Nene' (Puerto Rico), 'Thai Everbearing' (Thailand), 'Toledo' (Cuba), 'Tommy Atkins' (Florida) and 'Turpentine' (West Indies). SNPs in a selected subset of protein coding transcripts are currently being converted into Fluidigm assays for genotyping of mapping populations and germplasm collections. Using an alternate approach, SNPs (144) discovered by sequencing of candidate genes in 'Kensington Pride' have already been converted and used for genotyping.

Biological management of the invasive Vachellia nilotica ssp. indica (Benth.) Kyal. & Boatwr. (Fabales: Mimosoideae) in tropical Australia: stress-inducing potential of Anomalococcus indicus Ramakrishna (Insecta: Hemiptera: Coccoidea: Lecanodiaspididae), an agent of promise

Vachellia nilotica ssp. indica (hereafter, V. n. indica) is an important tree weed in Australia. Its dense populations induce undesirable changes in the vast areas of northern Australia. Because chemical and mechanical management options appear unviable for various reasons, biological management of this tree is considered a better option. Among the many trialled arthropods in Australian context, Anomalococcus indicus, a lecanodiaspid native to India, has been identified as a potent-candidate, since in India, its native terrain, it is the most widespread and occurs throughout the year. Severe infestations of A. indicus cause defoliation, wilting and death of branches, and occasionally the tree. Populations of A. indicus have been brought into Australia and are being tested for its host specificity under quarantine conditions. This article reports the physiological damage and stress it inflicts in the shoots of V. n. indica. Younger-nymphal instars of A. indicus feed on cortical-parenchyma cells of young stems, whereas the older instars and adults feed from the phloem of old stems. Two conspicuous responses of V. n. indica arising in response to the feeding action of A. indicus are changes in the cell-wall dynamics and irregular cell divisions. The feeding action of A. indicus elicits a sequence of reactions in the stem tissues of V. n. indica such as differentiation of thick-walled elements in the outer cortical parenchyma, differential thickening of cells with supernumerary layers of either suberin or lignin, proliferations of parenchyma and phloem, wall thickening and obliteration of inner lumen of phloem cells, and the sieve plates plugged with callosic deposits. The responses are the culminations of interaction between the virulence factor (one or more of the salivary proteins?) from A. indicus and the resistance factor in V. n. indica. We have analysed structural changes in the context of their functions, by comparing the feeding action of A. indicus with that of other hemipteroids. From the level of stress it induces, this study confirms that A. indicus has the potential to be an effective biological management of V. n. indica in Australia. © 2014 © 2014 Taylor & Francis and Aboricultural Association.

Identification of putative virulence-associated genes of Haemophilus parasuis through suppression subtractive hybridization.

Haemophilus parasuis is the causative agent of Glässer's disease. Up to now 15 serovars of H. parasuis have been identified, with significant differences existing in virulence between serovars. In this study, suppression subtractive hybridization (SSH) was used to identify the genetic difference between Nagasaki (H. parasuis serovar 5 reference strain, highly virulent) and SW114 (H. parasuis serovar 3 reference strain, non-virulent). A total of 191 clones were obtained from the SSH library. Using dot hybridization and PCR, 15 clones were identified containing fragments that were present in the Nagasaki genome while absent in the SW114 genome. Among these 15 fragments, three fragments (ssh1, ssh13, ssh15) encode cell surface-associated components; three fragments (ssh2, ssh5, ssh9) are associated with metabolism and stress response; one fragment (ssh8) is involved in assembly of fimbria and one fragment (ssh6) is a phage phi-105 ORF25-like protein. The remaining seven fragments are hypothetical proteins or unknown. Based on PCR analysis of the 15 serovar reference strains, eight fragments (ssh1, ssh2, ssh3, ssh6, ssh8, ssh10, ssh11 and ssh12) were found in three to five of most virulent serovars (1, 5, 10, 12, 13 and 14), zero to two in three moderately virulent serovars (2, 4 and 15), but absent in the low virulent serovar (8) and non-virulent serovars (3, 6, 7, 9 and 11). In vivo transcription fragments ssh1, ssh2, ssh8 and ssh12 were identified in total RNA samples extracted from experimental infected pig lung by RT-PCR. This study has provided some evidence of genetic differences between H. parasuis strains of different virulence.

Drafting the Kensington Pride Mango Genome

The mango industry in Australia is worth in excess of $150 million annually with the Kensington Pride (KP) cultivar capturing 60% of the domestic market. Valued by consumers for desirable taste and colour characteristics, KP has been used extensively as a parent in the Department of Agriculture and Fisheries’ (Queensland, Australia) mango breeding program with over 400 hybrid trees sharing KP as the male parent. In order to gain a better understanding of Australia’s most significant mango variety, Horticulture Innovation Australia had led an international collaboration between the Queensland Department of Agriculture and Fisheries (Australia), the International Crops Research Institute for the Semi-Arid Tropics (ICRISAT, India) and the Beijing Genomics Institute (China) to sequence the KP genome. Preliminary de novo assembly of illumina short read sequence data suggests that the KP genome is highly heterozygous and has an estimated genome size of 407 Mb. As refinements and additional sequence data are added to the assembly, a more complete picture of the mango genome will be elucidated.

Pasteurella multocida Expresses Two Lipopolysaccharide Glycoforms Simultaneously, but Only a Single Form Is Required for Virulence: Identification of Two Acceptor-Specific Heptosyl I Transferases

Lipopolysaccharide (LPS) is a critical virulence determinant in Pasteurella multocida and a major antigen responsible for host protective immunity. In other mucosal pathogens, variation in LPS or lipooligosaccharide structure typically occurs in the outer core oligosaccharide regions due to phase variation. P. multocida elaborates a conserved oligosaccharide extension attached to two different, simultaneously expressed inner core structures, one containing a single phosphorylated 3-deoxy-D-manno-octulosonic acid (Kdo) residue and the other containing two Kdo residues. We demonstrate that two heptosyltransferases, HptA and HptB, add the first heptose molecule to the Kdo1 residue and that each exclusively recognizes different acceptor molecules. HptA is specific for the glycoform containing a single, phosphorylated Kdo residue (glycoform A), while HptB is specific for the glycoform containing two Kdo residues (glycoform B). In addition, KdkA was identified as a Kdo kinase, required for phosphorylation of the first Kdo molecule. Importantly, virulence data obtained from infected chickens showed that while wild-type P. multocida expresses both LPS glycoforms in vivo, bacterial mutants that produced only glycoform B were fully virulent, demonstrating for the first time that expression of a single LPS form is sufficient for P. multocida survival in vivo. We conclude that the ability of P. multocida to elaborate alternative inner core LPS structures is due to the simultaneous expression of two different heptosyltransferases that add the first heptose residue to the nascent LPS molecule and to the expression of both a bifunctional Kdo transferase and a Kdo kinase, which results in the initial assembly of two inner core structures.

Decoration of Pasteurella multocida lipopolysaccharide with phosphocholine is important for virulence

Phosphocholine (PCho) is an important substituent of surface structures expressed by a number of bacterial pathogens. Its role in virulence has been investigated in several species, in which it has been shown to play a role in bacterial adhesion to mucosal surfaces, in resistance to antimicrobial peptides, or in sensitivity to complement-mediated killing. The lipopolysaccharide (LPS) structure of Pasteurella multocida strain Pm70, whose genome sequence is known, has recently been determined and does not contain PCho. However, LPS structures from the closely related, virulent P. multocida strains VP161 and X-73 were shown to contain PCho on their terminal galactose sugar residues. To determine if PCho was involved in the virulence of P. multocida, we used subtractive hybridization of the VP161 genome against the Pm70 genome to identify a four-gene locus (designated pcgDABC) which we show is required for the addition of the PCho residues to LPS. The proteins predicted to be encoded by pcgABC showed identity to proteins involved in choline uptake, phosphorylation, and nucleotide sugar activation of PCho. We constructed a P. multocida VP161 pcgC mutant and demonstrated that this strain produces LPS that lacks PCho on the terminal galactose residues. This pcgC mutant displayed reduced in vivo growth in a chicken infection model and was more sensitive to the chicken antimicrobial peptide fowlicidin-1 than the wild-type P. multocida strain

The complete validated mitochondrial genome of the silver gemfish Rexea solandri (Cuvier, 1832) (Perciformes, Gempylidae)

The silver gemfish Rexea solandri is an important economic resource but vulnerable to overfishing in Australian waters. The complete mitochondrial genome sequence is described from 1.6 million reads obtained via next generation sequencing. The total length of the mitogenome is 16,350 bp comprising 2 rRNA, 13 protein-coding genes, 22 tRNA and 2 non-coding regions. The mitogenome sequence was validated against sequences of PCR fragments and BLAST queries of Genbank. Gene order was equivalent to that found in marine fishes.

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.