Cytogenetics of the Javan file snake (Acrochordus javanicus) and the evolution of snake sex chromosomes

Advanced snakes (Caenophidia) are an important group including around 90% of the recent species of snakes. The basal splitting of the clade is still rather controversial, and it is not fully understood when the differentiation of sex chromosomes started in snake evolution. To help resolve these questions, we performed cytogenetic analysis on the Javan file snake, also known as the elephant trunk snake (Acrochordus javanicus) from the family Acrochordidae, which occupies an informative phylogenetic position. For the first time for acrochordids, we identified heteromorphic ZZ/ZW sex chromosomes with a highly heterochromatic W chromosome. These traits are likely synapomorphies of advanced snakes. In contrast to other caenophidian snakes, the Javan file snake lacks an accumulation of Bkm repeats and interstitial telomeric repeats on the W chromosome. This observation supports the sister group relationship between acrochordids and all other caenophidian snakes including the family Xenodermatidae and questions the suggested role of Bkm repeats in the formation of sex heterochromatin in snakes. The revealed partial gene content of the Z chromosome in acrochordids supports the hypothesis that the progressive degeneration of the W chromosome commenced in snakes before the basal split of Caenophidia, albeit its evolutionary rate in file snakes might be slower than in their sister lineage.     

Evolutionary stability of sex chromosomes in snakes

Amniote vertebrates possess various mechanisms of sex determination, but their variability is not equally distributed. The large evolutionary stability of sex chromosomes in viviparous mammals and birds was believed to be connected with their endothermy. However, some ectotherm lineages seem to be comparably conserved in sex determination, but previously there was a lack of molecular evidence to confirm this. Here, we document a stability of sex chromosomes in advanced snakes based on the testing of Z-specificity of genes using quantitative PCR (qPCR) across 37 snake species (our qPCR technique is suitable for molecular sexing in potentially all advanced snakes). We discovered that at least part of sex chromosomes is homologous across all families of caenophidian snakes (Acrochordidae, Xenodermatidae, Pareatidae, Viperidae, Homalopsidae, Colubridae, Elapidae and Lamprophiidae). The emergence of differentiated sex chromosomes can be dated back to about 60 Ma and preceded the extensive diversification of advanced snakes, the group with more than 3000 species. The Z-specific genes of caenophidian snakes are (pseudo)autosomal in the members of the snake families Pythonidae, Xenopeltidae, Boidae, Erycidae and Sanziniidae, as well as in outgroups with differentiated sex chromosomes such as monitor lizards, iguanas and chameleons. Along with iguanas, advanced snakes are therefore another example of ectothermic amniotes with a long-term stability of sex chromosomes comparable with endotherms.     

A workflow of massive identification and application of intron markers using snakes as a model

Relative to the commonly used mitochondrial and nuclear protein‐coding genes, the noncoding intron sequences are a promising source of informative markers that have the potential to resolve difficult phylogenetic nodes such as rapid radiations and recent divergences. Yet many issues exist in the use of intron markers, which prevent their extensive application as conventional markers. We used the diverse group of snakes as an example to try paving the way for massive identification and application of intron markers. We performed a series of bioinformatics screenings which identified appropriate introns between single‐copy and conserved exons from two snake genomes, adding particular constraints on sequence length variability and sequence variability. A total of 1,273 candidate intron loci were retrieved. Primers for nested polymerase chain reaction (PCR) were designed for over a hundred candidates and tested in 16 snake representatives. 96 intron markers were developed that could be amplified across a broad range of snake taxa with high PCR successful rates. The markers were then applied to 49 snake samples. The large number of amplicons was subjected to next‐generation sequencing (NGS). An analytic strategy was developed to accurately recover the amplicon sequences, and approximately, 76% of the marker sequences were recovered. The average p‐distances of the intron markers at interfamily, intergenus, interspecies, and intraspecies levels were .168, .052, .015, and .004, respectively, suggesting that they were useful to study snake relationships of different evolutionary depths. A snake phylogeny was constructed with the intron markers, which produced concordant results with robust support at both interfamily and intragenus levels. The intron markers provide a convenient way to explore the signals in the noncoding regions to address the controversies on the snake tree. Our improved strategy of genome screening is effective and can be applied to other animal groups. NGS coupled with appropriate sequence processing can greatly facilitate the extensive application of molecular markers.     

