When one phenotype is not enough: divergent evolutionary trajectories govern venom variation in a widespread rattlesnake species.

Understanding the origin and maintenance of phenotypic variation, particularly across a continuous spatial distribution, represents a key challenge in evolutionary biology. For this, animal venoms represent ideal study systems: they are complex, variable, yet easily quantifiable molecular phenotypes with a clear function. Rattlesnakes display tremendous variation in their venom composition, mostly through strongly dichotomous venom strategies, which may even coexist within a single species. Here, through dense, widespread population-level sampling of the Mojave rattlesnake, Crotalus scutulatus, we show that genomic structural variation at multiple loci underlies extreme geographical variation in venom composition, which is maintained despite extensive gene flow. Unexpectedly, neither diet composition nor neutral population structure explain venom variation. Instead, venom divergence is strongly correlated with environmental conditions. Individual toxin genes correlate with distinct environmental factors, suggesting that different selective pressures can act on individual loci independently of their co-expression patterns or genomic proximity. Our results challenge common assumptions about diet composition as the key selective driver of snake venom evolution and emphasize how the interplay between genomic architecture and local-scale spatial heterogeneity in selective pressures may facilitate the retention of adaptive functional polymorphisms across a continuous space.

The proposed new species, cacao red vein virus, and three previously recognized badnavirus species are associated with cacao swollen shoot disease.

BACKGROUND: Cacao swollen shoot virus (CSSV), Cacao swollen shoot CD virus (CSSCDV), and Cacao swollen shoot Togo A virus (CSSTAV) cause cacao swollen shoot disease (CSSD) in West Africa. During 2000-2003, leaf and shoot-swelling symptoms and rapid tree death were observed in cacao in Cote d'Ivoire and Ghana. Molecular tests showed positive infection in only ~50-60% of symptomatic trees, suggesting the possible emergence of an unknown badnavirus. METHODS: The DNA virome was determined from symptomatic cacao samples using Illumina-Hi Seq, and sequence accuracy was verified by Sanger sequencing. The resultant 14, and seven previously known, full-length badnaviral genomic and RT-RNase H sequences were analyzed by pairwise distance analysis to resolve species relationships, and by Maximum likelihood (ML) to reconstruct phylogenetic relationships. The viral coding and non-coding sequences, genome organization, and predicted conserved protein domains (CPDs) were identified and characterized at the species level. RESULTS: The 21 CSSD-badnaviral genomes and RT-RNase H sequences shared 70-100% and 72-100% identity, respectively. The RT-RNase H analysis predicted four species, based on an ≥80% species cutoff. The ML genome sequence tree resolved three well-supported clades, with ≥70% bootstrap, whereas, the RT-RNase H phylogeny was poorly resolved, however, both trees grouped CSSD isolates within one large clade, including the newly discovered Cacao red vein virus (CRVV) proposed species. The genome arrangement of the four species consists of four, five, or six predicted open reading frames (ORFs), and the CPDs have similar architectures. By comparison, two New World cacao-infecting badnaviruses encode four ORFs, and harbor CPDs like the West African species. CONCLUSIONS: Three previously recognized West African cacao-infecting badnaviral species were identified, and a fourth, previously unidentified species, CRVV, is described for the first time. The CRVV is a suspect causal agent of the rapid decline phenotype, however Koch's Postulates have not been proven. To reconcile viral evolutionary with epidemiology considerations, more detailed information about CSSD-genomic variability is essential. Also, the functional basis for the multiple genome arrangements and subtly distinct CPD architectures among cacao-infecting badnaviruses is poorly understood. New knowledge about functional relationships may help explain the diverse symptomatologies observed in affected cacao trees.

Incipient speciation with biased gene flow between two lineages of the Western Diamondback Rattlesnake (Crotalus atrox).

We used mitochondrial DNA sequence data from 151 individuals to estimate population genetic structure across the range of the Western Diamondback Rattlesnake (Crotalus atrox), a widely distributed North American pit viper. We also tested hypotheses of population structure using double-digest restriction site associated DNA (ddRADseq) data, incorporating thousands of nuclear genome-wide SNPs from 42 individuals. We found strong mitochondrial support for a deep divergence between eastern and western C. atrox populations, and subsequent intermixing of these populations in the Inter-Pecos region of the United States and Mexico. Our nuclear RADseq data also identify these two distinct lineages of C. atrox, and provide evidence for nuclear admixture of eastern and western alleles across a broad geographic region. We identified contrasting patterns of mitochondrial and nuclear genetic variation across this genetic fusion zone that indicate partially restricted patterns of gene flow, which may be due to either pre- or post-zygotic isolating mechanisms. The failure of these two lineages to maintain complete genetic isolation, and evidence for partially-restricted gene flow, imply that these lineages were in the early stages of speciation prior to secondary contact.

Using rattlesnake microsatellites to determine paternity in captive bushmasters (Lachesis muta).

