Ancestral Character Reconstruction Research Articles

Phylogeographical and evolutionary history of variola major virus; a question of timescales?

Aim of this study was to reconstruct the phylogeography of variola virus (VARV) in the XX century, using 47 VARV whole genome sequences available in public databases, through two different methods for ancestral character reconstruction: a frequently used Bayesian framework and a fast maximum-likelihood (ML) based method. The substitution rate of the whole VARV genome was estimated to be between 6.7×10-6 and 1.1×10-5 substitutions/site/year. Both ML and Bayesian methods gave similar trees topology, showing two distinct monophyletic groups: one (known as P1) including the great part of variola major and the second (P2) including West African and American (variola minor) isolates and close evolutionary rate estimations, between 6.73×10-6 and 1.1×10-5 for the whole genome. The phylogeographical reconstruction of P1 suggested that the common ancestor of the variola major circulating in the Old World between the 1940s and the 1970s most probably originated in the Far East in the first decades of the XX century, and then spread to Indian subcontinent in the 1920s. India represented a center of further spread of VARV to eastern Africa in the 1940s and to the Middle East in the 1960s. The phylogeographic scenario obtained by the maximum-likelihood based method was congruent with that obtained by Bayesian framework, but the analysis was faster indicating the usefulness of this method in the analyses of large viral genomes. Our results may help to explain the controversial reconstructions of the history of VARV obtained using long or short timescale for calibration.

A well-supported nuclear phylogeny of Poaceae and implications for the evolution of C4 photosynthesis

Poaceae (the grasses) includes rice, maize, wheat, and other crops, and is the most economically important angiosperm family. Poaceae is also one of the largest plant families, consisting of over 11 000 species with a global distribution that contributes to diverse ecosystems. Poaceae species are classified into 12 subfamilies, with generally strong phylogenetic support for their monophyly. However, many relationships within subfamilies, among tribes and/or subtribes, remain uncertain. To better resolve the Poaceae phylogeny, we generated 342 transcriptomic and seven genomic datasets; these were combined with other genomic and transcriptomic datasets to provide sequences for 357 Poaceae species in 231 genera, representing 45 tribes and all 12 subfamilies. Over 1200 low-copy nuclear genes were retrieved from these datasets, with several subsets obtained using additional criteria, and used for coalescent analyses to reconstruct a Poaceae phylogeny. Our results strongly support the monophyly of 11 subfamilies; however, the subfamily Puelioideae was separated into two non-sister clades, one for each of the two previously defined tribes, supporting a hypothesis that places each tribe in a separate subfamily. Molecular clock analyses estimated the crown age of Poaceae to be ∼101 million years old. Ancestral character reconstruction of C3/C4 photosynthesis supports the hypothesis of multiple independent origins of C4 photosynthesis. These origins are further supported by phylogenetic analysis of the ppc gene family that encodes the phosphoenolpyruvate carboxylase, which suggests that members of three paralogous subclades (ppc-aL1a, ppc-aL1b, and ppc-B2) were recruited as functional C4ppc genes. This study provides valuable resources and a robust phylogenetic framework for evolutionary analyses of the grass family.

Genome size evolution of the extant lycophytes and ferns

Ferns and lycophytes have remarkably large genomes. However, little is known about how their genome size evolved in fern lineages. To explore the origins and evolution of chromosome numbers and genome size in ferns, we used flow cytometry to measure the genomes of 240 species (255 samples) of extant ferns and lycophytes comprising 27 families and 72 genera, of which 228 species (242 samples) represent new reports. We analyzed correlations among genome size, spore size, chromosomal features, phylogeny, and habitat type preference within a phylogenetic framework. We also applied ANOVA and multinomial logistic regression analysis to preference of habitat type and genome size. Using the phylogeny, we conducted ancestral character reconstruction for habitat types and tested whether genome size changes simultaneously with shifts in habitat preference. We found that 2C values had weak phylogenetic signal, whereas the base number of chromosomes (x) had a strong phylogenetic signal. Furthermore, our analyses revealed a positive correlation between genome size and chromosome traits, indicating that the base number of chromosomes (x), chromosome size, and polyploidization may be primary contributors to genome expansion in ferns and lycophytes. Genome sizes in different habitat types varied significantly and were significantly correlated with habitat types; specifically, multinomial logistic regression indicated that species with larger 2C values were more likely to be epiphytes. Terrestrial habitat is inferred to be ancestral for both extant ferns and lycophytes, whereas transitions to other habitat types occurred as the major clades emerged. Shifts in habitat types appear be followed by periods of genomic stability. Based on these results, we inferred that habitat type changes and multiple whole-genome duplications have contributed to the formation of large genomes of ferns and their allies during their evolutionary history.