Distribution of genetic variation underlying adult migration timing in steelhead of the Columbia River basin

Fish migrations are energetically costly, especially when moving between freshwater and saltwater, but are a viable strategy for Pacific salmon and trout (Oncorhynchus spp.) due to the advantageous resources available at various life stages. Anadromous steelhead (O. mykiss) migrate vast distances and exhibit variation for adult migration phenotypes that have a genetic basis at candidate genes known as greb1L and rock1. We examined the distribution of genetic variation at 13 candidate markers spanning greb1L, intergenic, and rock1 regions versus 226 neutral markers for 113 populations (n = 9,471) of steelhead from inland and coastal lineages in the Columbia River. Patterns of population structure with neutral markers reflected genetic similarity by geographic region as demonstrated in previous studies, but candidate markers clustered populations by genetic variation associated with adult migration timing. Mature alleles for late migration had the highest frequency overall in steelhead populations throughout the Columbia River, with only 9 of 113 populations that had a higher frequency of premature alleles for early migration. While a single haplotype block was evident for the coastal lineage, we identified multiple haplotype blocks for the inland lineage. The inland lineage had one haplotype block that corresponded to candidate markers within the greb1L gene and immediately upstream in the intergenic region, and the second block only contained candidate markers from the intergenic region. Haplotype frequencies had similar patterns of geographic distribution as single markers, but there were distinct differences in frequency between the two haplotype blocks for the inland lineage. This may represent multiple recombination events that differed between lineages where phenotypic differences exist between freshwater entry versus arrival timing as indicated by Micheletti et al. (2018a). Redundancy analyses were used to model environmental effects on allelic frequencies of candidate markers, and significant variables were migration distance, temperature, isothermality, and annual precipitation. This study improves our understanding of the spatial distribution of genetic variation underlying adult migration timing in steelhead as well as associated environmental factors and has direct conservation and management implications.     

Multiple molecular approaches yield no evidence for sex‐determining genes in lake sturgeon (Acipenser fulvescens)

Common DNA‐based sexing assays have been widely used for the conservation and management of mammals and birds. However, many fishes do not have genetic sex determination and in those that do, the plasticity of the genes involved means that species‐specific assays are normally required. Such DNA‐sexing markers would be especially valuable in lake sturgeon (Acipenser fulvescens) because of their sexual monomorphism, delayed sexual maturity, and conservation status. We tried to identify genetic differences between male and female lake sturgeon using several different molecular genetic methods, including randomly amplified polymorphic DNA, representational difference analyses, subtractive hybridization, and a candidate gene approach. Ultimately, a number of genes were identified but none was sex‐specific. Although the ultimate mechanism of sex determination is yet unknown, it is possible that sex determination is environmental in lake sturgeon, especially since recent studies have also failed to identify sex determination genes in other sturgeon species.     

ZW,XY, and yet ZW: Sex chromosome evolution in snakes even more complicated