We used 10 microsatellite DNA markers originally described for two Crotalus and one Sistrurus species to infer paternity in a captive-hatched clutch of Lachesis muta. Although the dam was known, records listed two potential sires, which prevented the inclusion of those offspring in a captive breeding program. Samples were collected from both possible sires, the presumed dam and two offspring. Five of the ten markers were paternity informative and either unambiguously identified one of the two males as the sire of the offspring or excluded the other male. A sixth marker identified the sire in one of the two offspring. These results reveal rattlesnake DNA markers as a useful tool in paternity testing and captive breeding management of Lachesis muta. It also indicates that, within the American crotalid radiation, a selection of microsatellite DNA markers cross-amplify and provide useful genotypic information across species and genera.

Sentinel plot surveillance of cotton leaf curl disease in Pakistan- a case study at the cultivated cotton-wild host plant interface.

A sentinel plot case study was carried out to identify and map the distribution of begomovirus-betasatellite complexes in sentinel plots and commercial cotton fields over a four-year period using molecular and high-throughput DNA 'discovery' sequencing approaches. Samples were collected from 15 study sites in the two major cotton-producing areas of Pakistan. Whitefly- and leafhopper-transmitted geminiviruses were detected in previously unreported host plant species and locations. The most prevalent begomovirus was cotton leaf curl Kokhran virus-Burewala (CLCuKoV-Bu). Unexpectedly, a recently recognized recombinant, cotton leaf curl Multan virus-Rajasthan (CLCuMuV-Ra) was prevalent in five of 15 sites. cotton leaf curl Alabad virus (CLCuAlV) and cotton leaf curl Kokhran virus-Kokhran, 'core' members of CLCuD-begomoviruses that co-occurred with CLCuMuV in the 'Multan' epidemic were detected in one of 15 sentinel plots. Also identified were chickpea chlorotic dwarf virus and 'non-core' CLCuD-begomoviruses, okra enation leaf curl virus, squash leaf curl virus, and tomato leaf curl New Delhi virus. Cotton leaf curl Multan betasatellite (CLCuMuB) was the most prevalent CLCuD-betasatellite, and less commonly, two 'non-core' betasatellites. Recombination analysis revealed previously uncharacterized recombinants among helper virus-betasatellite complexes consisting of CLCuKoV, CLCuMuV, CLCuAlV and CLCuMuB. Population analyses provided early evidence for CLCuMuV-Ra expansion and displacement of CLCuKoV-Bu in India and Pakistan from 2012-2017. Identification of 'core' and non-core CLCuD-species/strains in cotton and other potential reservoirs, and presence of the now predominant CLCuMuV-Ra strain are indicative of ongoing diversification. Investigating the phylodynamics of geminivirus emergence in cotton-vegetable cropping systems offers an opportunity to understand the driving forces underlying disease outbreaks and reconcile viral evolution with epidemiological relationships that also capture pathogen population shifts.

No evidence of male-biased sexual selection in a snake with conventional Darwinian sex roles.

Decades of research on sexual selection have demonstrated that 'conventional' Darwinian sex roles are common in species with anisogamous gametes, with those species often exhibiting male-biased sexual selection. Yet, mating system characteristics such as long-term sperm storage and polyandry have the capacity to disrupt this pattern. Here, these ideas were explored by quantifying sexual selection metrics for the western diamond-backed rattlesnake (Crotalus atrox). A significant standardized sexual selection gradient was not found for males (ß SS = 0.588, p = 0.199) or females (ß SS = 0.151, p = 0.664), and opportunities for sexual selection (Is ) and selection (I) did not differ between males (Is = 0.069, I = 0.360) and females (Is = 0.284, I = 0.424; both p > 0.05). Furthermore, the sexes did not differ in the maximum intensity of precopulatory sexual selection (males: s' max = 0.155, females: s' max = 0.080; p > 0.05). Finally, there was no evidence that male snout-vent length, a trait associated with mating advantage, is a target of sexual selection (p > 0.05). These results suggest a lack of male-biased sexual selection in this population. Mating system characteristics that could erode male-biased sexual selection, despite the presence of conventional Darwinian sex roles, are discussed.

Genetic variability, community structure, and horizontal transfer of endosymbionts among three Asia II-Bemisia tabaci mitotypes in Pakistan.

