Analyses of Twelve New Whole Genome Sequences of Cassava Brown Streak Viruses and Ugandan Cassava Brown Streak Viruses from East Africa: Diversity,Supercomputing and Evidence for Further Speciation

Cassava brown streak disease is caused by two devastating viruses, Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV) which are frequently found infecting cassava, one of sub-Saharan Africa’s most important staple food crops. Each year these viruses cause losses of up to $100 million USD and can leave entire families without their primary food source, for an entire year. Twelve new whole genomes, including seven of CBSV and five of UCBSV were uncovered in this research, doubling the genomic sequences available in the public domain for these viruses. These new sequences disprove the assumption that the viruses are limited by agro-ecological zones, show that current diagnostic primers are insufficient to provide confident diagnosis of these viruses and give rise to the possibility that there may be as many as four distinct species of virus. Utilizing NGS sequencing technologies and proper phylogenetic practices will rapidly increase the solution to sustainable cassava production.     

Diversity,Distribution and Effects on Cassava Cultivars of Cassava Brown Streak Viruses in Malawi

Cassava brown streak disease (CBSD) has emerged as a major threat to cassava (Manihot esculenta) in eastern and southern Africa. CBSD was first reported in Malawi in the 1950s, but little data on the distribution and epidemiology of the disease are available. A diagnostic survey was therefore conducted in Malawi to determine the distribution, incidence and diversity of viruses causing the disease, and to characterize its effects on local cassava cultivars. Diagnostic tests confirmed the presence of cassava brown streak viruses (CBSVs) in 90% of leaf samples from symptomatic plants. Average CBSD foliar severity was 2.5, although this varied significantly between districts. Both Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV) (genus Ipomovirus, family Potyviridae) were detected from sampled plants. UCBSV was widespread, whereas CBSV was detected only in the two most northerly districts. The average abundance of the whitefly vector (Bemisia tabaci) was 0.4 per plant, a low value that was partly attributable to the fact that the survey was conducted during the cool part of the year known to be unfavourable for B. tabaci whiteflies. Spearman's correlation analyses showed a positive correlation between CBSD foliar incidence and CBSD severity and between CBSD severity and CBSD stem incidence. Of the 31 cassava varieties encountered, 20–20 was most severely affected, whilst Mtutumusi was completely unaffected. Although data from this study do not indicate a significant CBSD deterioration in Malawi, strengthened management efforts are required to reduce the current impact of the disease.     

Exploiting the Combination of Natural and Genetically Engineered Resistance to Cassava Mosaic and Cassava Brown Streak Viruses Impacting Cassava Production in Africa

Cassava brown streak disease (CBSD) and cassava mosaic disease (CMD) are currently two major viral diseases that severely reduce cassava production in large areas of Sub-Saharan Africa. Natural resistance has so far only been reported for CMD in cassava. CBSD is caused by two virus species, Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV). A sequence of the CBSV coat protein (CP) highly conserved between the two virus species was used to demonstrate that a CBSV-CP hairpin construct sufficed to generate immunity against both viral species in the cassava model cultivar (cv. 60444). Most of the transgenic lines showed high levels of resistance under increasing viral loads using a stringent top-grafting method of inoculation. No viral replication was observed in the resistant transgenic lines and they remained free of typical CBSD root symptoms 7 month post-infection. To generate transgenic cassava lines combining resistance to both CBSD and CMD the hairpin construct was transferred to a CMD-resistant farmer-preferred Nigerian landrace TME 7 (Oko-Iyawo). An adapted protocol allowed the efficient Agrobacterium-based transformation of TME 7 and the regeneration of transgenic lines with high levels of CBSV-CP hairpin-derived small RNAs. All transgenic TME 7 lines were immune to both CBSV and UCBSV infections. Further evaluation of the transgenic TME 7 lines revealed that CBSD resistance was maintained when plants were co-inoculated with East African cassava mosaic virus (EACMV), a geminivirus causing CMD. The innovative combination of natural and engineered virus resistance in farmer-preferred landraces will be particularly important to reducing the increasing impact of cassava viral diseases in Africa.     

A Modified RNA-Seq Approach for Whole Genome Sequencing of RNA Viruses from Faecal and Blood Samples