Snakes are historically important in the formulation of several central concepts on the evolution of sex chromosomes. For over 50 years, it was believed that all snakes shared the same ZZ/ZW sex chromosomes, which are homomorphic and poorly differentiated in “basal” snakes such as pythons and boas, while heteromorphic and well differentiated in “advanced” (caenophidian) snakes. Recent molecular studies revealed that differentiated sex chromosomes are indeed shared among all families of caenophidian snakes, but that boas and pythons evolved likely independently male heterogamety (XX/XY sex chromosomes). The historical report of heteromorphic ZZ/ZW sex chromosomes in a boid snake was previously regarded as ambiguous. In the current study, we document heteromorphic ZZ/ZW sex chromosomes in a boid snake. A comparative approach suggests that these heteromorphic sex chromosomes evolved very recently and that they are poorly differentiated at the sequence level. Interestingly, two snake lineages with confirmed male heterogamety possess homomorphic sex chromosomes, but heteromorphic sex chromosomes are present in both snake lineages with female heterogamety. We point out that this phenomenon is more common across squamates. The presence of female heterogamety in non‐caenophidian snakes indicates that the evolution of sex chromosomes in this lineage is much more complex than previously thought, making snakes an even better model system for the evolution of sex chromosomes.     

An Analysis of Snake-Scent Application Behaviour in Siberian Chipmunks (Eutamias sibiricus asiaticus)

Some rodents assault live snakes by mobbing. Another snake-directed behaviour is shown by Siberian chipmunks, Eutamias sibiricus asiaticus, when encountering a dead snake. They approach the carcass cautiously, gnaw the snake skin, and chew and apply the gnawed bits to their body fur. We have termed this behaviour “Snake-Scent Application” (SSA). SSA behaviour is also elicited by snake urine and feces. Chemicals in snake skin and rectal and cloacal sacs release SSA. Snake urine applied to dead mice tends to suppress ingestion by snakes.     

Ticks on snakes: The ecological correlates of ectoparasite infection in free‐ranging snakes in tropical Australia

Parasites profoundly influence the lives of their hosts, yet the dynamics of host–parasite interactions are poorly understood – especially in reptiles. We examined the ecological correlates of parasitism by ixodid ticks in an assemblage of 10 snake species in tropical Australia. In total, we recorded 3803 ticks on 1841 individual snakes of six species (no ticks were found on the other species). Molecular analyses confirmed the tropical reptile tick (Amblyomma fimbriatum: Ixodidae) to be the most common snake tick at our study site, with inter‐ and intraspecific variation in tick prevalence and intensity. Tick attachment sites were random on most snake species, but both male and female ticks congregated on the heads of the colubrid snake Boiga irregularis and the python Simalia amethistina. In these same species, tick loads were higher on snakes captured in woodland than in rainforest. Females of two python species (Aspidites melanocephalus and S. amethistina) had higher tick loads than did males. In B. irregularis, individuals captured in the dry season had higher tick loads than those captured in the wet season. In most parasitized snake species, larger individuals had greater tick loads. Data from snake recaptures confirmed individual tick burdens frequently varied, with little correlation between tick loads on the same snake at successive captures (except for B. irregularis). Finally, tick intensity was not correlated with (and thus, presumably did not influence) the body condition of any snake species in our study. Use of specific types of refuge sites may strongly influence tick loads on snakes in this system.     

Recent rapid speciation and ecomorph divergence in Indo‐Australian sea snakes

The viviparous sea snakes (Hydrophiinae) are a young radiation of at least 62 species that display spectacular morphological diversity and high levels of local sympatry. To shed light on the mechanisms underlying sea snake diversification, we investigated recent speciation and eco‐morphological differentiation in a clade of four nominal species with overlapping ranges in Southeast Asia and Australia. Analyses of morphology and stomach contents identified the presence of two distinct ecomorphs: a ‘macrocephalic’ ecomorph that reaches >2 m in length, has a large head and feeds on crevice‐dwelling eels and gobies; and a ‘microcephalic’ ecomorph that rarely exceeds 1 m in length, has a small head and narrow fore‐body and hunts snake eels in burrows. Mitochondrial sequences show a lack of reciprocal monophyly between ecomorphs and among putative species. However, individual assignment based on newly developed microsatellites separated co‐distributed specimens into four significantly differentiated clusters corresponding to morphological species designations, indicating limited recent gene flow and progress towards speciation. A coalescent species tree (based on mitochondrial and nuclear sequences) and isolation‐migration model (mitochondrial and microsatellite markers) suggest between one and three transitions between ecomorphs within the last approximately 1.2 million to approximately 840 000 years. In particular, the macrocephalic ‘eastern’ population of Hydrophis cyanocinctus and microcephalic H. melanocephalus appear to have diverged very recently and rapidly, resulting in major phenotypic differences and restriction of gene flow in sympatry. These results highlight the viviparous sea snakes as a promising system for speciation studies in the marine environment.     