An Extraordinary Rosette and Resurrection New Spikemoss, Selaginella iridescens (Selaginellaceae) from Hainan Island, China

Selaginella iridescens (Selaginellaceae), a new species of spikemoss from Hainan Island, China, is described based on morphological and molecular evidence. The new species morphologically resembles S. pseudotamariscina from Vietnam which was recently recognized and segregated from S. tamariscina and S. pulvinata. Some characters of S. iridescens, including the iridescent leaves and long arista (0.35–1.20 mm long), the widely membranaceous margin of axillary leaves (ca. 2/3), dorsal leaves sulcate extending to the top, and slender main stem, dorsoventrally prostrate, distinguish it from S. pseudotamariscina. Phylogenetic results based on rbcL of 28 Selaginella species indicate that S. iridescens is sister to S. pseudotamariscina, and distant from the S. tamariscina-S.pulvinata clade. The ancestral character reconstruction result reveals that the rosette is apomorphic and has evolved independently at least six times in Selaginella.

Cuban history of CRF19 recombinant subtype of HIV-1.

CRF19 is a recombinant form of HIV-1 subtypes D, A1 and G, which was first sampled in Cuba in 1999, but was already present there in 1980s. CRF19 was reported almost uniquely in Cuba, where it accounts for ∼25% of new HIV-positive patients and causes rapid progression to AIDS (∼3 years).We analyzed a large data set comprising ∼350 pol and env sequences sampled in Cuba over the last 15 years and ∼350 from Los Alamos database. This data set contained both CRF19 (∼315), and A1, D and G sequences. We performed and combined analyses for the three A1, G and D regions, using fast maximum likelihood approaches, including: (1) phylogeny reconstruction, (2) spatio-temporal analysis of the virus spread, and ancestral character reconstruction for (3) transmission mode and (4) drug resistance mutations (DRMs). We verified these results with a Bayesian approach. This allowed us to acquire new insights on the CRF19 origin and transmission patterns. We showed that CRF19 recombined between 1966 and 1977, most likely in Cuban community stationed in Congo region. We further investigated CRF19 spread on the Cuban province level, and discovered that the epidemic started in 1970s, most probably in Villa Clara, that it was at first carried by heterosexual transmissions, and then quickly spread in the 1980s within the “men having sex with men” (MSM) community, with multiple transmissions back to heterosexuals. The analysis of the transmission patterns of common DRMs found very few resistance transmission clusters.Our results show a very early introduction of CRF19 in Cuba, which could explain its local epidemiological success. Ignited by a major founder event, the epidemic then followed a similar pattern as other subtypes and CRFs in Cuba. The reason for the short time to AIDS remains to be understood and requires specific surveillance, in Cuba and elsewhere.

Looking at the past to infer into the future: Thermal traits track environmental change in Liolaemidae.

The diversity of habitats generated by the Andes uplift resulted a mosaic of heterogeneous environments in South America for species to evolve a variety of ecological and physiological specializations. Species in the lizard family Liolaemidae occupy a myriad of habitats in the Andes. Here, we analyze the tempo and mode of evolution in the thermal biology of liolaemids. We assessed whether there is evidence of local adaptation (lability) or conservatism (stasis) in thermal traits. We tested the hypothesis that abiotic factors (e.g., geography, climate) rather than intrinsic factors (egg-laying [oviparous] or live-bearing [viviparous], substrate affinity) explain variation in field active body temperature (Tb ), preferred temperature (Tp ), hours of restriction of activity, and potential hours of activity. Although most traits exhibited high phylogenetic signal, we found variation in thermal biology was shaped by geography, climate, and ecological diversity. Ancestral character reconstruction showed shifts in Tb tracked environmental change in the past ∼20,000 years. Thermal preference is 3°C higher than Tb , yet exhibited a lower rate of evolution than Tb and air temperature. Viviparous Liolaemus have lower Tb s than oviparous species, whereas Tp is high for both modes of reproduction, a key difference that results in a thermal buffer for viviparous species to cope with global warming. The rapid increase in environmental temperatures expected in the next 50-80 years in combination with anthropogenic loss of habitats are projected to cause extirpations and extinctions in oviparous species.

