A centromere satellite concomitant with extensive karyotypic diversity across the Peromyscus genus defies predictions of molecular drive.

A common feature of eukaryotic centromeres is the presence of large tracts of tandemly arranged repeats, known as satellite DNA. However, these centromeric repeats appear to experience rapid evolution under forces such as molecular drive and centromere drive, seemingly without consequence to the integrity of the centromere. Moreover, blocks of heterochromatin within the karyotype, including the centromere, are hotspots for chromosome rearrangements that may drive speciation events by contributing to reproductive isolation. However, the relationship between the evolution of heterochromatic sequences and the karyotypic dynamics of these regions remains largely unknown. Here, we show that a single conserved satellite DNA sequence in the order Rodentia of the genus Peromyscus localizes to recurrent sites of chromosome rearrangements and heterochromatic amplifications. Peromyscine species display several unique features of chromosome evolution compared to other Rodentia, including stable maintenance of a strict chromosome number of 48 among all known species in the absence of any detectable interchromosomal rearrangements. Rather, the diverse karyotypes of Peromyscine species are due to intrachromosomal variation in blocks of repeated DNA content. Despite wide variation in the copy number and location of repeat blocks among different species, we find that a single satellite monomer maintains a conserved sequence and homogenized tandem repeat structure, defying predictions of molecular drive. The conservation of this satellite monomer results in common, abundant, and large blocks of chromatin that are homologous among chromosomes within one species and among diverged species. Thus, such a conserved repeat may have facilitated the retention of polymorphic chromosome variants within individuals and intrachromosomal rearrangements between species-both factors that have previously been hypothesized to contribute towards the extremely wide range of ecological adaptations that this genus exhibits.

Adaptation and conservation insights from the koala genome.

The koala, the only extant species of the marsupial family Phascolarctidae, is classified as 'vulnerable' due to habitat loss and widespread disease. We sequenced the koala genome, producing a complete and contiguous marsupial reference genome, including centromeres. We reveal that the koala's ability to detoxify eucalypt foliage may be due to expansions within a cytochrome P450 gene family, and its ability to smell, taste and moderate ingestion of plant secondary metabolites may be due to expansions in the vomeronasal and taste receptors. We characterized novel lactation proteins that protect young in the pouch and annotated immune genes important for response to chlamydial disease. Historical demography showed a substantial population crash coincident with the decline of Australian megafauna, while contemporary populations had biogeographic boundaries and increased inbreeding in populations affected by historic translocations. We identified genetically diverse populations that require habitat corridors and instituting of translocation programs to aid the koala's survival in the wild.

Genome of the Tasmanian tiger provides insights into the evolution and demography of an extinct marsupial carnivore.

The Tasmanian tiger or thylacine (Thylacinus cynocephalus) was the largest carnivorous Australian marsupial to survive into the modern era. Despite last sharing a common ancestor with the eutherian canids ~160 million years ago, their phenotypic resemblance is considered the most striking example of convergent evolution in mammals. The last known thylacine died in captivity in 1936 and many aspects of the evolutionary history of this unique marsupial apex predator remain unknown. Here we have sequenced the genome of a preserved thylacine pouch young specimen to clarify the phylogenetic position of the thylacine within the carnivorous marsupials, reconstruct its historical demography and examine the genetic basis of its convergence with canids. Retroposon insertion patterns placed the thylacine as the basal lineage in Dasyuromorphia and suggest incomplete lineage sorting in early dasyuromorphs. Demographic analysis indicated a long-term decline in genetic diversity starting well before the arrival of humans in Australia. In spite of their extraordinary phenotypic convergence, comparative genomic analyses demonstrated that amino acid homoplasies between the thylacine and canids are largely consistent with neutral evolution. Furthermore, the genes and pathways targeted by positive selection differ markedly between these species. Together, these findings support models of adaptive convergence driven primarily by cis-regulatory evolution.

DAX1/NR0B1 was expressed during mammalian gonadal development and gametogenesis before it was recruited to the eutherian X chromosome.

The nuclear receptor subfamily 0, group B, member 1 (NR0B1) gene is an orphan nuclear receptor that is X-linked in eutherian mammals and plays a critical role in the establishment and function of the hypothalamic-pituitary-adrenal-gonadal axis. Duplication or overexpression of NR0B1 in eutherian males causes male to female sex reversal, and mutation and deletions of NR0B1 cause testicular defects. Thus, gene dosage is critical for the function of NR0B1 in normal gonadogenesis. However, NR0B1 is autosomal in all noneutherian vertebrates, including marsupials and monotreme mammals, and two active copies of the gene are compatible with both male and female gonadal development. In the current study, we examined the evolution and expression of autosomal NR0B1 during gonadal development in a marsupial (the tammar wallaby) as compared to the role of its X-linked orthologues in a eutherian (the mouse). We show that NR0B1 underwent rapid evolutionary change when it relocated from its autosomal position in the nonmammalian vertebrates, monotremes, and marsupials to an X-linked location in eutherian mammals. Despite the acquisition of a novel genomic location and a unique N-terminal domain, NR0B1 protein distribution was remarkably similar between mice and marsupials both throughout gonadal development and during gamete formation. A conserved accumulation of NR0B1 protein was observed in developing oocytes, where its function appears to be critical in the early embryo, prior to zygotic genome activation. Together these findings suggest that NR0B1 had a conserved role in gonadogenesis that existed long before it moved to the X chromosome and despite undergoing significant evolutionary change.

