Mode and Tempo of 3D Genome Evolution in Drosophila.

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.

3D genome evolution and reorganization in the Drosophila melanogaster species group.

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.

Ancient mtDNA from the extinct Indian cheetah supports unexpectedly deep divergence from African cheetahs.

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.

Thermal selection as a driver of marine ecological speciation.

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.

Mitochondrial DNA is unsuitable to test for isolation by distance.

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.

Complete mitochondrial genome sequence of the Indian clouded leopard (Neofelis nebulosa).

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).