A new infrageneric classification of the pantropical genusChamaecrista(Fabaceae: Caesalpinioideae) based on a comprehensive molecular phylogenetic analysis and morphology

AbstractChamaecrista with > 330 species, six sections, three subsections and 39 series has had a long and complex taxonomic history. The genus is monophyletic, but most of its traditional infrageneric categories are not. To test the monophyly of sections, subsections and series of Chamaecrista, we used two molecular phylogenetic approaches. The first (Broad) based on two DNA regions (ITS and trnL-F) includes a comprehensive sampling of Chamaecrista spp. and infrageneric taxa. The second (Multilocus) is based on four molecular regions (ITS, ETS, trnL-F and trnE-T) for a smaller but representative sampling. We performed ancestral character reconstructions to identify morphological characters that could serve as synapomorphies for major clades. Both molecular approaches support Chamaecrista and sections Apoucouita, Grimaldia and Xerocalyx as monophyletic, but sections Chamaecrista, Caliciopsis and Absus and most of the series are not monophyletic. The four main clades recovered are all characterized by a combination of morphological characters: a clade of tree species with cauliflorous inflorescences (including species of section Apoucouita); a mostly Brazilian campo rupestre clade (including all species of subsections Adenophyllum, Baseophyllum and Otophyllum); a clade of mostly herbaceous/shrubby species with solitary flowers or fascicles (including sections Chamaecrista, Caliciopsis and Xerocalyx and extra-American species) and a clade (with three main subclades) of species with viscous indumentum (including section Grimaldia and section Absus subsection Absus). We propose a new infrageneric classification for Chamaecrista supported by molecular phylogenetic analyses and morphology, recognizing the four main clades as sections Apoucouita, Baseophyllum, Chamaecrista and Absus, the last with three subsections (Absus, Viscosa and Zygophyllum), but we do not recognize any previously circumscribed series. Our taxonomic treatment includes descriptions of and a key to the newly defined infrageneric taxa and an updated species list for the genus under the new classification.

Nuclear phylotranscriptomics and phylogenomics support numerous polyploidization events and hypotheses for the evolution of rhizobial nitrogen-fixing symbiosis in Fabaceae

Fabaceae are the third largest angiosperm family, with 765 genera and ∼19 500 species. They are important both economically and ecologically, and global Fabaceae crops are intensively studied in part for their nitrogen-fixing ability. However, resolution of the intrasubfamilial Fabaceae phylogeny and divergence times has remained elusive, precluding a reconstruction of the evolutionary history of symbiotic nitrogen fixation in Fabaceae. Here, we report a highly resolved phylogeny using >1500 nuclear genes from newly sequenced transcriptomes and genomes of 391 species, along with other datasets, for a total of 463 legumes spanning all 6 subfamilies and 333 of 765 genera. The subfamilies are maximally supported as monophyletic. The clade comprising subfamilies Cercidoideae and Detarioideae is sister to the remaining legumes, and Duparquetioideae and Dialioideae are successive sisters to the clade of Papilionoideae and Caesalpinioideae. Molecular clock estimation revealed an early radiation of subfamilies near the K/Pg boundary, marked by mass extinction, and subsequent divergence of most tribe-level clades within ∼15 million years. Phylogenomic analyses of thousands of gene families support 28 proposed putative whole-genome duplication/whole-genome triplication events across Fabaceae, including those at the ancestors of Fabaceae and five of the subfamilies, and further analyses supported the Fabaceae ancestral polyploidy. The evolution of rhizobial nitrogen-fixing nodulation in Fabaceae was probed by ancestral character reconstruction and phylogenetic analyses of related gene families and the results support the hypotheses of one or two switch(es) to rhizobial nodulation followed by multiple losses. Collectively, these results provide a foundation for further morphological and functional evolutionary analyses across Fabaceae.

Niche partitioning among three snail-eating snakes revealed by dentition asymmetry and prey specialisation.