Targeted resistance gene mapping in soybean using modified AFLPs

The soybean [Glycine max (Merr.) L.] linkage group F contains a vital region of clustered genes for resistance to numerous pathogens including the soybean mosaic virus resistance gene, Rsv1. In order to develop new genetic markers that map to this gene cluster, we employed a targeted approach that utilizes the speed and high-throughput of AFLP, but modified it to incorporate sequence information from the highly conserved nucleotide binding site (NBS) region of cloned disease resistance genes. By using a labeled degenerate primer corresponding to the p-loop portion of the NBS region of resistance genes, such as N, L6, and Rps2, we were able to quickly amplify numerous polymorphic bands between parents of a population segregating for resistance to Rsv1. Of these polymorphic bands, bulk segregant analysis revealed four markers that were closely linked to Rsv1. These markers were cloned and used as probes for RFLP analysis. The four clones mapped to within a 6-cM region surrounding Rsv1, the closest being 0.4 cM away from the gene. Sequence analysis showed that all four clones contain the p-loop sequence corresponding to the degenerate primer and that one of the four clones contains an open reading frame sequence which when translated is related to the NBS region of other cloned disease resistance genes. The rapid identification of four markers closely linked to Rsv1 in soybean demonstrates the utility of this method for generating markers tightly linked to important plant disease resistance genes. Received: 25 September 1999 / Accepted: 3 November 1999     

Expression of tomato reference genes using established primer sets: Stability across experimental set‐ups

Reliable reference genes are critical for relative quantification using quantitative real‐time PCR (qPCR). Ten tomato genes (Solanum lycopersicum) and their respective primer sets, which have been used over the last 6 years as references in expression studies, were evaluated for their performance using leaf tissue samples grown under semi‐controlled conditions and infected with grey mould (Botrytis cinerea) or late blight (Phytophthora infestans). The target genes coding for U6 snRNA‐associated Sm‐like protein LSm7, calcineurin B‐like protein and V‐type proton ATPase were the most stable expressed of all the genes tested in three experimental repetitions. Evaluation of candidate reference genes with geNorm and NormFinder softwares yielded the lowest mean values for their respective primer sets LSM7, SlCBL1 and SlATPase, suggesting stable expression. However, SlATPase primer set revealed a comparably high intra‐group variation and was thus not considered further. In follow‐up experiments with P. infestans, the geNorm and NormFinder values of primer sets LSM7 and SlCBL1 were even lower, indicating the stability of their expression also under these conditions. Primer efficiency differed by ‐18 to +5 percentage points from values presented in the literature. Our findings show that a reference primer set which delivers the best results in one system may be outperformed by another under different experimental conditions, thus recommending a reassessment of both expression stability and qPCR efficiency whenever the biological or technical experimental set‐up is changed. On the basis of our results, we recommend the use of LSM7 and SlCBL1 as reference primer sets for gene expression studies on plant tissue derived from open or semi‐controlled conditions.     

Reappraising plastid markers of the red algae for phylogenetic community ecology in the genomic era

