Systematic analyses reveal uniqueness and origin of the CFEM domain in fungi

CFEM domain commonly occurs in fungal extracellular membrane proteins. To provide insights for understanding putative functions of CFEM, we investigate the evolutionary dynamics of CFEM domains by systematic comparative genomic analyses among diverse animals, plants, and more than 100 fungal species, which are representative across the entire group of fungi. We here show that CFEM domain is unique to fungi. Experiments using tissue culture demonstrate that the CFEM-containing ESTs in some plants originate from endophytic fungi. We also find that CFEM domain does not occur in all fungi. Its single origin dates to the most recent common ancestors of Ascomycota and Basidiomycota, instead of multiple origins. Although the length and architecture of CFEM domains are relatively conserved, the domain-number varies significantly among different fungal species. In general, pathogenic fungi have a larger number of domains compared to other species. Domain-expansion across fungal genomes appears to be driven by domain duplication and gene duplication via recombination. These findings generate a clear evolutionary trajectory of CFEM domains and provide novel insights into the functional exchange of CFEM-containing proteins from cell-surface components to mediators in host-pathogen interactions.

Allopolyploidization is not so simple: evidence from the origin of the tribe Cyprinini (Teleostei: Cypriniformes)

The identification of allopolyploidization events benefits from molecular dating and divergence assessments of progenitor genomes. Information on gene duplications only, either orthologs or paralogs, provides incomplete information. Analyses of mitochondrial DNA yield insights into matrilineal history, which may differ from patrilineal evolution. Two important food and pet fishes, the common carp (Cyprinus carpio) and goldfish (Carassius sp.), appear to have experienced allotetraploidization sometime from 12 to 20 million years ago (Ma). However, much work is necessary to detail the initial polyploidization event. Herein, we use this group of fishes as a model system to investigate competing scenarios for allopolyploidization. We analyze both the nuclear genes encoding growth hormone (GH), recombination activating protein 1 (RAG1) and HOXA2B gene, and the maternal heredited 12 concatenated mitochondrial protein-coding gene in 19 species of cyprinids and use two species in Balitoridae as outgroup taxa. Our analyses clarify the phylogenetic position of the paternal and maternal ancestors for the common carp and goldfish. The estimation of matrilineal divergence (10.71-12.42 Ma) is significantly younger than the dates of the parental ancestor divergedthat obtained by nuclear genes (16.62-19.64 Ma). Analyses of both genomes date the allopolyploidization event of the common ancestor of Cy. carpio and Ca sp. to about 10.71-12.42 Ma, which is most likely represented by maternal divergent time. The divergence of the two copies of the nuclear genes which was more ancient than the maternal markers might have been included the divergence of the progenitors' genome divergence when the allopolyploidization event occurred. Thus, the scenarios of allopolyploidzation for this group of fish can be suggested as the following: the matrilineal common ancestor of species in tribe Cyprinini might have doubled its genome by mating with a paternal ancestor in the subfamily Cyprininae, which was a sister-group that diverged around 4.20-8.93 Ma. Our work provides new evidence for the divergence dates of allopolyploidization within the Cyprinini, and documents the necessity of considering both matrilineal and patrilineal histories when investigating allopolyploidization.

Mitochondrial DNA variation of Nigerian domestic helmeted guinea fowl

We analyzed genetic diversity of 215 mitochondrial DNA (mtDNA) D-loop sequences from seven populations of domesticated helmeted guinea fowl (Numida meleagris) in Nigeria and compared that with results of samples collected in Kenya (n = 4) and China (n = 22). In total, 241 sequences were assigned to 22 distinct haplotypes. Haplotype diversity in Nigeria was 0.693 ± 0.022. The network grouped most matrilines into two main haplogroups: A and B. There was an absence of a geographic signal, and two haplotypes dominated across all locations with the exception of the Kebbi population in the northwest of the country; AMOVA also confirmed this observation (FST  = 0.035). The low genetic diversity may be a result of recent domestication, whereas the lack of maternal genetic structure likely suggests the extensive genetic intermixing within the country. Additionally, the differentiation of the Kebbi population may be due to a certain demographic history and/or artificial selection that shaped its haplotype profile. The current data do not permit us to make further conclusions; therefore, more research evidence from genetics and archaeology is still required.