Endosymbionts associated with the whitefly Bemisia tabaci cryptic species are known to contribute to host fitness and environmental adaptation. The genetic diversity and population complexity were investigated for endosymbiont communities of B. tabaci occupying different micro-environments in Pakistan. Mitotypes of B. tabaci were identified by comparative sequence analysis of the mitochondria cytochrome oxidase I (mtCOI) gene sequence. Whitefly mitotypes belonged to the Asia II-1, -5, and -7 mitotypes of the Asia II major clade. The whitefly-endosymbiont communities were characterized based on 16S ribosomal RNA operational taxonomic unit (OTU) assignments, resulting in 43 OTUs. Most of the OTUs occurred in the Asia II-1 and II-7 mitotypes (r 2 = .9, p < .005), while the Asia II-5 microbiome was less complex. The microbiome OTU groups were mitotype-specific, clustering with a basis in phylogeographical distribution and the corresponding ecological niche of their whitefly host, suggesting mitotype-microbiome co-adaptation. The primary endosymbiont Portiera was represented by a single, highly homologous OTU (0%-0.67% divergence). Two of six Arsenophonus OTUs were uniquely associated with Asia II-5 and -7, and one occurred exclusively in Asia II-1, two only in Asia II-5, and one in both Asia II-1 and -7. Four other secondary endosymbionts, Cardinium, Hemipteriphilus, Rickettsia, and Wolbachia OTUs, were found at ≤29% frequencies. The most prevalent Arsenophonus OTU was found in all three Asia II mitotypes (55% frequency), whereas the same strain of Cardinium and Wolbachia was found in both Asia II-1 and -5, and a single Hemipteriphilus OTU occurred in Asia II-1 and -7. This pattern is indicative of horizontal transfer, suggestive of a proximity between mitotypes sufficient for gene flow at overlapping mitotype ecological niches.

Biochemical and biological activities of the venom of a new species of pitviper from Vietnam, Triceratolepidophis sieversorum.

Biochemical and biological activities of a venom sample from a recently discovered new genus and species of pitviper from Vietnam, Triceratolepidophis sieversorum, were assayed and compared with those of five other viperid snakes (Bothrops asper, Crotalus atrox, Protobothrops flavoviridis, Trimeresurus insularis, and Vipera ammodytes). The venom had high casein hydrolysis, arginine ester hydrolysis and haemorrhagic activities, lacked measurable phosphodiesterase and L-amino acid oxidase activities, and had no procoagulant activity on either bovine fibrinogen or human plasma. Other enzymatic activities (phospholipase A(2), kallikrein) were moderate. The approximate i.p. LD(50) (mice) of the venom is about 5-6 mg/kg.

Mating systems, reproductive success, and sexual selection in secretive species: a case study of the western diamond-backed rattlesnake, Crotalus atrox.

Long-term studies of individual animals in nature contribute disproportionately to our understanding of the principles of ecology and evolution. Such field studies can benefit greatly from integrating the methods of molecular genetics with traditional approaches. Even though molecular genetic tools are particularly valuable for species that are difficult to observe directly, they have not been widely adopted. Here, we used molecular genetic techniques in a 10-year radio-telemetric investigation of the western diamond-backed rattlesnake (Crotalus atrox) for an analysis of its mating system and to measure sexual selection. Specifically, we used microsatellite markers to genotype 299 individuals, including neonates from litters of focal females to ascertain parentage using full-pedigree likelihood methods. We detected high levels of multiple paternity within litters, yet found little concordance between paternity and observations of courtship and mating behavior. Larger males did not father significantly more offspring, but we found evidence for size-specific male-mating strategies, with larger males guarding females for longer periods in the mating seasons. Moreover, the spatial proximity of males to mothers was significantly associated with reproductive success. Overall, our field observations alone would have been insufficient to quantitatively measure the mating system of this population of C. atrox, and we thus urge more widespread adoption of molecular tools by field researchers studying the mating systems and sexual selection of snakes and other secretive taxa.

Spatio-temporal patterns of genetic change amongst populations of cassava Bemisia tabaci whiteflies driving virus pandemics in East and Central Africa.

The greatest current threat to cassava in sub-Saharan Africa, is the continued expansion of plant virus pandemics being driven by super-abundant populations of the whitefly vector, Bemisia tabaci. To track the association of putatively genetically distinct populations of B. tabaci with pandemics of cassava mosaic disease (CMD) and cassava brown streak disease (CBSD), a comprehensive region-wide analysis examined the phylogenetic relationships and population genetics of 642 B. tabaci adults sampled from cassava in six countries of East and Central Africa, between 1997 and 2010, using a mitochondrial DNA cytochrome oxidase I marker (780 bases). Eight phylogenetically distinct groups were identified, including one, designated herein as 'East Africa 1' (EA1), not previously described. The three most frequently occurring groups comprised >95% of all samples. Among these, the Sub-Saharan Africa 2 (SSA2) group diverged by c. 8% from two SSA1 sub-groups (SSA1-SG1 and SSA1-SG2), which themselves were 1.9% divergent. During the 14-year study period, the group associated with the CMD pandemic expansion shifted from SSA2 to SSA1-SG1. Population genetics analyses of SSA1, using Tajima's D, Fu's Fs and Rojas' R2 statistics confirmed a temporal transition in SSA1 populations from neutrally evolving at the outset, to rapidly expanding from 2000 to 2003, then back to populations more at equilibrium after 2004. Based on available evidence, hybrid introgression appears to be the most parsimonious explanation for the switch from SSA2 to SSA1-SG1 in whitefly populations driving cassava virus pandemics in East and Central Africa.