Selection of appropriate genetic markers to quantify phylogenetic diversity is crucial for community ecology studies. Yet, systematic evaluation of marker genes for this purpose is scarcely done. Recently, the combined effort of phycologists has produced a rich plastid genome resource with taxonomic representation spanning all of the major lineages of the red algae (Rhodophyta). In this proof‐of‐concept study, we leveraged this resource by developing and applying a phylogenomic strategy to seek candidate plastid markers suitable for phylogenetic community analysis. We ranked the core genes of 107 published plastid genomes based on various sequence‐derived properties and their tree distance to plastid genome phylogenies. The resulting ranking revealed that the most widely used marker, rbcL, is not necessarily the optimal marker, while other promising markers might have been overlooked. We designed and tested PCR primers for several candidate marker genes, and successfully amplified one of them, rpoC1, in a taxonomically broad set of red algal specimens. We suggest that our general marker identification methodology and the rpoC1 primers will be useful to the phycological community for investigating the biodiversity and community ecology of the red algae.     

Sex Identification of Four Penguin Species Using Locus‐Specific PCR

Traditional methods for sex identification are not applicable to sexually monomorphic species, leading to difficulties in the management of their breeding programs. To identify sex in sexually monomorphic birds, molecular methods have been established. Two established primer pairs (2550F/2718R and p8/p2) amplify the CHD1 gene region from both the Z and W chromosomes. Here, we evaluated the use of these primers for sex identification in four sexually monomorphic penguin species: king penguins (Aptenodytes patagonicus), rockhopper penguins (Eudyptes chrysocome), gentoo penguins (Pygoscelis papua), and Magellanic penguins (Spheniscus magellanicus). For all species except rockhopper penguins, primer pair 2550F/2718R resulted in two distinct CHD1Z and CHD1W PCR bands, allowing for sex identification. For rockhopper penguins, only primer pair p8/p2 yielded different CHD1Z and CHD1W bands, which were faint and similar in size making them difficult to distinguish. As a result, we designed a new primer pair (PL/PR) that efficiently determined the gender of individuals from all four penguin species. Sequencing of the PCR products confirmed that they were from the CHD1 gene region. Primer pair PL/PR can be evaluated for use in sexing other penguin species, which will be crucial for the management of new penguin breeding programs. Zoo Biol 32:257–261, 2013. © 2012 Wiley Periodicals, Inc.     

Candidate gene polymorphisms for behavioural adaptations during urbanization in blackbirds

Successful urban colonization by formerly rural species represents an ideal situation in which to study adaptation to novel environments. We address this issue using candidate genes for behavioural traits that are expected to play a role in such colonization events. We identified and genotyped 16 polymorphisms in candidate genes for circadian rhythms, harm avoidance and migratory and exploratory behaviour in 12 paired urban and rural populations of the blackbird Turdus merula across the Western Palaearctic. An exonic microsatellite in the SERT gene, a candidate gene for harm avoidance behaviour, exhibited a highly significant association with habitat type in an analysis conducted across all populations. Genetic divergence at this locus was consistent in 10 of the 12 population pairs; this contrasts with previously reported stochastic genetic divergence between these populations at random markers. Our results indicate that behavioural traits related to harm avoidance and associated with the SERT polymorphism experience selection pressures during most blackbird urbanization events. These events thus appear to be influenced by homogeneous adaptive processes in addition to previously reported demographic founder events.     

CRISPR/Cas9‐mediated mutagenesis of the white gene in the tephritid pest Bactrocera tryoni

The Queensland fruit fly, Bactrocera tryoni (Froggatt), is a polyphagous horticultural pest in Australia that is capable of causing significant damage to more than 100 different host fruits and vegetables. Chemical applications and ecological control strategies, such as the sterile insect technique (SIT), are commonly used to suppress established populations and eradicate invasive outbreaks following migration. The recently published B. tryoni draft genome provides new opportunities to identify candidate genes for targeted genome modification in order to generate advanced genetic strains for management using sterile insect strategies. Here, we demonstrate CRISPR/Cas‐mediated mutagenesis in B. tryoni through generating a series of frame‐shift mutations in the ATP‐dependent binding cassette transporter, white, causing a classic white‐eye phenotype. This work establishes methods for CRISPR/Cas genome editing in tephritids and demonstrates its potential for developing genetic sexing strains which could be used for SIT‐based pest control.     