DNA barcoding and the identification of tree frogs (Amphibia: Anura: Rhacophoridae)

The DNA barcoding gene COI (cytochrome c oxidase subunit I) effectively identifies many species. Herein, we barcoded 172 individuals from 37 species belonging to nine genera in Rhacophoridae to test if the gene serves equally well to identify species of tree frogs. Phenetic neighbor joining and phylogenetic Bayesian inference were used to construct phylogenetic trees, which resolved all nine genera as monophyletic taxa except for Rhacophorus, two new matrilines for Liuixalus, and Polypedates leucomystax species complex. Intraspecific genetic distances ranged from 0.000 to 0.119 and interspecific genetic distances ranged from 0.015 to 0.334. Within Rhacophorus and Kurixalus, the intra- and interspecific genetic distances did not reveal an obvious barcode gap. Notwithstanding, we found that COI sequences unambiguously identified rhacophorid species and helped to discover likely new cryptic species via the synthesis of genealogical relationships and divergence patterns. Our results supported that COI is an effective DNA barcoding marker for Rhacophoridae.

A genome draft of the legless anguid lizard, Ophisaurus gracilis

BACKGROUND

Transition from a lizard-like to a snake-like body form is one of the most important transformations in reptilian evolution. The increasing number of sequenced reptilian genomes is enabling a deeper understanding of vertebrate evolution, although the genetic basis of the loss of limbs in reptiles remains enigmatic. Here we report genome sequencing, assembly, and annotation for the Asian glass lizard Ophisaurus gracilis, a limbless lizard species with an elongated snake-like body form. Addition of this species to the genome repository will provide an excellent resource for studying the genetic basis of limb loss and trunk elongation.

FINDINGS

O. gracilis genome sequencing using the Illumina HiSeq2000 platform resulted in 274.20 Gbp of raw data that was filtered and assembled to a final size of 1.78 Gbp, comprising 6,717 scaffolds with N50 = 1.27 Mbp. Based on the k-mer estimated genome size of 1.71 Gbp, the assembly appears to be nearly 100% complete. A total of 19,513 protein-coding genes were predicted, and 884.06 Mbp of repeat sequences (approximately half of the genome) were annotated. The draft genome of O. gracilis has similar characteristics to both lizard and snake genomes.

CONCLUSIONS

We report the first genome of a lizard from the family Anguidae, O. gracilis. This supplements currently available genetic and genomic resources for amniote vertebrates, representing a major increase in comparative genome data available for squamate reptiles in particular.

DomeTree: a canonical toolkit for mitochondrial DNA analyses in domesticated animals

