Search

In the present study, we describe an interspecific kleptoparasitic interaction between two sympatric mammalian carnivores in the high altitudinal Trans-Himalaya region of Himachal Pradesh, India. The study was based on the inferences drawn from the circumstantial evidence (direct and indirect) noticed in the study area in Pin Valley National Park. The inferences from the analysis of the evidence suggested the interaction between a Snow Leopard Panthera uncia, a Red Fox Vulpes vulpes, and a donkey. The arrangement of evidence in a sequential manner suggested that a donkey was killed by a Snow Leopard and a Red Fox stole the food from the carrion of the Snow Leopard’s prey. The Red Fox was killed by the Snow Leopard, which was caught while stealing. The present study represents an example of kleptoparasitic interaction between the Snow Leopard and the Red Fox. This study also proves that such interactions may cost the life of a kleptoparasite and supports the retaliation behaviour of Snow Leopards.

Topologically associating domains (TADs) are thought to play an important role in preventing gene misexpression by spatially constraining enhancer–promoter contacts. The deleterious nature of gene misexpression implies that TADs should, therefore, be conserved among related species. Several early studies comparing chromosome conformation between species reported high levels of TAD conservation; however, more recent studies have questioned these results. Furthermore, recent work suggests that TAD reorganization is not associated with extensive changes in gene expression. Here, we investigate the evolutionary conservation of TADs among 11 species of Drosophila. We use Hi-C data to identify TADs in each species and employ a comparative phylogenetic approach to derive empirical estimates of the rate of TAD evolution. Surprisingly, we find that TADs evolve rapidly. However, we also find that the rate of evolution depends on the chromatin state of the TAD, with TADs enriched for developmentally regulated chromatin evolving significantly slower than TADs enriched for broadly expressed, active chromatin. We also find that, after controlling for differences in chromatin state, highly conserved TADs do not exhibit higher levels of gene expression constraint. These results suggest that, in general, most TADs evolve rapidly and their divergence is not associated with widespread changes in gene expression. However, higher levels of evolutionary conservation and gene expression constraints in TADs enriched for developmentally regulated chromatin suggest that these TAD subtypes may be more important for regulating gene expression, likely due to the larger number of long-distance enhancer–promoter contacts associated with developmental genes.

The KwaZulu-Natal sardine run, popularly known as the “greatest shoal on Earth,” is a mass migration of South African sardines from their temperate core range into the subtropical Indian Ocean. It has been suggested that this represents the spawning migration of a distinct subtropical stock. Using genomic and transcriptomic data from sardines collected around the South African coast, we identified two stocks, one cool temperate (Atlantic) and the other warm temperate (Indian Ocean). Unexpectedly, we found that sardines participating in the sardine run are primarily of Atlantic origin and thus prefer colder water. These sardines separate from the warm-temperate stock and move into temporarily favorable Indian Ocean habitat during brief cold-water upwelling periods. Once the upwelling ends, they find themselves trapped in physiologically challenging subtropical habitat and subject to intense predation pressure. This makes the sardine run a rare example of a mass migration that has no apparent fitness benefits.

AimBiogeographical transition zones are areas of overlap between the faunas of adjacent biogeographical entities. Particularly, the well-defined transition zones along linear coastlines are interesting natural laboratories to study dispersal and incipient speciation. Few studies have explored whether marine biogeographical transition zones harbour biodiversity that is distinct from that of the biogeographical entities they separate. The Wild Coast in eastern South Africa is a poorly-studied transition zone between region’s warm-temperate and subtropical faunas, and is generally considered to be an area of faunal overlap.LocationThe South African portion of the Western Indian OceanMethodsSequences of the DNA barcoding marker COI were generated from 306 estuarine sandprawns (Callichirus kraussi) collected at 13 sites. Genetic structure and evolutionary history were assessed using a haplotype network and a Bayesian discrete phylogeographic analysis.ResultTwo populations were identified whose ranges are centred on the Wild Coast, a rare one in the northern portion and a more common one in the central and southern portion of this biogeographical transition zone. These populations are not closely related to each other, but descend from subtropical and warm-temperate sister populations, respectively.ConclusionsThis is the first study to indicate that the Wild Coast marine biogeographical transition zone is not merely an area of faunal overlap, and one of very few studies to have discovered genetically unique populations within a marine biogeographical transition zone. The Wild Coast may harbour additional unique biodiversity that remains to be discovered, including rare species that require protection. More research is required to understand how this environmentally dynamic marine biogeographical transition zone differs from the adjacent biogeographical provinces.