To date, very large scale sequencing of many clinically important RNA viruses has been complicated by their high population molecular variation, which creates challenges for polymerase chain reaction and sequencing primer design. Many RNA viruses are also difficult or currently not possible to culture, severely limiting the amount and purity of available starting material. Here, we describe a simple, novel, high-throughput approach to Norovirus and Hepatitis C virus whole genome sequence determination based on RNA shotgun sequencing (also known as RNA-Seq). We demonstrate the effectiveness of this method by sequencing three Norovirus samples from faeces and two Hepatitis C virus samples from blood, on an Illumina MiSeq benchtop sequencer. More than 97% of reference genomes were recovered. Compared with Sanger sequencing, our method had no nucleotide differences in 14,019 nucleotides (nt) for Noroviruses (from a total of 2 Norovirus genomes obtained with Sanger sequencing), and 8 variants in 9,542 nt for Hepatitis C virus (1 variant per 1,193 nt). The three Norovirus samples had 2, 3, and 2 distinct positions called as heterozygous, while the two Hepatitis C virus samples had 117 and 131 positions called as heterozygous. To confirm that our sample and library preparation could be scaled to true high-throughput, we prepared and sequenced an additional 77 Norovirus samples in a single batch on an Illumina HiSeq 2000 sequencer, recovering >90% of the reference genome in all but one sample. No discrepancies were observed across 118,757 nt compared between Sanger and our custom RNA-Seq method in 16 samples. By generating viral genomic sequences that are not biased by primer-specific amplification or enrichment, this method offers the prospect of large-scale, affordable studies of RNA viruses which could be adapted to routine diagnostic laboratory workflows in the near future, with the potential to directly characterize within-host viral diversity.     

The Genomics of Speciation in Drosophila: Diversity, Divergence, and Introgression Estimated Using Low-Coverage Genome Sequencing

In nature, closely related species may hybridize while still retaining their distinctive identities. Chromosomal regions that experience reduced recombination in hybrids, such as within inversions, have been hypothesized to contribute to the maintenance of species integrity. Here, we examine genomic sequences from closely related fruit fly taxa of the Drosophila pseudoobscura subgroup to reconstruct their evolutionary histories and past patterns of genic exchange. Partial genomic assemblies were generated from two subspecies of Drosophila pseudoobscura (D. ps.) and an outgroup species, D. miranda. These new assemblies were compared to available assemblies of D. ps. pseudoobscura and D. persimilis, two species with overlapping ranges in western North America. Within inverted regions, nucleotide divergence among each pair of the three species is comparable, whereas divergence between D. ps. pseudoobscura and D. persimilis in non-inverted regions is much lower and closer to levels of intraspecific variation. Using molecular markers flanking each of the major chromosomal inversions, we identify strong crossover suppression in F₁ hybrids extending over 2 megabase pairs (Mbp) beyond the inversion breakpoints. These regions of crossover suppression also exhibit the high nucleotide divergence associated with inverted regions. Finally, by comparison to a geographically isolated subspecies, D. ps. bogotana, our results suggest that autosomal gene exchange between the North American species, D. ps. pseudoobscura and D. persimilis, occurred since the split of the subspecies, likely within the last 200,000 years. We conclude that chromosomal rearrangements have been vital to the ongoing persistence of these species despite recent hybridization. Our study serves as a proof-of-principle on how whole genome sequencing can be applied to formulate and test hypotheses about species formation in lesser-known non-model systems.     

Characterization of Foot-And-Mouth Disease Viruses (FMDVs) from Ugandan Cattle Outbreaks during 2012-2013: Evidence for Circulation of Multiple Serotypes