The level of dentition asymmetry in snail-eating snakes may reflect their prey choice and feeding efficiency on asymmetric land snails. The three species of Pareas snakes (Squamata: Pareidae) in Taiwan, which form partially sympatric distribution on the island, provide a potential case to test the hypothesis of niche partitioning and character displacement with regard to dentition asymmetry and specialisation in feeding behaviour. In this study, behavioural experiments confirmed that P. formosensis feeds exclusively on slugs, whereas P. atayal and P. komaii consumed both. However, P. atayal more efficiently preys on land snails than P. komaii, exhibiting a shorter handling time and fewer mandibular retractions. Micro-CT and ancestral character reconstruction demonstrated the lowest asymmetry in P. formosensis (the slug specialist), the highest dentition asymmetry in P. atayal (the land snail specialist) and flexibility in P. komaii (the niche switcher): increased dentition asymmetry when sympatrically distributed with the slug eater (character displacement), and decreased asymmetry when living alone (ecological niche release). Ecological niche modelling showed that the distribution of P. formosensis is associated with the presence of slugs, while that of P. atayal could be explained by the land snails. Combining the results from morphology, phylogeny, behavioural experiments and ecological niche modelling, we showed that competition in the sympatric region might have facilitated character displacement among congeners, while the absence of competition in allopatric region has led to ecological niche release.

Photosynthesis metabolism in the Compositae: Current knowledge and future directions

AbstractCompositae accounts for ca. 10% of all flowering plants, being found in nearly all types of habitats worldwide. The family is particularly diverse in markedly dry and seasonal habitats, raising questions about what processes led to diversification in these challenging environments. The presence of C4 and CAM metabolism in some Compositae taxa has been known since the 1970s, and although some of these taxa have been extensively studied, such as Flaveria, there has been no systematic effort in gathering and reviewing data about photosynthesis metabolism in the family in the last few decades. In the present paper, we gathered data from more than 50 articles, spanning several different methods, raising information on photosynthesis metabolism for more than 400 Compositae species. We also present newly acquired carbon isotope data for 66 species, mostly from previously unsampled tribes. We analyzed photosynthesis metabolism in the tribe Tageteae using a phylogeny and ancestral character reconstruction. C4 photosynthesis in Compositae is restricted to two tribes: Tageteae, where it had two independent origins, and Coreopsideae, where it seems to have arisen only once. CAM metabolism is found in succulent members of tribe Senecioneae, and in Astereae and Eupatorieae, with an unknown number of evolutionary origins. We discuss the evolutionary implications of the results and propose future directions for studying photosynthesis metabolism in the family.

Molecular phylogeny and biogeography ofPabstiella(Pleurothallidinae: Orchidaceae) highlight the importance of the Atlantic Rainforest for speciation in the genus

AbstractPabstiella consists of c. 130 epiphytic species in the Neotropics. We present a phylogenetic analysis based on nrITS, matK and trnH-psbA sequences from 59 species of the genus and 40 Pleurothallidinae and two Laeliinae and one Bletiinae as an outgroup, using maximum likelihood, Bayesian inference and maximum parsimony. We also performed molecular dating, biogeographical analyses and ancestral morphological character reconstruction. Our results confirm the monophyly of Pabstiella with strong support. Ten clades are inferred and are herein proposed as sections. Pabstiella originated in the Andes and the Atlantic Rainforest in the Late Miocene (c. 7.93 Mya) in an epoch when these biomes were probably connected. A main vicariance event divided an early-diverging lineage that inhabited the Andes from an Atlantic Rainforest lineage that diversified in this region during the Pliocene and Pleistocene, mainly in the Serra do Mar in south-eastern Brazil. Our findings also suggest that the Atlantic Rainforest may have played an important role in the origin of subtribe Pleurothallidinae. The morphological character reconstruction showed high levels of homoplasy, with few recognized synapomorphies associated with stems and petals. Other characters related to the habit and stems were identified as important in the evolutionary history of the genus.

Phylogeography of the marine pathogen, Vibrio vulnificus, revealed the ancestral scenarios of its evolution.