Mitochondrial DNA (mtDNA) is widely used in various genetic studies of domesticated animals. Many applications require comprehensive knowledge about the phylogeny of mtDNA variants. Herein, we provide the most up-to-date mtDNA phylogeny (i.e. haplogroup tree or matrilineal genealogy) and a standardized hierarchical haplogroup nomenclature system for domesticated cattle, dogs, goats, horses, pigs, sheep, yaks and chickens. These high-resolution mtDNA haplogroup trees based on 1240 complete or near-complete mtDNA genome sequences are available in open resource DomeTree (http://www.dometree.org). In addition, we offer the software MitoToolPy (http://www.mitotool.org/mp.html) to facilitate the mtDNA data analyses. We will continuously and regularly update DomeTree and MitoToolPy.

The complete mitochondrial genome of the black-tailed brush lizard Urosaurus nigricaudus (Reptilia, Squamata, Phrynosomatidae)

Previous studies using mitochondrial DNA (mtDNA) genes suggest the black-tailed brush lizard, Urosaurus nigricaudus, which is a small-sized lizard from the peninsula of Baja California, Mexico, has 4 deeply isolated mtDNA lineages with sequence divergence ranging from 4% to 11.2%. We present its complete mitochondrial genome. This genome is 17,298 bp long and comprises 2 rRNAs, 22 tRNAs, 13 protein-coding genes, 1 L-strand origin of replication and 1 control region. The overall nucleotide content is A = 34.2%; C = 26.8%; G = 13.5%; T = 25.5%. The gene organization and features agree with the general vertebrate organization and that found in other lizards. The control region is 1909 bp long and is located between tRNA^Pro and tRNA^Phe.

The endemic insular and peninsular species Chaetodipus spinatus (Mammalia, Heteromyidae) breaks patterns for Baja California

The Baja California peninsula is the second longest, most geographically isolated peninsula on Earth. Its physiography and the presence of many surrounding islands has facilitated studies of the underlying patterns and drivers of genetic structuring for a wide spectrum of organisms. Chaetodipus spinatus is endemic to the region and occurs on 12 associated islands, including 10 in the Gulf of California and two in the Pacific Ocean. This distribution makes it a model species for evaluating natural historical barriers. We test hypotheses associated with the relationship between the range of the species, patterns in other species, and its relationship to Pleistocene-Holocene climatic changes. We analyzed sequence data from mtDNA genes encoding cytochrome b (Cytb) and cytochrome c oxidase subunits I (COI) and III (COIII) in 26 populations including all 12 islands. The matrilineal genealogy, statistical parsimony network and Bayesian skyline plot indicated an origin of C. spinatus in the southern part of the peninsula. Our analyses detected several differences from the common pattern of peninsular animals: no mid-peninsula break exists, Isla Carmen hosts the most divergent population, the population on an ancient southern Midriff island does not differ from peninsular populations, and a mtDNA peninsular discordance occurs near Loreto.

The evolution of mitochondrial genomes in modern frogs (Neobatrachia): nonadaptive evolution of mitochondrial genome reorganization

Background

Although mitochondrial (mt) gene order is highly conserved among vertebrates, widespread gene rearrangements occur in anurans, especially in neobatrachians. Protein coding genes in the mitogenome experience adaptive or purifying selection, yet the role that selection plays on genomic reorganization remains unclear. We sequence the mitogenomes of three species of Glandirana and hot spots of gene rearrangements of 20 frog species to investigate the diversity of mitogenomic reorganization in the Neobatrachia. By combing these data with other mitogenomes in GenBank, we evaluate if selective pressures or functional constraints act on mitogenomic reorganization in the Neobatrachia. We also look for correlations between tRNA positions and codon usage.

Results

Gene organization in Glandirana was typical of neobatrachian mitogenomes except for the presence of pseudogene trnS (AGY). Surveyed ranids largely exhibited gene arrangements typical of neobatrachian mtDNA although some gene rearrangements occurred. The correlation between codon usage and tRNA positions in neobatrachians was weak, and did not increase after identifying recurrent rearrangements as revealed by basal neobatrachians. Codon usage and tRNA positions were not significantly correlated when considering tRNA gene duplications or losses. Change in number of tRNA gene copies, which was driven by genomic reorganization, did not influence codon usage bias. Nucleotide substitution rates and d_N/d_S ratios were higher in neobatrachian mitogenomes than in archaeobatrachians, but the rates of mitogenomic reorganization and mt nucleotide diversity were not significantly correlated.

Conclusions

No evidence suggests that adaptive selection drove the reorganization of neobatrachian mitogenomes. In contrast, protein-coding genes that function in metabolism showed evidence for purifying selection, and some functional constraints appear to act on the organization of rRNA and tRNA genes. As important nonadaptive forces, genetic drift and mutation pressure may drive the fixation and evolution of mitogenomic reorganizations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-691) contains supplementary material, which is available to authorized users.

DNA barcodes and species distribution models evaluate threats of global climate changes to genetic diversity: a case study from Nanorana parkeri (Anura: Dicroglossidae)