To investigate the foot-and-mouth disease virus (FMDV) serotypes circulating in Uganda’s cattle population, both serological and virological analyses of samples from outbreaks that occurred during 2012–2013 were performed. Altogether, 79 sera and 60 oropharyngeal fluid (OP)/tissue/oral swab samples were collected from herds with reported FMD outbreaks in seven different Ugandan districts. Overall, 61/79 (77%) of the cattle sera were positive for antibodies against FMDV by PrioCHECK FMDV NS ELISA and solid phase blocking ELISA detected titres ≥ 80 for serotypes O, SAT 1, SAT 2 and SAT 3 in 41, 45, 30 and 45 of these 61 seropositive samples, respectively. Virus neutralisation tests detected the highest levels of neutralising antibodies (titres ≥ 45) against serotype O in the herds from Kween and Rakai districts, against SAT 1 in the herd from Nwoya district and against SAT 2 in the herds from Kiruhura, Isingiro and Ntungamo districts. The isolation of a SAT 2 FMDV from Isingiro was consistent with the detection of high levels of neutralising antibodies against SAT 2; sequencing (for the VP1 coding region) indicated that this virus belonged to lineage I within this serotype, like the currently used vaccine strain. From the Wakiso district 11 tissue/swab samples were collected; serotype A FMDV, genotype Africa (G-I), was isolated from the epithelial samples. This study shows that within a period of less than one year, FMD outbreaks in Uganda were caused by four different serotypes namely O, A, SAT 1 and SAT 2. Therefore, to enhance the control of FMD in Uganda, there is need for efficient and timely determination of outbreak virus strains/serotypes and vaccine matching. The value of incorporating serotype A antigen into the imported vaccines along with the current serotype O, SAT 1 and SAT 2 strains should be considered.     

Plasmid Flux in Escherichia coli ST131 Sublineages,Analyzed by Plasmid Constellation Network (PLACNET), a New Method for Plasmid Reconstruction from Whole Genome Sequences

Bacterial whole genome sequence (WGS) methods are rapidly overtaking classical sequence analysis. Many bacterial sequencing projects focus on mobilome changes, since macroevolutionary events, such as the acquisition or loss of mobile genetic elements, mainly plasmids, play essential roles in adaptive evolution. Existing WGS analysis protocols do not assort contigs between plasmids and the main chromosome, thus hampering full analysis of plasmid sequences. We developed a method (called plasmid constellation networks or PLACNET) that identifies, visualizes and analyzes plasmids in WGS projects by creating a network of contig interactions, thus allowing comprehensive plasmid analysis within WGS datasets. The workflow of the method is based on three types of data: assembly information (including scaffold links and coverage), comparison to reference sequences and plasmid-diagnostic sequence features. The resulting network is pruned by expert analysis, to eliminate confounding data, and implemented in a Cytoscape-based graphic representation. To demonstrate PLACNET sensitivity and efficacy, the plasmidome of the Escherichia coli lineage ST131 was analyzed. ST131 is a globally spread clonal group of extraintestinal pathogenic E. coli (ExPEC), comprising different sublineages with ability to acquire and spread antibiotic resistance and virulence genes via plasmids. Results show that plasmids flux in the evolution of this lineage, which is wide open for plasmid exchange. MOB_F12/IncF plasmids were pervasive, adding just by themselves more than 350 protein families to the ST131 pangenome. Nearly 50% of the most frequent γ–proteobacterial plasmid groups were found to be present in our limited sample of ten analyzed ST131 genomes, which represent the main ST131 sublineages.     

Two New Species of Liuixalus (Rhacophoridae,Anura): Evidence from Morphological and Molecular Analyses

Due to small body sizes, superficial similarities in morphologies, and obscure activity behaviors, the phylogeny and taxonomy of species in the genus Liuixalus were very troublesome. Some species might comprise a complex of cryptic species. To investigate the species of group, we constructed the matrilineal genealogy of the genus using 16s rRNA mitochondrial DNA sequences. Analyses recovered six well supported matrilines that involved L. romeri, L. ocellatus, L. hainanus, L. calcarius, Liuixalus shiwandashansp. nov. and Liuixalus jinxiuensissp. nov., though the historical relationships among them remained unresolved. Currently, Liuixalus included 4 species, distributed eastwards from northern Vietnam to Hong Kong, China. Based on genealogical and morphological distinctiveness, we described Liuixalus jinxiuensissp. nov. from the type locality Mt. Dayao, Jinxiu, China and Liuixalus shiwandashansp. nov. from the type locality Mt. Shiwanda, China. A combination of morphological measurements, genetic, bioacoustic and osteological analysis was served to diagnose the new taxa.     

Molecular Analysis of H7 Avian Influenza Viruses from Australia and New Zealand: Genetic Diversity and Relationships from 1976 to 2007