Vibrio vulnificus is the leading cause of seafood‐associated deaths worldwide. Despite the growing knowledge about the population structure of V. vulnificus, the evolutionary history and the ancestral relationships of strains isolated from various regions around the world have not been determined. Using the largest collection of sequence and isolate data of V. vulnificus to date, we applied ancestral character reconstruction to study the phylogeography of V. vulnificus. Multilocus sequence typing data from 10 housekeeping genes were used for the inference of ancestral states and reconstruction of the evolutionary history. The findings showed that the common ancestor of all V. vulnificus populations originated from East Asia, and later evolved into two main clusters that spread with time and eventually evolved into distinct populations in different parts of the world. While we found no meaningful insights concerning the evolution of V. vulnificus populations in the Middle East; however, we were able to reconstruct the ancestral scenarios of its evolution in East Asia, North America, and Western Europe.

Comparative phylogenomics reveal complex evolution of life history strategies in a clade of bivalves with parasitic larvae (Bivalvia: Unionoida: Ambleminae).

Freshwater mussels are a species-rich group with biodiversity patterns strongly shaped by a life history strategy that includes an obligate parasitic larval stage. In this study, we set out to reconstruct the life history evolution and systematics in a clade of freshwater mussels adapted to parasitizing a molluscivorous host fish. Anchored hybrid enrichment and ancestral character reconstruction revealed a complex pattern of life history evolution with host switching and multiple instances of convergence, including reduction in size of larvae, increased fecundity, and growth during encapsulation. Our phylogenomic analyses also recovered non-monophyly of taxa exhibiting multiple traits used as the basis for previous taxonomic hypotheses. Taxa with axe-head shaped glochidia were resolved as paraphyletic, but our results strongly suggest the complex morphology is an adaptation to reduce larval size, with reduction in size further accentuated in taxa previously assigned to Leptodea. To more accurately reflect the evolutionary history of this group, we make multiple systematic changes, including the description of a new genus, Atlanticoncha gen. nov., and the synonymy of the genus Leptodea under Potamilus. Our findings contribute to the growing body of literature showing that cladistic hypotheses based solely on morphological characters, including larval morphology, can be flawed in freshwater mussels.

Comparative transcriptomics of the venoms of continental and insular radiations of West African cones.

The transcriptomes of the venom glands of 13 closely related species of vermivorous cones endemic to West Africa from genera Africonus and Varioconus were sequenced and venom repertoires compared within a phylogenetic framework using one Kalloconus species as outgroup. The total number of conotoxin precursors per species varied between 108 and 221. Individuals of the same species shared about one-fourth of the total conotoxin precursors. The number of common sequences was drastically reduced in the pairwise comparisons between closely related species, and the phylogenetical signal was totally eroded at the inter-generic level (no sequence was identified as shared derived), due to the intrinsic high variability of these secreted peptides. A common set of four conotoxin precursor superfamilies (T, O1, O2 and M) was expanded in all studied cone species, and thus, they are considered the basic venom toolkit for hunting and defense in the West African vermivorous cone snails. Maximum-likelihood ancestral character reconstructions inferred shared conotoxin precursors preferentially at internal nodes close to the tips of the phylogeny (between individuals and between closely related species) as well as in the common ancestor of Varioconus. Besides the common toolkit, the two genera showed significantly distinct catalogues of conotoxin precursors in terms of type of superfamilies present and the abundance of members per superfamily, but had similar relative expression levels indicating functional convergence. Differential expression comparisons between vermivorous and piscivorous cones highlighted the importance of the A and S superfamilies for fish hunting and defense.

Molecular phylogenetic analyses reveal a deep dichotomy in the conifer-inhabiting genus Trisetacus (Eriophyoidea: Nalepellidae), with the two lineages differing in their female genital morphology and host associations.