Multiplex PCR assay for rapid identification of three endangered snake species of India

Species identification has been the core issue in all approaches of conservation of endangered wild life. In this regard molecular techniques for species authentication have proved indispensable. A novel multiplex PCR assay for the identification of three Indian snake species Python morulus, Ptyas mucosus, and Naja naja is successfully demonstrated using 16S rRNA gene. Three reverse primers and a common forward primer were designed to generate three different size species-specific PCR fragments. Absence of any PCR amplification in non-target species proves the specificity of the primers. These four primers were combined in a multiplex assay to enable identification of three snake species in a single reaction. The assay described here shows its utility in identifying unknown snake specimen and in case of samples yielding low quality DNA. This multiplex PCR technique using novel primers is an unprecedented approach offered for forensic identification of exhibits originating from three Indian snake species. It is expected that this endeavor will help strengthening conservation efforts for these species.     

A multiplex PCR assay for molecular sexing of the naked mole-rat (Heterocephalus glaber)

For molecular sexing of the naked mole-rat (Heterocephalus glaber), we designed a PCR primer set to amplify part of the Y-linked DBY gene. When this primer set was applied to the samples of known sex with the 16S rRNA gene (16S rDNA) primers as control, PCR products were successfully obtained as two DNA bands in males, a male-specific 163 bp DBY band and a 446 bp band of 16S rDNA shared with females, whereas females showed only the common band. This result shows that this multiplex PCR assay is useful for sex identification of H. glaber.     

Scale‐dependent Orientation in Movement Paths: A Case Study of an African Viper

Decisions relating to the orientation of movement by animals and how this translates into movement patterns can occur at multiple spatial scales simultaneously, but this interaction is poorly understood for many groups of animals. Using the tracks left by moving snakes in their sandy habitat, we studied the movement paths of the African snake Bitis schneideri (Namaqua dwarf adder) for evidence of broad‐scale directional persistence and short‐range avoidance of exposure. Although snakes clearly displayed directional persistence, they preferentially moved to nearby shrubs, thereby minimizing exposure to solar and thermal radiation and/or predation. Thus, snakes made decisions relating to orientation at a minimum of two scales, the interaction of which resulted in snakes moving ≈17% (mean straightness index = 0.85) further than the simple broad‐scale straight‐line distance. We assert that the actual path chosen by moving snakes represents a trade‐off of various costs and risks that include risk of predation, exposure to the elements, time and energy expenditure. Our study highlights the need for cognizance of the possibility of the scale dependence of orientation and movement in studies of snake movement, and adds to a growing literature demonstrating previously unrecognized behavioural complexity in snakes.     

Molecular sexing of birds: A comparative review of polymerase chain reaction (PCR)-based methods

Accurate identification of sex in birds is important for the management and conservation of avian wildlife in several ways, namely in the development of population, behavioral and ecological studies, as well as in the improvement of ex situ captive breeding programs. In general, nestlings, juveniles and adult birds of a wide number of sexually monomorphic species cannot be sexed based on phenotypic traits. The development of molecular methodologies for avian sexing overcame these difficulties, allowing a reliable gender differentiation for these species. The polymerase chain reaction (PCR)-based methods have been widely applied in molecular sexing of birds, using a large diversity of sex-linked markers. During the last 15 yrs, there was a continuous improvement in the PCR-based protocols for bird sexing, increasing the accuracy, speed and high-throughput applicability of these techniques. The recent advances in real-time PCR platforms and whole genome analysis methods provided new resources for the detection and analysis of novel specific markers and protocols. This review presents a comparative guide of classical and recent advances in PCR-based methods for avian molecular sexing, highlighting its strengths and limitations. Future research opportunities in this field are also addressed.