Full-genome sequencing of 11 Australian and 1 New Zealand avian influenza A virus isolate (all subtype H7) has enabled comparison of the sequences of each of the genome segments to those of other subtype H7 avian influenza A viruses. The inference of phylogenetic relationships for each segment has been used to develop a model of the natural history of these viruses in Australia. Phylogenetic analysis of the hemagglutinin segment indicates that the Australian H7 isolates form a monophyletic clade. This pattern is consistent with the long-term, independent evolution that is, in this instance, associated with geographic regions. On the basis of the analysis of the other H7 hemagglutinin sequences, three other geographic regions for which similar monophyletic clades have been observed were confirmed. These regions are Eurasia plus Africa, North America, and South America. Analysis of the neuraminidase sequences from the H7N1, H7N3, and H7N7 genomes revealed the same region-based relationships. This pattern of independent evolution of Australian isolates is supported by the results of analysis of each of the six remaining genomic segments. These results, in conjunction with the occurrence of five different combinations of neuraminidase subtypes (H7N2, H7N3, H7N4, H7N6, H7N7) among the 11 Australian isolates, suggest that the maintenance host(s) is nearly exclusively associated with Australia. The single lineage of Australian H7 hemagglutinin sequences, despite the occurrence of multiple neuraminidase types, suggests the existence of a genetic pool from which a variety of reassortants arise rather than the presence of a small number of stable viral clones. This pattern of evolution is likely to occur in each of the regions mentioned above.The emergence of highly pathogenic avian influenza viruses of subtype H5N1 as a potential human pandemic disease threat has focused attention on the roles that wild birds play in the maintenance and distribution of avian influenza viruses (18, 22). Moreover, the H5 and H7 subtypes of avian influenza A virus are major causes of economic loss in poultry production through disease. In Australia, there have been five documented outbreaks of H7 subtype avian influenza A virus disease, with evidence of adaptation to the poultry host being provided by sequence data supporting the presence of high-pathogenicity avian influenza virus (HPAI) isolates in poultry. Waterfowl (Anseriformes order, particularly ducks, geese, and swans) and the waders and gulls (Charadriiformes order, particularly gulls, terns, and waders) have been found to be the major global natural reservoirs of influenza A viruses. Transmission of avian influenza viruses from wild birds to production poultry and geographic spread are dependent upon the migratory behavior of the wild bird reservoir hosts. Members of the Anseriformes and Charadriiformes orders undertake both irregular and regular transcontinental and intercontinental migrations. During these migrations, large numbers of birds congregate at aquatic feeding locations, providing ideal sites for cross-species transmission of avian influenza viruses. A variety of mechanisms have been observed whereby influenza A viruses adapt rapidly. These include genetic shifts facilitated through genome segment reassortment, as well as genetic drift through the insertion, deletion, and substitution of nucleotides. The error-prone RNA replication and a lack of error correction are the causes of drift. In vivo, this results in viral genetic diversity within any viral sample, or a quasispecies, thus providing a pool of closely related variant viruses from enabling events, such as viral adaptation to new hosts (25). Long-term sampling of water birds in North America and Europe has started to elucidate the ecology and biology of the avian influenza A virus types in the natural reservoirs in these regions (8, 18, 22). There is a suggestion that two superfamilies, the Eurasia (which in the context of this paper includes Europe, Asia, and Africa) and the Americas superfamilies, exist; however, the extent of overlap and the rate of transfer of influenza viruses between these two regions are not well-defined. Recent studies suggest that intercontinental virus exchange is slow and limited (17), while a detailed analysis of the differences between H7 hemagglutinin (HA) segments circulating in Europe and China showed that the H7 hemagglutinin segments shared a recent common ancestor and limited sequence divergence on a background of multiple reassortant virus genotypes between 1999 and 2005 (7).Avian influenza A viruses of the Oceania region (Australia, New Zealand, and southwest Pacific) have been far less well studied (3). Australia and New Zealand are at the southern extremity of a number of major bird migration pathways. Waders in the Charadriidae family migrate to south and southeastern Australia and New Zealand from their summer breeding grounds in Arctic regions of Siberia and Alaska, where they freely mix with the same or other species which migrate into the shared breeding grounds of Eurasia and the Americas (30). Pelagic seabirds of the Procellariformes order breed on and around Australian and New Zealand coasts during the southern hemisphere summer and migrate to maritime regions of the northern Pacific associated with Japan, Russia, and Alaska. Some move as far as the west coasts of North and South America (28). Unlike North and South America and Europe, where regular migrations of ducks, geese, swans, etc., are established, the members of the Anatidae family (ducks, etc.) in Australia and New Zealand are mainly endemic residents (30). However, within Australia, ducks undertake long-distance movements in response to water availability. Movements of waterfowl from northern Australia to nearby areas of Southeast Asia are believed to occur but are limited, as suggested by Wallace''s Line (19). Generally, these waterfowl movements have not been well studied (30). The risks to Australian poultry production systems by movement of H5N1 via migratory shorebirds and nomadic wildfowl have been assessed to be low using risk-based analysis techniques (9, 10).Regular and extensive surveillance sampling of migratory birds has been undertaken in North America and northern Europe (17, 18). The findings have shed significant insights into the ecology of the viruses and their hosts (8, 17). In contrast, surveillance sampling of wild birds in Asia and Oceania has been spasmodic and sparse, until the recent emergence of H5N1 highly pathogenic avian influenza virus as a poultry and human disease threat. Spasmodic and small-scale outbreaks of highly pathogenic avian influenza virus have occurred in Australian poultry production flocks located in the southeastern region of the continent. These poultry production areas are concentrated close to large human population centers (26, 33, 34). Each of the Australian outbreaks has been rapidly controlled by slaughter of infected flocks. All have been caused by avian influenza viruses of the H7 subtype, which appear to have entered production poultry from water birds, possibly wild ducks, via contaminated water supplies used on the poultry farms. Disease has occurred on five occasions: 1976 (H7N7), 1985 (H7N7), 1992 (H7N3), 1994 (H7N3), and 1997 (H7N4) (13, 14, 26, 27, 31, 34). National on-farm biosecurity measures have been focused on reducing the likelihood of future outbreaks. The availability of avian influenza virus isolates from poultry and wild birds associated with these outbreaks, along with a small number of subtype H7 avian influenza viruses isolated from wild ducks during recent national surveillance programs in Australia and New Zealand, provided the opportunity to explore the relationships of Australian and New Zealand subtype H7 avian influenza virus isolates with viruses circulating elsewhere in the world.     