We analyzed the phylogenetic relationships of the genus Trisetacus using two genes [cytochrome c oxidase subunit I (COI) and D1-D2 region of 28S rDNA (D1-D2 28S)], a representive taxon sampling (nearly 40% of known diversity), and a large set of close and distant outgroups. Our analyses suggest the presence of a dichotomy between Trisetacus associated with Cupressaceae and Pinaceae. The following smaller molecular clades were found: Pin-1 (bud mites, twig sheath mites, bark gall mites, and endoparasitic mites from pinaceans), Pin-2 (needle sheath mites from pines), Pin-2a (putative Nearctic group of needle sheath mites), Pin-2b (putative Palearctic group of needle sheath mites), Cup-1 and 2 (bud, cone, seed mites and mites living under bark scales from cupressaceans). The monophyly of the recently proposed subgenus Brevithecus nested within clade Cup-2 was confirmed. Ancestral character reconstruction analyses recovered: (1) Pinaceae as the ancestral hosts of Nalepellidae and Trisetacus, (2) repetitive reductions of the spermathecal tube independently occurred in two lineages of Trisetacus from Cupressaceae, and (3) several mite habitats on host (galls, cones, twig sheaths, seeds, inside leaves, and under scales) are evolutionarily derived states, whereas living in buds or needle sheaths are ancestral states for Trisetacus clades Cup and Pin. Using confocal microscopy, we identified six basic types of the female internal genitalia of Trisetacus based on shapes of the spermatheca and spermathecal tube. These genitalic types are strongly correlated with lineages recovered by molecular phylogenetic analyses, suggesting that the female genital morphology is both evolutionarily conserved and is a factor influencing macroevolutionary patterns in this group of mites.

Parallel Evolution of Bower-Building Behavior in Two Groups of Bowerbirds Suggested by Phylogenomics.

The bowerbirds in New Guinea and Australia include species that build the largest and perhaps most elaborately decorated constructions outside of humans. The males use these courtship bowers, along with their displays, to attract females. In these species, the mating system is polygynous and the females alone incubate and feed the nestlings. The bowerbirds also include 10 species of the socially monogamous catbirds in which the male participates in most aspects of raising the young. How the bower-building behavior evolved has remained poorly understood, as no comprehensive phylogeny exists for the family. It has been assumed that the monogamous catbird clade is sister to all polygynous species. We here test this hypothesis using a newly developed pipeline for obtaining homologous alignments of thousands of exonic and intronic regions from genomic data to build a phylogeny. Our well-supported species tree shows that the polygynous, bower-building species are not monophyletic. The result suggests either that bower-building behavior is an ancestral condition in the family that was secondarily lost in the catbirds, or that it has arisen in parallel in two lineages of bowerbirds. We favor the latter hypothesis based on an ancestral character reconstruction showing that polygyny but not bower-building is ancestral in bowerbirds, and on the observation that Scenopoeetes dentirostris, the sister species to one of the bower-building clades, does not build a proper bower but constructs a court for male display. This species is also sexually monomorphic in plumage despite having a polygynous mating system. We argue that the relatively stable tropical and subtropical forest environment in combination with low predator pressure and rich food access (mostly fruit) facilitated the evolution of these unique life-history traits. [Adaptive radiation; bowerbirds; mating system, sexual selection; whole genome sequencing.]

Coping with the extremes: comparative osteology of the tepui-associated toad Oreophrynella and its bearing on the evolution of osteological novelties in the genus

AbstractThe only study of the osteology of the toad genus Oreophrynella dates back to 1971 and was based on a single species. Here, we use high-resolution X-ray microcomputed tomography to analyse the osteology of all nine described Oreophrynella species, which are compared with representatives of other bufonid lineages. Oreophrynella is unique among bufonids in having opposable digits. Osteological synapomorphies confirmed for the genus are as follows: presence of parietal fontanelles and exposed frontoparietal fontanelle, absence of quadratojugal, five presacral vertebrae, distally enlarged terminal phalanges and urostyle greatly expanded into flanges. Ancestral character reconstruction indicates that arboreal habits in some Oreophrynella species are likely to have evolved after the evolution of opposable digits. Opposable digits, in combination with an extension of the interdigital integument and the relative length/orientation of the digits, are likely to be adaptations to facilitate life on rocky tepui summits and an exaptation to arboreality. Cranial simplification in Oreophrynella, in the form of cranial fontanelles and absence of the quadratojugal, is possibly driven by a reduction of developmental costs, increase in flexibility and reduction of body weight. Cranial simplification combined with the shortening of the vertebral column and the shift towards a partly firmisternal girdle might be adaptations to the peculiar tumbling behaviour displayed by Oreophrynella.

Mitogenomic phylogeny of the Naticidae (Gastropoda: Littorinimorpha) reveals monophyly of the Polinicinae