Topologically associating domains, or TADs, are functional units that organize chromosomes into 3D structures of interacting chromatin. TADs play an important role in regulating gene expression by constraining enhancer-promoter contacts; there is evidence that deletion of TAD boundaries leads to aberrant expression of neighboring genes. While the mechanisms of TAD formation have been well-studied, current knowledge on the extent of TAD conservation across species is inconclusive. Due to the integral role TADs play in gene regulation, their structure and organization is expected to be conserved during evolution. However, more recent research suggests that TAD structures diverge relatively rapidly. We use Hi-C chromosome conformation capture to measure evolutionary conservation of whole TADs and TAD boundary elements between D. melanogaster and D. triauraria, two early-branching species from the melanogaster species group which diverged ~15 million years ago. We find that 75% of TAD boundaries are orthologous while only 25% of TAD domains are conserved and these are enriched for Polycomb-repressed chromatin. Our results show that TADs have been reorganized since the common ancestor of D. melanogaster and D. triauraria, yet the sequence elements that specify TAD boundaries remain highly conserved. We propose that evolutionary divergence in 3D genome organization results from shuffling of conserved boundary elements across chromosomes, breaking old TADs and creating new TAD architectures. This result supports the existence of distinct TAD subtypes: some may be evolutionarily flexible while others remain highly conserved due to their importance in restricting gene-regulatory element interactions.

The Indian cheetah was hunted to extinction by the mid-20th century. While analysis of 139 bp of mitochondrial DNA (mtDNA) has confirmed that the Indian cheetah was part of the Asiatic subspecies (Acinonyx jubatus venaticus), the detailed relationships between cheetah populations remains unclear due to limited genetic data. We clarify these relationships by studying larger fragments of cheetah mtDNA, both from an Indian cheetah museum specimen and two African cheetah, one modern and one historic, imported into India at different times. Our results suggest that the most recent common ancestor of cheetah mtDNA is approximately twice as ancient as currently recognised. The Indian and Southeast African (Acinonyx jubatus jubatus) cheetah mtDNA diverged approximately 72 kya, while the Southeast and Northeast African (Acinonyx jubatus soemmeringii) cheetah mtDNA diverged around 139 kya. Additionally, the historic African cheetah sampled from India proved to have an A. j. jubatus haplotype, suggesting a hitherto unrecognised South African route of cheetah importation into India in the 19th century. Together, our results provide a deeper understanding of the relationships between cheetah subspecies, and have important implications for the conservation of A. j. venaticus and potential reintroduction of cheetahs into India.

Intraspecific genetic structure in widely distributed marine species often mirrors the boundaries between temperature-defined bioregions. This suggests that the same thermal gradients that maintain distinct species assemblages also drive the evolution of new biodiversity. Ecological speciation scenarios are often invoked to explain such patterns, but the fact that adaptation is usually only identified when phylogenetic splits are already evident makes it impossible to rule out the alternative scenario of allopatric speciation with subsequent adaptation. We integrated large-scale genomic and environmental datasets along one of the world's best-defined marine thermal gradients (the South African coastline) to test the hypothesis that incipient ecological speciation is a result of divergence linked to the thermal environment. We identified temperature-associated gene regions in a coastal fish species that is spatially homogeneous throughout several temperature-defined biogeographic regions based on selectively neutral markers. Based on these gene regions, the species is divided into geographically distinct regional populations. Importantly, the ranges of these populations are delimited by the same ecological boundaries that define distinct infraspecific genetic lineages in co-distributed marine species, and biogeographic disjunctions in species assemblages. Our results indicate that temperature-mediated selection represents an early stage of marine ecological speciation in coastal regions that lack physical dispersal barriers.

Supplementary methods, figures and tables