A New Tool for the Rapid Cloning of Amplified and Hypermethylated Human DNA Sequences from Restriction Landmark Genome Scanning Gels

Restriction landmark genome scanning (RLGS) is an effective genome-scanning technique capable of identifying DNA amplification and aberrant DNA methylation. Previously published methods for the cloning of human DNA fragments from RLGS gels have been successful only for high-copy-number fragments (repetitive elements or DNA amplifications). We present here the first technique capable of efficiently cloning single-copy human DNA fragments ("spots") identified in RLGS profiles. This technique takes advantage of a plasmid-based, human genomic DNA, NotI/EcoRV boundary library. The library is arrayed in microtiter plates. When clones from a single plate are pooled and mixed with genomic DNA, the resultant RLGS gel is a normal profile with a defined set of spots showing enhanced intensity for that particular plate. This was performed for a set of 32 plates as well as their pooled rows and columns. Thus, we have mapped individual RLGS spots to exact plate, row, and column addresses in the library and have thereby obtained immediate access to these clones. The feasibility of the technique is demonstrated in examples of cloning methylated DNA fragments identified in human breast tumor and testicular tumor RLGS profiles and in the cloning of an amplified DNA fragment identified in a human medulloblastoma RLGS profile.     

Processes Underpinning Development and Maintenance of Diversity in Rice in West Africa: Evidence from Combining Morphological and Molecular Markers