AbstractThe Naticidae is a species‐rich family of predatory marine gastropods with substantial interspecific morphological diversity. The classification of the Naticidae has been traditionally based on morphology data, but the phylogenetic relationships within the family are debated due to conflicting molecular results, especially regarding the monophyly of subfamilies Polinicinae and Naticinae. To further resolve the phylogenetic controversies within the Naticidae, we undertake a phylogenetic approach using 14 newly sequenced complete or nearly complete (only lacking a control region) mitochondrial genomes. Both the maximum likelihood and Bayesian inference analyses supported monophyly of the Polinicinae, but paraphyly of the Naticinae due to the placement of the enigmatic genus Notocochlis. The ancestral character reconstruction suggests that the operculum, a character that currently defines the two subfamilies, evolved from an ancestor with a calcareous operculum in the evolutionary history of naticids. In addition, the chronogram estimates that naticids was originated in late Triassic (about 227 million years ago), consistent with previous hypotheses. Our study highlights the importance of using complete mitochondrial genomes while reconstructing phylogenetic relationships within the Naticidae. The evolution scenario of the naticid operculum contributes new insights into the classification of Naticidae.

Phylogenetic relationships between fungus-associated Neotropical species of the genera Hirtodrosophila, Mycodrosophila and Zygothrica (Diptera, Drosophilidae), with insights into the evolution of breeding sites usage

The Neotropical region harbors an astonishing diversity of species, but still encompasses the least studied biogeographic region of the world. These properties apply for different taxonomic groups, and can be exemplified by drosophilids. In fact, high levels of cryptic diversity have recently been discovered for Neotropical species of the Zygothrica genus group, but relationships among these species, or them and other Drosophilidae species still remains to be addressed. Therefore, the aim of this study was to evaluate the phylogenetic relationships between fungus-associated Neotropical species of the genera Hirtodrosophila, Mycodrosophila and Zygothrica, which together with Paramycodrosophila and Paraliodrosophila compose the Zygothrica genus group. For this, fragments of the mitochondrial cytochrome oxidase subunits I (COI) and II (COII) genes, and the nuclear alpha methyldopa (Amd) and dopa decarboxylase (Ddc) genes were newly characterized for 43 Neotropical specimens of fungus-associated drosophilids, and analyzed in the context of 51 additional Drosophilinae sequences plus one Steganinae outgroup. Based on the resulting phylogeny, the evolution of breeding sites usage was also evaluated through ancestral character reconstructions. Our results revealed the Zygothrica genus group as a monophyletic lineage of Drosophila that branches after the subgenera Sophophora and Drosophila. Within this lineage, Mycodrosophila species seem to encompass the early offshoot, followed by a grade of Hirtodrosophila species, with derived branches mostly occupied by representatives of Zygothrica. This genus, in particular, was subdivided into five major clades, two of which include species of Hirtodrosophila, whose generic status needs to be reevatuated. According to our results, the use of fungi as breeding sites encompasses a symplesiomorphy for the Zygothrica genus group, since one of the recovered clades is currently specialized in using flowers as breeding sites whereas a sole species presents a reversal to the use of fruits of a plant of Gentianales. So, in general, this study supports the paraphyly of Drosophila in relation to fungus-associated Neotropical species of Drosophilidae, providing the first molecular insights into the phylogenetic patterns related to the evolution of this diverse group of species and some of its characteristic traits.

Studying 3D genome evolution using genomic sequence.

The three dimensions (3D) genome is essential to numerous key processes such as the regulation of gene expression and the replication-timing program. In vertebrates, chromatin looping is often mediated by CTCF, and marked by CTCF motif pairs in convergent orientation. Comparative high-throughput sequencing technique (Hi-C) recently revealed that chromatin looping evolves across species. However, Hi-C experiments are complex and costly, which currently limits their use for evolutionary studies over a large number of species. Here, we propose a novel approach to study the 3D genome evolution in vertebrates using the genomic sequence only, e.g. without the need for Hi-C data. The approach is simple and relies on comparing the distances between convergent and divergent CTCF motifs by computing a ratio we named the 3D ratio or '3DR'. We show that 3DR is a powerful statistic to detect CTCF looping encoded in the human genome sequence, thus reflecting strong evolutionary constraints encoded in DNA and associated with the 3D genome. When comparing vertebrate genomes, our results reveal that 3DR which underlies CTCF looping and topologically associating domain organization evolves over time and suggest that ancestral character reconstruction can be used to infer 3DR in ancestral genomes. The R code is available at https://github.com/morphos30/PhyloCTCFLooping. Supplementary data are available at Bioinformatics online.