Topologically associating domains, or TADs, are functional units that organize chromosomes into 3D structures of interacting chromatin. TADs play an important role in regulating gene expression by constraining enhancer-promoter contacts and there is evidence that deletion of TAD boundaries leads to aberrant expression of neighboring genes. While the mechanisms of TAD formation have been well-studied, current knowledge on the patterns of TAD evolution across species is limited. Due to the integral role TADs play in gene regulation, their structure and organization is expected to be conserved during evolution. However, more recent research suggests that TAD structures diverge relatively rapidly. We use Hi-C chromosome conformation capture to measure evolutionary conservation of whole TADs and TAD boundary elements between D. melanogaster and D. triauraria, two early-branching species from the melanogaster species group which diverged ∼15 million years ago. We find that the majority of TADs have been reorganized since the common ancestor of D. melanogaster and D. triauraria, via a combination of chromosomal rearrangements and gain/loss of TAD boundaries. TAD reorganization between these two species is associated with a localized effect on gene expression, near the site of disruption. By separating TADs into subtypes based on their chromatin state, we find that different subtypes are evolving under different evolutionary forces. TADs enriched for broadly expressed, transcriptionally active genes are evolving rapidly, potentially due to positive selection, whereas TADs enriched for developmentally-regulated genes remain conserved, presumably due to their importance in restricting gene-regulatory element interactions. These results provide novel insight into the evolutionary dynamics of TADs and help to reconcile contradictory reports related to the evolutionary conservation of TADs and whether changes in TAD structure affect gene expression., Author summary Chromosomes are organized inside the three-dimensional nucleus in a hierarchical fashion, from the wrapping of DNA around nucleosomes to the presence of compartments of active and inactive chromatin. Here, we focus on domains of interacting chromatin known as topologically associating domains, or TADs. TADs are dozens to hundreds of kilobases in size and tend to contain genes with similar epigenetic states. While much research has been done on the mechanisms of TAD formation, less is known about the role of TADs in genome evolution. The research that does exist is contradictory, with some studies claiming that TADs are highly conserved and others reporting that TADs diverge rapidly. We compared the TAD locations between two Drosophila species and found that the majority have been reorganized. We also identified TAD subtypes that differ in terms of their evolutionary conservation: TADs that contain developmental genes are less likely to be reorganized while TADs that contain broadly expressed genes evolve rapidly. We found that TAD reorganization has an effect on gene expression, but only for genes located near the site of disruption. These results provide new insight into the evolution of 3D genome organization and help to reconcile previous contradictions in the literature.

Tests for isolation by distance (IBD) are the most commonly used method of assessing spatial genetic structure. Many studies have exclusively used mitochondrial DNA (mtDNA) sequences to test for IBD, but this marker is often in conflict with multilocus markers. Here, we report a review of the literature on IBD, with the aims of determining (a) whether significant IBD is primarily a result of lumping spatially discrete populations, and (b) whether microsatellite datasets are more likely to detect IBD when mtDNA does not. We also provide empirical data from four species in which mtDNA failed to detect IBD by comparing these with microsatellite and SNP data. Our results confirm that IBD is mostly found when distinct regional populations are pooled, and this trend disappears when each is analysed separately. Discrepancies between markers were found in almost half of the studies reviewed, and microsatellites were more likely to detect IBD when mtDNA did not. Our empirical data rejected the lack of IBD in the four species studied, and support for IBD was particularly strong for the SNP data. We conclude that mtDNA sequence data are often not suitable to test for IBD, and can be misleading about species’ true dispersal potential. The observed failure of mtDNA to reliably detect IBD, in addition to being a single-locus marker, is likely a result of a selection-driven reduction in genetic diversity obscuring spatial genetic differentiation.

Southern Africa is a marine biodiversity hotspot that not only comprises faunal elements from the Atlantic and Indian Oceans, but also large numbers of endemic species. Using mitochondrial and nuclear DNA sequence data, we explored whether genetic structure in the endemic coastal goby Psammogobius knysnaensis, a species whose range straddles both biomes, is linked to the boundary between the two oceans. Subtle genetic structure was identified between Atlantic and Indian Ocean populations, with genetic diversity being lower in the Atlantic, and particularly on the west coast. Our results point to partial isolation between the populations associated with each biome that is most likely driven by the region’s oceanography, but unlike in other species, there is no evidence for distinct regional evolutionary lineages that are likely adapted to the environmental conditions prevalent in each region. The exclusive presence of P. knysnaensis in sheltered habitats (estuaries and lagoons) may protect this species from the severe impacts of cold water upwelling on the west coast.

The complete mitochondrial genome of sequence 16,859 bp of Indian clouded leopard (Neofelis nebulosa) has been sequenced using next generation sequencing technology Torrent PGM platform. The complete mitochondrial genome sequence of clouded leopard consists of 13 protein-coding, 22 tRNA, and two rRNA genes and a control region (CR). The mitochondrial genome is relatively similar to other felid mitochondrial genomes with respect to gene arrangement, composition, tRNA structures and skews of AT/GC bases to be typical of those reported for other mammals. The nucleotide composition of the genome shows that there is more A–T% than G–C% on the positive strand as revealed by positive AT and CG skews. The base composition of the mitochondrial genome of clouded leopard is as follows: A, 5362 bp (31.8%); C, 4560bp (27.0%); G, 2475 bp (14.6%); T, 4462 bp (26.4%). Most of the genes have ATG initiation codons, except ND1, ND2, ND3, ND4, ND6, and CYTB (ATA start codon).

Supplementary methods, figures and tables