We assessed the interplay of artificial and natural selection in rice adaptation in low-input farming systems in West Africa. Using 20 morphological traits and 176 molecular markers, 182 farmer varieties of rice (Oryza spp.) from 6 West African countries were characterized. Principal component analysis showed that the four botanical groups (Oryza sativa ssp. indica, O. sativa ssp. japonica, O. glaberrima, and interspecific farmer hybrids) exhibited different patterns of morphological diversity. Regarding O. glaberrima, morphological and molecular data were in greater conformity than for the other botanical groups. A clear difference in morphological features was observed between O. glaberrima rices from the Togo hills and those from the Upper Guinea Coast, and among O. glaberrima rices from the Upper Guinea Coast. For the other three groups such clear patterns were not observed. We argue that this is because genetic diversity is shaped by different environmental and socio-cultural selection pressures. For O. glaberrima, recent socio-cultural selection pressures seemed to restrict genetic diversity while this was not observed for the other botanical groups. We also show that O. glaberrima still plays an important role in the selection practices of farmers and resulting variety development pathways. This is particularly apparent in the case of interspecific farmer hybrids where a relationship was found between pericarp colour, panicle attitude and genetic diversity. Farmer varieties are the product of long and complex trajectories of selection governed by local human agency. In effect, rice varieties have emerged that are adapted to West African farming conditions through genotype × environment × society interactions. The diversity farmers maintain in their rice varieties is understood to be part of a risk-spreading strategy that also facilitates successful and often serendipitous variety innovations. We advocate, therefore, that farmers and farmer varieties should be more effectively involved in crop development.     

New Findings on Evolution of Metal Homeostasis Genes: Evidence from Comparative Genome Analysis of Bacteria and Archaea

In order to examine the natural history of metal homeostasis genes in prokaryotes, open reading frames with homology to characterized P_IB-type ATPases from the genomes of 188 bacteria and 22 archaea were investigated. Major findings were as follows. First, a high diversity in N-terminal metal binding motifs was observed. These motifs were distributed throughout bacterial and archaeal lineages, suggesting multiple loss and acquisition events. Second, the CopA locus separated into two distinct phylogenetic clusters, CopA1, which contained ATPases with documented Cu(I) influx activity, and CopA2, which contained both efflux and influx transporters and spanned the entire diversity of the bacterial domain, suggesting that CopA2 is the ancestral locus. Finally, phylogentic incongruences between 16S rRNA and P_IB-type ATPase gene trees identified at least 14 instances of lateral gene transfer (LGT) that had occurred among diverse microbes. Results from bootstrapped supported nodes indicated that (i) a majority of the transfers occurred among proteobacteria, most likely due to the phylogenetic relatedness of these organisms, and (ii) gram-positive bacteria with low moles percent G+C were often involved in instances of LGT. These results, together with our earlier work on the occurrence of LGT in subsurface bacteria (J. M. Coombs and T. Barkay, Appl. Environ. Microbiol. 70:1698-1707, 2004), indicate that LGT has had a minor role in the evolution of P_IB-type ATPases, unlike other genes that specify survival in metal-stressed environments. This study demonstrates how examination of a specific locus across microbial genomes can contribute to the understanding of phenotypes that are critical to the interactions of microbes with their environment.     

A New Marine Reptile (Sauropterygia) from New Zealand: Further Evidence for A Late Cretaceous Austral Radiation of Cryptoclidid Plesiosaurs

Kaiwhekea katiki gen. et sp. nov. represents the first described cryptoclidid plesiosaurian from New Zealand. It is one of the largest cryptoclidids known, at a length of over 6.5 m, and represents the third reported genus of austral Late Cretaceous cryptoclidids. Kaiwhekea katiki is from siltstones of the Katiki Formation, upper Haumurian Stage (Cenomanian–Maastrichtian; c. 69–70 Ma) of coastal Otago, South Island, New Zealand. In the Late Cretaceous, the locality lay close to the polar circle. The holotype and only known specimen is an articulated skeleton with skull, preserved mostly as natural molds, but which lacks the forelimbs and pectoral girdle. The skull is relatively large and possesses several distinct characters, including a substantial, deep, jugal. There are about 43 upper and 42 lower teeth in each jaw quadrant; all are homodont, slim, and slightly recurved, lacking prominent ornament. Kaiwhekea probably took single soft-bodied prey. Based on cranial structure, it clearly belongs with the Cryptoclididae, but is not certainly close to the southern Late Cretaceous cryptoclidids Morturneria (Seymour Island, Antarctica) and Aristonectes (Chile, Argentina).