Anthropogenic global climate changes are one of the greatest threats to biodiversity. Distribution modeling can predict the effects of climate changes and potentially their effects on genetic diversity. DNA barcoding quickly identifies patterns of genetic diversity. As a case study, we use DNA barcodes and distribution models to predict threats under climate changes in the frog Nanorana parkeri, which is endemic to the Qinghai-Tibetan Plateau. Barcoding identifies major lineages W and E. Lineage W has a single origin in a refugium and Lineage E derives from three refugia. All refugia locate in river valleys and each greatly contributes to the current level of intraspecific genetic diversity. Species distribution models suggest that global climate changes will greatly influence N. parkeri, especially in the level of genetic diversity, because two former refugia will fail to provide suitable habitat. Our pipeline provides a novel application of DNA barcoding and has important implications for the conservation of biodiversity in southern areas of the Qinghai-Tibetan Plateau.

Genetic assessments and parentage analysis of captive Bolson tortoises (Gopherus flavomarginatus) inform their "rewilding" in New Mexico

The Bolson tortoise (Gopherus flavomarginatus) is the first species of extirpated megafauna to be repatriated into the United States. In September 2006, 30 individuals were translocated from Arizona to New Mexico with the long-term objective of restoring wild populations via captive propagation. We evaluated mtDNA sequences and allelic diversity among 11 microsatellite loci from the captive population and archived samples collected from wild individuals in Durango, Mexico (n = 28). Both populations exhibited very low genetic diversity and the captive population captured roughly 97.5% of the total wild diversity, making it a promising founder population. Genetic screening of other captive animals (n = 26) potentially suitable for reintroduction uncovered multiple hybrid G. flavomarginatus×G. polyphemus, which were ineligible for repatriation; only three of these individuals were verified as purebred G. flavomarginatus. We used these genetic data to inform mate pairing, reduce the potential for inbreeding and to monitor the maintenance of genetic diversity in the captive population. After six years of successful propagation, we analyzed the parentage of 241 hatchlings to assess the maintenance of genetic diversity. Not all adults contributed equally to successive generations. Most yearly cohorts of hatchlings failed to capture the diversity of the parental population. However, overlapping generations of tortoises helped to alleviate genetic loss because the entire six-year cohort of hatchlings contained the allelic diversity of the parental population. Polyandry and sperm storage occurred in the captives and future management strategies must consider such events.

The updated phylogenies of the phasianidae based on combined data of nuclear and mitochondrial DNA

The phylogenetic relationships of species in the Phasianidae, Order Galliformes, are the object of intensive study. However, convergent morphological evolution and rapid species radiation result in much ambiguity in the group. Further, matrilineal (mtDNA) genealogies conflict with trees based on nuclear DNA retrotransposable elements. Herein, we analyze 39 nearly complete mitochondrial genomes (three new) and up to seven nuclear DNA segments. We combine these multiple unlinked, more informative genetic markers to infer historical relationships of the major groups of phasianids. The nuclear DNA tree is largely congruent with the tree derived from mt genomes. However, branching orders of mt/nuclear trees largely conflict with those based on retrotransposons. For example, Gallus/Bambusicola/Francolinus forms the sister-group of Coturnix/Alectoris in the nuclear/mtDNA trees, yet the tree based on retrotransposable elements roots the former at the base of the tree and not with the latter. Further, while peafowls cluster with Gallus/Coturnix in the mt tree, they root at the base of the phasianids following Gallus in the tree based on retrotransposable elements. The conflicting branch orders in nuclear/mtDNA and retrotransposons-based trees in our study reveal the complex topology of the Phasianidae.

DNA barcoding of Vietnamese bent-toed geckos (Squamata: Gekkonidae: Cyrtodactylus) and the description of a new species