Redefining the human oral mycobiome with improved practices in amplicon-based taxonomy: discovery of Malassezia as a prominent commensal.

Fungi are a large, complex group, increasingly recognized as emerging threats. Their roles as modifiers of health mandate accurate portrayals of fungal communities in humans. As an entry point into the airways and gastrointestinal tract, fungi in the mouth are relevant to several biocompartments. We have revised current practices in sequence-based taxonomy assignments and employed the improvements to address the question of the fungal genera present in the healthy human mouth. The human oral mycobiome was surveyed using massively parallel, high throughput sequencing of internal transcribed spacer 1 (ITS1) amplicons from saliva following robust extraction methods. Taxonomy was assigned by comparison to a curated reference dataset, followed by filtering with an empirically determined BLAST E-value match statistic (10(-42)). Nomenclature corrections further refined results by conjoining redundant names for a single fungal genus. Following these curation steps, about two-thirds of the initially identified genera were eliminated. In comparison with the one similar metagenomic study and several earlier culture-based ones, our findings change the current conception of the oral mycobiome, especially with the discovery of the high prevalence and abundance of the genus Malassezia. Previously identified as an important pathogen of the skin, and recently reported as the predominant fungal genus at the nostril and backs of the head and ear, this is the first account of Malassezia in the human mouth. Findings from this study were in good agreement with others on the existence of many consensus members of the core mycobiome, and on unique patterns for individual subjects. This research offered a cautionary note about unconditional acceptance of lengthy lists of community members produced by automated assignments, provided a roadmap for enhancing the likely biological relevance of sequence-based fungal surveys, and built the foundation for understanding the role of fungi in health and disease of the oral cavity.

Enhancing genome assemblies by integrating non-sequence based data.

INTRODUCTION: Many genome projects were underway before the advent of high-throughput sequencing and have thus been supported by a wealth of genome information from other technologies. Such information frequently takes the form of linkage and physical maps, both of which can provide a substantial amount of data useful in de novo sequencing projects. Furthermore, the recent abundance of genome resources enables the use of conserved synteny maps identified in related species to further enhance genome assemblies. METHODS: The tammar wallaby (Macropus eugenii) is a model marsupial mammal with a low coverage genome. However, we have access to extensive comparative maps containing over 14,000 markers constructed through the physical mapping of conserved loci, chromosome painting and comprehensive linkage maps. Using a custom Bioperl pipeline, information from the maps was aligned to assembled tammar wallaby contigs using BLAT. This data was used to construct pseudo paired-end libraries with intervals ranging from 5-10 MB. We then used Bambus (a program designed to scaffold eukaryotic genomes by ordering and orienting contigs through the use of paired-end data) to scaffold our libraries. To determine how map data compares to sequence based approaches to enhance assemblies, we repeated the experiment using a 0.5× coverage of unique reads from 4 KB and 8 KB Illumina paired-end libraries. Finally, we combined both the sequence and non-sequence-based data to determine how a combined approach could further enhance the quality of the low coverage de novo reconstruction of the tammar wallaby genome. RESULTS: Using the map data alone, we were able order 2.2% of the initial contigs into scaffolds, and increase the N50 scaffold size to 39 KB (36 KB in the original assembly). Using only the 0.5× paired-end sequence based data, 53% of the initial contigs were assigned to scaffolds. Combining both data sets resulted in a further 2% increase in the number of initial contigs integrated into a scaffold (55% total) but a 35% increase in N50 scaffold size over the use of sequence-based data alone. CONCLUSIONS: We provide a relatively simple pipeline utilizing existing bioinformatics tools to integrate map data into a genome assembly which is available at http://www.mcb.uconn.edu/fac.php?name=paska. While the map data only contributed minimally to assigning the initial contigs to scaffolds in the new assembly, it greatly increased the N50 size. This process added structure to our low coverage assembly, greatly increasing its utility in further analyses.

Phylogenetic diversification patterns and divergence times in ground beetles (Coleoptera: Carabidae: Harpalinae).

BACKGROUND: Harpalinae is a species rich clade of carabid beetles with many unusual morphological forms and ecological interactions. How this diversity evolved has been difficult to reconstruct, perhaps because harpalines underwent a rapid burst of diversification early in their evolutionary history. Here we investigate the tempo of evolution in harpalines using molecular divergence dating techniques and explore the rates of lineage accumulation in harpalines and their sister group. RESULTS: According to molecular divergence date estimates, harpalines originated in the mid Cretaceous but did not diversify extensively until the late Cretaceous or early Paleogene about 32 million years after their origin. In a relatively small window of time, harpalines underwent rapid speciation. Harpalines have a relative high net diversification rate and increased cladogenesis in some regions of the clade. We did not see a significant decrease in diversification rate through time in the MCCR test, but a model of diversification with two shift points to lower diversification rates fit the harpaline lineage accumulation through time the best. CONCLUSIONS: Our results indicate harpalines are significantly more diverse and have higher diversification than their sistergroup. Instead of an immediate burst of explosive diversification, harpalines may have had a long "fuse" before major lineages diversified during the early Paleogene when other taxa such as mammals, birds, and some flowering plants were also rapidly diversifying.