Species of bent-toed gecko (Cyrtodactylus) in Vietnam have been described at a rate of nearly four species per year since 2007 mostly based on morphological data. A tool that guides species delimitation will accelerate the rate of documentation, and at a time when the recognition of species greatly benefits conservation. We use DNA barcoding using COI (550 bp) to re-examine the levels of genetic divergence and taxonomic status of 21 described species of Vietnamese bent-toed geckos. Tree-based analyses resolve all sampled species and identify potential undescribed taxa. Kimura 2-parameter genetic distances between the described species average 21.0±4.2% and range from 4.3% to 28.7%. Further, our analyses discover two potentially new species from Vietnam, two from Laos and one from China. Herein we describe the new species Cyrtodactylus puhuensis sp. nov. from Vietnam on the basis of both genetics and morphology. Genetically, it differs from the remaining species by an average K2P distance of 24.0±1.8%. Morphologically, the new species is diagnosed by its medium-size (snout-vent length 79.24 mm and tail length 82.59 mm, for the single known individual), in having a series of moderately enlarged transverse subcaudals and a series of moderately enlarged femoral scales that extend from precloacal scales, in possessing femoral scales without pores, with males having five precloacal pores, and in exhibiting 8 supralabials, 10 infralabials, 23 narrow subdigital lamellae on its fourth toe, and 36 transverse ventrals.

Repeated functional convergent effects of NaV1.7 on acid insensitivity in hibernating mammals

Hibernating mammals need to be insensitive to acid in order to cope with conditions of high CO2; however, the molecular basis of acid tolerance remains largely unknown. The African naked mole-rat (Heterocephalus glaber) and hibernating mammals share similar environments and physiological features. In the naked mole-rat, acid insensitivity has been shown to be conferred by the functional motif of the sodium ion channel NaV1.7. There is now an opportunity to evaluate acid insensitivity in other taxa. In this study, we tested for functional convergence of NaV1.7 in 71 species of mammals, including 22 species that hibernate. Our analyses revealed a functional convergence of amino acid sequences, which occurred at least six times independently in mammals that hibernate. Evolutionary analyses determined that the convergence results from both parallel and divergent evolution of residues in the functional motif. Our findings not only identify the functional molecules responsible for acid insensitivity in hibernating mammals, but also open new avenues to elucidate the molecular underpinnings of acid insensitivity in mammals.

Conservation genetics of the endangered Spotted Turtle (Clemmys guttata) illustrate the risks of “bottleneck tests”

Studies of population genetics in turtles have suggested that turtles do not experience genetic impacts of bottlenecks as strongly as expected. However, recent studies cast doubt on two commonly used tests implemented in the program BOTTLENECK, suggesting that these findings should be re-evaluated. The Spotted Turtle (Clemmys guttata (Schneider, 1792)) is endangered both globally and within Canada, but genetic data required to develop effective recovery strategies are unavailable. Here, we conducted the first study of population genetic structure in C. guttata. We then used multiple small populations of C. guttata as replicates to test whether the commonly used program BOTTLENECK could detect the genetic signature of bottlenecks in our study populations, which are all thought to have experienced significant declines in the past 2–3 generations (75 years). Turtles (n = 256) were genotyped at 11 microsatellite loci. A suite of Bayesian population genetics analyses and a principal coordinates analysis identified a minimum of 6 distinct genetic populations and a maximum of 10 differentiated subpopulations across the sampled Canadian range of C. guttata, which corresponded to demographically independent units. BOTTLENECK failed to detect population declines. A literature review found that bottleneck tests in 17 of 18 previous genetic studies of tortoises and freshwater turtles were based on suboptimal sampling, potentially confounding their results. High retention of genetic diversity (allelic richness and heterozygosity) in isolated populations of C. guttata and other turtle species is encouraging for species recovery, but conclusions about the prevalence of genetic bottlenecks in such populations should be re-examined.

Evolutionary perspective on hepatitis B virus with an expanded sampling strategy

To investigate the role hepatitis B e antigen (HBeAg) plays in the evolution of hepatitis B virus (HBV), we sequenced the basic core promoter (BCP) and precore (preC) regions of 348 clones total from ten HBV Chinese patients. Eleven mutations were more frequent in HBeAg-negative patients than in HBeAg-positive patients. Further, the sequencing of dozens of variants was found to be necessary to obtain mutation profiles. Phylogenetic and median-joining network analyses suggested that variants from each patient had a single common ancestor (monophyly). Higher haplotype and nucleotide diversities were identified in HBeAg-negative patients. Analysis of dN/dS suggested that viruses experiencing a stronger immune response had lower haplotype diversity. Because HBeAg seroconversion was associated with viral diversity it served as an indicator of HBV evolution. Significantly, this study indicated a larger sampling of variants from each patient was required to understand effectively the properties of HBV.

Obligate mutualism within a host drives the extreme specialization of a fig wasp genome

Background

Fig pollinating wasps form obligate symbioses with their fig hosts. This mutualism arose approximately 75 million years ago. Unlike many other intimate symbioses, which involve vertical transmission of symbionts to host offspring, female fig wasps fly great distances to transfer horizontally between hosts. In contrast, male wasps are wingless and cannot disperse. Symbionts that keep intimate contact with their hosts often show genome reduction, but it is not clear if the wide dispersal of female fig wasps will counteract this general tendency. We sequenced the genome of the fig wasp Ceratosolen solmsi to address this question.

Results

The genome size of the fig wasp C. solmsi is typical of insects, but has undergone dramatic reductions of gene families involved in environmental sensing and detoxification. The streamlined chemosensory ability reflects the overwhelming importance of females finding trees of their only host species, Ficus hispida, during their fleeting adult lives. Despite long-distance dispersal, little need exists for detoxification or environmental protection because fig wasps spend nearly all of their lives inside a largely benign host. Analyses of transcriptomes in females and males at four key life stages reveal that the extreme anatomical sexual dimorphism of fig wasps may result from a strong bias in sex-differential gene expression.

Conclusions

Our comparison of the C. solmsi genome with other insects provides new insights into the evolution of obligate mutualism. The draft genome of the fig wasp, and transcriptomic comparisons between both sexes at four different life stages, provide insights into the molecular basis for the extreme anatomical sexual dimorphism of this species.

Phylogeny of the Cyrtodactylus irregularis species complex (Squamata: Gekkonidae) from Vietnam with the description of two new species

The number of described species of bent-toed geckos of the Cyrtodactylus irregularis species complex in Vietnam has increased from one to eight in the last six years. We combined morphological and molecular analyses to explore phylogenetic relationships among all described species in the group. The phylogeny required the description of two new species, Cyrtodactylus phuocbinhensis sp. nov. and Cyrtodactylus taynguyenensis sp. nov. Further, the tree resolved two additional undescribed clades that may also be new species. The species C. bugiamapensis and C. ziegleri were found to require redefinition. Cyrtodactylus phuocbinhensis sp. nov. is characterized by a series of enlarged femoral scales separated from preanal scales while Cyrtodactylus taynguyenensis sp. nov. does not possess enlarged femoral scales. Both new species are distinguished from other congeners by a combination of the following characters: small subcaudal scales, not transversely enlarged; presence (C. phuocbinhensis sp. nov.) or absence (C. taynguyenensis sp. nov.) of enlarged femoral scales; number of preanal pores; and dorsal pattern. Genetic distances between described species and new species were 16.5% and 2.0% in COI and RPL35, respectively, for C. phuocbinhensis sp. nov., and these distances were 18.8% and 2.2% for C. taynguyenensis sp. nov., respectively.

The evolution and expression of the moth visual opsin family

Because visual genes likely evolved in response to their ambient photic environment, the dichotomy between closely related nocturnal moths and diurnal butterflies forms an ideal basis for investigating their evolution. To investigate whether the visual genes of moths are associated with nocturnal dim-light environments or not, we cloned long-wavelength (R), blue (B) and ultraviolet (UV) opsin genes from 12 species of wild-captured moths and examined their evolutionary functions. Strong purifying selection appeared to constrain the functions of the genes. Dark-treatment altered the levels of mRNA expression in Helicoverpa armigera such that R and UV opsins were up-regulated after dark-treatment, the latter faster than the former. In contrast, B opsins were not significantly up-regulated. Diel changes of opsin mRNA levels in both wild-captured and lab-reared individuals showed no significant fluctuation within the same group. However, the former group had significantly elevated levels of expression compared with the latter. Consequently, environmental conditions appeared to affect the patterns of expression. These findings and the proportional expression of opsins suggested that moths potentially possessed color vision and the visual system played a more important role in the ecology of moths than previously appreciated. This aspect did not differ much from that of diurnal butterflies.