Potential role of the heat shock protein 90 (hsp90) in buffering mutations to favour cyclical parthenogenesis in the peach potato aphid Myzus persicae (Aphididae, Hemiptera)

Heat-shock proteins 90 (hsp90s) are a class of molecules able to stabilize a network of 'client' proteins that are involved in several processes. Furthermore, recent studies indicated that mutations in the hsp90-encoding gene induce a wide range of phenotypic abnormalities, which have been interpreted as an increased sensitivity of different developmental pathways to hidden/cryptic mutations. In order to verify the role of hsp90 in aphids, we amplified and sequenced the hsp90 gene in 17 lineages of the peach potato aphid Myzus persicae (Sulzer, 1776) looking for the presence of mutations. In particular, we compared lineages with different reproductive modes (obligate vs. cyclical parthenogenesis), propensity to develop winged females and karyotype stability. Differently from the cyclical parthenogenetic lineages that possessed functional hsp90 genes, the seven analysed asexual lineages showed severe mutations (including frameshift and non-sense mutations). In vivo functional assays with the hsp90-inhibitor geldanamycin showed that some lineages with cyclical parthenogenesis may lose their ability to induce sexuales in the absence of active hsp90 revealing the presence of cryptic mutations in their genomes. As a whole, our data suggest that hsp90 could play in aphids a role in buffering hidden/cryptic mutations that disrupt cyclical parthenogenesis.

A key role for sex chromosomes in the regulation of parthenogenesis in the brown alga Ectocarpus

Although evolutionary transitions from sexual to asexual reproduction are frequent in eukaryotes, the genetic bases of these shifts remain largely elusive. Here, we used classic quantitative trait analysis, combined with genomic and transcriptomic information to dissect the genetic basis of asexual, parthenogenetic reproduction in the brown alga Ectocarpus. We found that parthenogenesis is controlled by the sex locus, together with two additional autosomal loci, highlighting the key role of the sex chromosome as a major regulator of asexual reproduction. We identify several negative effects of parthenogenesis on male fitness, and different fitness effects of parthenogenetic capacity depending on the life cycle generation. Although allele frequencies in natural populations are currently unknown, we discuss the possibility that parthenogenesis may be under both sex-specific selection and generation/ploidally-antagonistic selection, and/or that the action of fluctuating selection on this trait may contribute to the maintenance of polymorphisms in populations. Importantly, our data provide the first empirical illustration, to our knowledge, of a trade-off between the haploid and diploid stages of the life cycle, where distinct parthenogenesis alleles have opposing effects on sexual and asexual reproduction and may help maintain genetic variation. These types of fitness trade-offs have profound evolutionary implications in natural populations and may structure life history evolution in organisms with haploid-diploid life cycles.

A Single SNP Turns a Social Honey Bee (Apis mellifera) Worker into a Selfish Parasite

The evolution of altruism in complex insect societies is arguably one of the major transitions in evolution and inclusive fitness theory plausibly explains why this is an evolutionary stable strategy. Yet, workers of the South African Cape honey bee (Apis mellifera capensis) can reverse to selfish behavior by becoming social parasites and parthenogenetically producing female offspring (thelytoky). Using a joint mapping and population genomics approach, in combination with a time-course transcript abundance dynamics analysis, we show that a single nucleotide polymorphism at the mapped thelytoky locus (Th) is associated with the iconic thelytokous phenotype. Th forms a linkage group with the ecdysis-triggering hormone receptor (Ethr) within a nonrecombining region under strong selection in the genome. A balanced detrimental allele system plausibly explains why the trait is specific to A. m. capensis and cannot easily establish itself into genomes of other honey bee subspecies.

Translocation of a parthenogenesis gene candidate to an alternate carrier chromosome in apomictic Brachiaria humidicola

BACKGROUND

The apomictic reproductive mode of Brachiaria (syn. Urochloa) forage species allows breeders to faithfully propagate heterozygous genotypes through seed over multiple generations. In Brachiaria, reproductive mode segregates as single dominant locus, the apospory-specific genomic region (ASGR). The AGSR has been mapped to an area of reduced recombination on Brachiaria decumbens chromosome 5. A primer pair designed within ASGR-BABY BOOM-like (BBML), the candidate gene for the parthenogenesis component of apomixis in Pennisetum squamulatum, was diagnostic for reproductive mode in the closely related species B. ruziziensis, B. brizantha, and B. decumbens. In this study, we used a mapping population of the distantly related commercial species B. humidicola to map the ASGR and test for conservation of ASGR-BBML sequences across Brachiaria species.

RESULTS

Dense genetic maps were constructed for the maternal and paternal genomes of a hexaploid (2n = 6x = 36) B. humidicola F₁ mapping population (n = 102) using genotyping-by-sequencing, simple sequence repeat, amplified fragment length polymorphism, and transcriptome derived single nucleotide polymorphism markers. Comparative genomics with Setaria italica provided confirmation for x = 6 as the base chromosome number of B. humidicola. High resolution molecular karyotyping indicated that the six homologous chromosomes of the sexual female parent paired at random, whereas preferential pairing of subgenomes was observed in the apomictic male parent. Furthermore, evidence for compensated aneuploidy was found in the apomictic parent, with only five homologous linkage groups identified for chromosome 5 and seven homologous linkage groups of chromosome 6. The ASGR mapped to B. humidicola chromosome 1, a region syntenic with chromosomes 1 and 7 of S. italica. The ASGR-BBML specific PCR product cosegregated with the ASGR in the F₁ mapping population, despite its location on a different carrier chromosome than B. decumbens.

CONCLUSIONS

The first dense molecular maps of B. humidicola provide strong support for cytogenetic evidence indicating a base chromosome number of six in this species. Furthermore, these results show conservation of the ASGR across the Paniceae in different chromosomal backgrounds and support postulation of the ASGR-BBML as candidate genes for the parthenogenesis component of apomixis.

Aberrant Stamen Development is Associated with Parthenocarpic Fruit Set Through Up-Regulation of Gibberellin Biosynthesis in Tomato

Parthenocarpy, a process in which fruit set occurs without fertilization, leads to the production of seedless fruit. A number of floral homeotic mutants with abnormal stamen development exhibit parthenocarpic fruit set. Flower development is thought to repress ovary growth before anthesis. However, the mechanism of parthenocarpic fruit development caused by aberrant flower formation is poorly understood. To investigate the molecular mechanism of parthenocarpic fruit development in floral homeotic mutants, we performed functional analysis of Tomato APETALA3 (TAP3) by loss-of-function approaches. Organ-specific promoter was used to induce organ-specific loss of function in stamen and ovary/fruit. We observed increased cell expansion in tap3 mutants and TAP3-RNAi lines during parthenocarpic fruit growth. These were predominantly accompanied by the up-regulation of GA biosynthesis genes, including SlGA20ox1, SlGA20ox2, and SlGA20ox3, as well as reduced expression of the GA-inactivating gene SlGA2ox1 and the auxin signaling gene SlARF7 involved in a crosstalk between GA and auxin. These transcriptional profiles are in agreement with the GA levels in these lines. These results suggest that stamen development negatively regulates fruit set by repressing the GA biosynthesis.

DNA damage and micronuclei in parthenogenetic and bisexual Darevskia rock lizards from the areas with different levels of soil pollution

Natural species are widely used as indicator organisms to estimate of the impact of environmental pollution. Here we present the results of first study of a reliability of parthenogenetic Darevskia аrmeniaca and bisexual Darevskia raddei rock lizards as sentinels for monitoring of environmental genotoxicity. The comet assay and micronucleus test were applied to the lizards sampled in six areas in Armenia and Artsakh with different levels of soil contamination. The results obtained showed a clear relationship between the pollution level of lizards' habitats and the frequency of DNA damage in the comet assay. Low baseline frequency of micronuclei in D. аrmeniaca and D. raddei, however, makes this parameter ineffective for environmental genotoxicity evaluation. The parthenogenetic lizards D. аrmeniaca showed higher sensitivity toward genotoxic pollutions compared with bisexual D. raddei living in the same environment. The correlations between soil content of heavy metals Cr, Cu, Zn, Mo, Pb and DNA damage in D. аrmeniaca and between Cu, As, Mo, Pb and DNA damage in D. raddei were revealed. Overall, the lizards D. raddei and D. аrmeniaca appeared to be sensitive species in detecting soil pollution in natural environment. The application of the comet assay in Darevskia lizard species can be considered as a more appropriate method than a micronucleus test. The use of parthenogenetic lizards D. аrmeniaca as bioindicator will permit to assess the environmental genotoxicity independent of the genetic polymorphism of bisexual species.

The parthenocarpic gene Pat-k is generated by a natural mutation of SlAGL6 affecting fruit development in tomato (Solanum lycopersicum L.)

BACKGROUND

Parthenocarpy is a desired trait in tomato because it can overcome problems with fruit setting under unfavorable environmental conditions. A parthenocarpic tomato cultivar, 'MPK-1', with a parthenocarpic gene, Pat-k, exhibits stable parthenocarpy that produces few seeds. Because 'MPK-1' produces few seeds, seedlings are propagated inefficiently via cuttings. It was reported that Pat-k is located on chromosome 1. However, the gene had not been isolated and the relationship between the parthenocarpy and low seed set in 'MPK-1' remained unclear. In this study, we isolated Pat-k to clarify the relationship between parthenocarpy and low seed set in 'MPK-1'.

RESULTS

Using quantitative trait locus (QTL) analysis for parthenocarpy and seed production, we detected a major QTL for each trait on nearly the same region of the Pat-k locus on chromosome 1. To isolate Pat-k, we performed fine mapping using an F₄ population following the cross between a non-parthenocarpic cultivar, 'Micro-Tom' and 'MPK-1'. The results showed that Pat-k was located in the 529 kb interval between two markers, where 60 genes exist. By using data from a whole genome re-sequencing and genome sequence analysis of 'MPK-1', we could identify that the SlAGAMOUS-LIKE 6 (SlAGL6) gene of 'MPK-1' was mutated by a retrotransposon insertion. The transcript level of SlAGL6 was significantly lower in ovaries of 'MPK-1' than a non-parthenocarpic cultivar. From these results, we could conclude that Pat-k is SlAGL6, and its down-regulation in 'MPK-1' causes parthenocarpy and low seed set. In addition, we observed abnormal micropyles only in plants homozygous for the 'MPK-1' allele at the Pat-k/SlAGL6 locus. This result suggests that Pat-k/SlAGL6 is also related to ovule formation and that the low seed set in 'MPK-1' is likely caused by abnormal ovule formation through down-regulation of Pat-k/SlAGL6.

CONCLUSIONS

Pat-k is identical to SlAGL6, and its down-regulation causes parthenocarpy and low seed set in 'MPK-1'. Moreover, down-regulation of Pat-k/SlAGL6 could cause abnormal ovule formation, leading to a reduction in the number of seeds.

The gene encoding the insulin-like androgenic gland hormone in an all-female parthenogenetic crayfish

Male sexual differentiation in crustaceans is controlled by the androgenic gland (AG), a unique male endocrine organ that, in decapods, is located at the base of the 5th pereiopod. In these animals, the insulin-like androgenic gland hormone (IAG) is the major factor secreted from the AG to induce masculinization and maintain male characteristics. It has, however, recently been proposed that this hormone also plays a role in growth and ovarian development in females. In this study, we tested such a possibility by searching for the IAG gene in the marbled crayfish, a parthenogenetic animal that reproduces asexually to form an all-female genetic clone. Based on the phylogenetic relationship between the marbled crayfish and Procambarus fallax, a gonochoristic species of the same North American Cambaridae family, we searched for the IAG gene in the marbled crayfish and then fully sequenced it. The open reading frame of the gene was found to be completely identical in the two species, and their introns shared over 94% identity. It was also found that, in addition to its expression at the base of the 5th pereiopod and in the testes of male P. fallax crayfish, IAG was expressed in the muscle tissue of P. fallax males and females and even of the parthenogenetic marbled crayfish. These findings provide new insight into possible functions of IAG, in addition to its role as a masculinization-inducing factor, and also constitute the basis for a discussion of the evolutionary relationship between the above two species.

Facultative parthenogenesis validated by DNA analyses in the green anaconda (Eunectes murinus)

In reptiles, the mode of reproduction is typically sexual. However, facultative parthenogenesis occurs in some Squamata, such as Komodo dragon (Varanus komodoensis) and Burmese python (Python bivittatus). Here, we report facultative parthenogenesis in the green anaconda (Eunectes murinus). We found two fully developed female neonates and 17 undeveloped eggs in the oviduct of a female anaconda isolated from other individuals for eight years and two months at Ueno Zoo, Japan. To clarify the zygosity of the neonates, we analyzed 18 microsatellite markers of which 16 were informative. We observed only maternal alleles and no paternal alleles for all 16 markers. To examine the possibility of the long-term sperm storage, we estimated allele frequencies in a putative parental stock by genotyping five unrelated founders. If all founders, including the mother, are originated from a single Mendelian population, then the probability that the neonates were produced by sexual reproduction with an unrelated male via long-term sperm storage was infinitesimally small (2.31E-32 per clutch). We also examined samples from two additional offspring that the mother delivered eight years before her death. We consistently observed paternal alleles in these elder offspring, indicating that the mother had switched from sexual reproduction to asexual reproduction during the eight years of isolation. This is the first case of parthenogenesis in Eunectes to be validated by DNA analysis, and suggests that facultative parthenogenesis is widespread in the Boidae.

Cryptic gene pools in the Hypericum perforatum-H. maculatum complex: diploid persistence versus trapped polyploid melting

BACKGROUND AND AIMS

In Central Europe Hypericum perforatum and Hypericum maculatum show significant hybridization and introgression as a consequence of Pleistocene range fluctuations, and their gene pools are merging on higher ploidy levels. This paper discusses whether polyploid hybrid gene pools are trapped in the ecological climatic niche space of their diploid ancestors, and tests the idea of geographical parthenogenesis.

METHODS

DNA sequence information of nuclear ribosomal DNA and plastid loci, ploidy level estimates and ecological niche modelling are used to characterize the various diploid and polyploid gene pools and unravel spatio-temporal patterns of gene flow among them.

KEY RESULTS

On the diploid level, the three gene pools are clearly distinct between and within species of H. perforatum (two gene pools) and H. maculatum, and their divergence dates back to the first half of the Pleistocene. All polyploids in Central Europe show high levels of past and contemporary gene flow between all three gene pools. The correlation of genetic and geographical distances breaks down if the latter is larger than 250 km, indicating recent and ongoing gene flow. The two species are ecologically differentiated, but in particular hybrids among all three gene pools do not show significant niche differences compared to their parental gene pools, except for some combinations with H. maculatum.

CONCLUSIONS

Inter- and intraspecific gene flow between inter- and intra-species gene pools is limited on the diploid level, and the geographical distribution of the diploids largely reflects Pleistocene evolutionary history. Secondary contact promoted hybridization and introgression on the polyploid level, enabling offspring to escape the diploid gene pools. However, the hybrid polyploids do not show significant niche differences compared to their diploid progenitors. It is concluded that the observed absence of niche divergence has precluded further differentiation and geographical partitioning of new polyploid lineages being effectively separated from the parental lines. The predominantly apomictic reproducing polyploids are trapped in the polyploid gene pool and the ecological climatic niche space of their diploid ancestors.

The origin of multiple clones in the parthenogenetic lizard species Darevskia rostombekowi

The all-female Caucasian rock lizard Darevskia rostombekowi and other unisexual species of this genus reproduce normally via true parthenogenesis. Typically, diploid parthenogenetic reptiles exhibit some amount of clonal diversity. However, allozyme data from D. rostombekowi have suggested that this species consists of a single clone. Herein, we test this hypothesis by evaluating variation at three variable microsatellite loci for 42 specimens of D. rostombekowi from four populations in Armenia. Analyses based on single nucleotide polymorphisms of each locus reveal five genotypes or presumptive clones in this species. All individuals are heterozygous at the loci. The major clone occurs in 24 individuals and involves three populations. Four rare clones involve one or several individuals from one or two populations. Most variation owes to parent-specific single nucleotide polymorphisms, which occur as heterozygotes. This result fails to reject the hypothesis of a single hybridization founder event that resulted in the initial formation of one major clone. The other clones appear to have originated via post-formation microsatellite mutations of the major clone.

Genome analysis of Diploscapter coronatus: insights into molecular peculiarities of a nematode with parthenogenetic reproduction

BACKGROUND

Sexual reproduction involving the fusion of egg and sperm is prevailing among eukaryotes. In contrast, the nematode Diploscapter coronatus, a close relative of the model Caenorhabditis elegans, reproduces parthenogenetically. Neither males nor sperm have been observed and some steps of meiosis are apparently skipped in this species. To uncover the genomic changes associated with the evolution of parthenogenesis in this nematode, we carried out a genome analysis.

RESULTS

We obtained a 170 Mbp draft genome in only 511 scaffolds with a N50 length of 1 Mbp. Nearly 90% of these scaffolds constitute homologous pairs with a 5.7% heterozygosity on average and inversions and translocations, meaning that the 170 Mbp sequences correspond to the diploid genome. Fluorescent staining shows that the D. coronatus genome consists of two chromosomes (2n = 2). In our genome annotation, we found orthologs of 59% of the C. elegans genes. However, a number of genes were missing or very divergent. These include genes involved in sex determination (e.g. xol-1, tra-2) and meiosis (e.g. the kleisins rec-8 and coh-3/4) giving a possible explanation for the absence of males and the second meiotic division. The high degree of heterozygosity allowed us to analyze the expression level of individual alleles. Most of the homologous pairs show very similar expression levels but others exhibit a 2-5-fold difference.

CONCLUSIONS

Our high-quality draft genome of D. coronatus reveals the peculiarities of the genome of parthenogenesis and provides some clues to the genetic basis for parthenogenetic reproduction. This draft genome should be the basis to elucidate fundamental questions related to parthenogenesis such as its origin and mechanisms through comparative analyses with other nematodes. Furthermore, being the closest outgroup to the genus Caenorhabditis, the draft genome will help to disclose many idiosyncrasies of the model C. elegans and its congeners in future studies.

Asexual but Not Clonal: Evolutionary Processes in Automictic Populations

Many parthenogenetically reproducing animals produce offspring not clonally but through different mechanisms collectively referred to as automixis. Here, meiosis proceeds normally but is followed by a fusion of meiotic products that restores diploidy. This mechanism typically leads to a reduction in heterozygosity among the offspring compared to the mother. Following a derivation of the rate at which heterozygosity is lost at one and two loci, depending on the number of crossovers between loci and centromere, a number of models are developed to gain a better understanding of basic evolutionary processes in automictic populations. Analytical results are obtained for the expected neutral genetic variation, effective population size, mutation-selection balance, selection with overdominance, the spread of beneficial mutations, and selection on crossover rates. These results are complemented by numerical investigations elucidating how associative overdominance (two off-phase deleterious mutations at linked loci behaving like an overdominant locus) can in some cases maintain heterozygosity for prolonged times, and how clonal interference affects adaptation in automictic populations. These results suggest that although automictic populations are expected to suffer from the lack of gene shuffling with other individuals, they are nevertheless, in some respects, superior to both clonal and outbreeding sexual populations in the way they respond to beneficial and deleterious mutations. Implications for related genetic systems such as intratetrad mating, clonal reproduction, selfing, as well as different forms of mixed sexual and automictic reproduction are discussed.

Molecular characterization and functional analyses of a diapause hormone receptor-like gene in parthenogenetic Artemia

In arthropods, mature females under certain conditions produce and release encysted gastrula embryos that enter diapause, a state of obligate dormancy. The process is presumably regulated by diapause hormone (DH) and diapause hormone receptor (DHR) that were identified in the silkworm, Bombyx mori and other insects. However, the molecular structure and function of DHR in crustaceans remains unknown. Here, a DHR-like gene from parthenogenetic Artemia (Ar-DHR) was isolated and sequenced. The cDNA sequence consists of 1410bp with a 1260-bp open reading frame encoding a protein consisting of 420 amino acid residues. The results of real-time PCR (qRT-PCR) and Western blot analysis showed that the mRNA and protein of Ar-DHR were mainly expressed at the diapause stage. Furthermore, we found that Ar-DHR was located on the cell membrane of the pre-diapause cyst but in the cytoplasm of the diapause cyst by analysis of immunofluorescence. In vivo knockdown of Ar-DHR by RNA interference (RNAi) and antiserum neutralization consistently inhibited diapause cysts formation. The results indicated that Ar-DHR plays an important role in the induction and maintenance of embryonic diapause in Artemia. Thus, our findings provide an insight into the regulation of diapause formation in Artemia and the function of Ar-DHR.

Rapid breeding of parthenocarpic tomato plants using CRISPR/Cas9

Parthenocarpy in horticultural crop plants is an important trait with agricultural value for various industrial purposes as well as direct eating quality. Here, we demonstrate a breeding strategy to generate parthenocarpic tomato plants using the CRISPR/Cas9 system. We optimized the CRISPR/Cas9 system to introduce somatic mutations effectively into SlIAA9-a key gene controlling parthenocarpy-with mutation rates of up to 100% in the T0 generation. Furthermore, analysis of off-target mutations using deep sequencing indicated that our customized gRNAs induced no additional mutations in the host genome. Regenerated mutants exhibited morphological changes in leaf shape and seedless fruit-a characteristic of parthenocarpic tomato. And the segregated next generation (T1) also showed a severe phenotype associated with the homozygous mutated genome. The system developed here could be applied to produce parthenocarpic tomato in a wide variety of cultivars, as well as other major horticultural crops, using this precise and rapid breeding technique.

The nuclear 28S gene fragment D3 as species marker in oribatid mites (Acari, Oribatida) from German peatlands

To make oribatid mites an applicable tool in monitoring programs it is necessary to find a molecular species marker that allows distinct, rapid and easy species identification. In previous studies, the common barcoding sequence COI showed to be too variable to serve as species marker in oribatid mites. The aim of the present study is to evaluate the potential use of the D3 region of the nuclear 28S rDNA gene for species identification. Therefore, we generated a reference DNA library of 28S D3 to identify specimens of the Oribatida from Germany, with focus on species occurring in peatlands being one of the most endangered habitats in Europe. New DNA sequences were obtained from 325 individuals and 64 species (58 genera, 34 families). By adding 28S D3-sequences from GenBank we altogether analysed 385 sequences from 89 German species, 32 of them restricted to peatlands and further 42 occurring in peatlands occasionally, representing 46 and 33% of the oribatids in German peatlands, respectively. P-distances were measured between species within families as well as for intraspecific divergence. 28S D3 showed low intraspecific genetic p-distances between 0 and 0.5%, interspecific distances within families varied between 0 and 9.7%. Most species pairs within families were further separated by one to four indels in addition to substitutions. Altogether, 93% of all analysed species are clearly delineated by 28S D3. Our study emphasises that 28S D3 rDNA is a useful barcode for the identification of oribatid mite specimens and represents an important step in building-up a comprehensive barcode library to allow metabarcoding analyses of environmental peatland samples for Oribatida in Germany as well as in Central Europe.

Population Structure, Diversity and Reproductive Mode of the Grape Phylloxera (Daktulosphaira vitifoliae) across Its Native Range

Grape Phylloxera, Daktulosphaira vitifoliae, is a gall-forming insect that feeds on the leaves and roots of many Vitis species. The roots of the cultivated V. vinifera cultivars and hybrids are highly susceptible to grape phylloxera feeding damage. The native range of this insect covers most of North America, and it is particularly abundant in the eastern and central United States. Phylloxera was introduced from North America to almost all grape-growing regions across five of the temperate zone continents. It devastated vineyards in each of these regions causing large-scale disruptions to grape growers, wine makers and national economies. In order to understand the population diversity of grape phylloxera in its native range, more than 500 samples from 19 States and 34 samples from the introduced range (northern California, Europe and South America) were genotyped with 32 simple sequence repeat markers. STRUCTURE, a model based clustering method identified five populations within these samples. The five populations were confirmed by a neighbor-joining tree and principal coordinate analysis (PCoA). These populations were distinguished by their Vitis species hosts and their geographic locations. Samples collected from California, Europe and South America traced back to phylloxera sampled in the northeastern United States on V. riparia, with some influence from phylloxera collected along the Atlantic Coast and Central Plains on V. vulpina. Reproductive statistics conclusively confirmed that sexual reproduction is common in the native range and is combined with cyclical parthenogenesis. Native grape phylloxera populations were identified to be under Hardy-Weinberg equilibrium. The identification of admixed samples between many of these populations indicates that shared environments facilitate sexual reproduction between different host associated populations to create new genotypes of phylloxera. This study also found that assortative mating might occur across the sympatric range of the V. vulpina west and V. cinerea populations.

Sex Determination, Sex Chromosomes, and Karyotype Evolution in Insects

Insects harbor a tremendous diversity of sex determining mechanisms both within and between groups. For example, in some orders such as Hymenoptera, all members are haplodiploid, whereas Diptera contain species with homomorphic as well as male and female heterogametic sex chromosome systems or paternal genome elimination. We have established a large database on karyotypes and sex chromosomes in insects, containing information on over 13000 species covering 29 orders of insects. This database constitutes a unique starting point to report phylogenetic patterns on the distribution of sex determination mechanisms, sex chromosomes, and karyotypes among insects and allows us to test general theories on the evolutionary dynamics of karyotypes, sex chromosomes, and sex determination systems in a comparative framework. Phylogenetic analysis reveals that male heterogamety is the ancestral mode of sex determination in insects, and transitions to female heterogamety are extremely rare. Many insect orders harbor species with complex sex chromosomes, and gains and losses of the sex-limited chromosome are frequent in some groups. Haplodiploidy originated several times within insects, and parthenogenesis is rare but evolves frequently. Providing a single source to electronically access data previously distributed among more than 500 articles and books will not only accelerate analyses of the assembled data, but also provide a unique resource to guide research on which taxa are likely to be informative to address specific questions, for example, for genome sequencing projects or large-scale comparative studies.

The Genetic Basis of Natural Variation in Drosophila (Diptera: Drosophilidae) Virgin Egg Retention

Drosophila melanogaster is able to thrive in harsh northern climates through adaptations in life-history traits and physiological mechanisms that allow for survival through the winter. We examined the genetic basis of natural variation in one such trait, female virgin egg retention, which was previously shown to vary clinally and seasonally. To further our understanding of the genetic basis and evolution of virgin egg retention, we performed a genome-wide association study (GWAS) using the previously sequenced Drosophila Genetic Reference Panel (DGRP) mapping population. We found 29 single nucleotide polymorphisms (SNPs) associated with virgin egg retention and assayed 6 available mutant lines, each harboring a mutation in a candidate gene, for effects on egg retention time. We found that four out of the six mutant lines had defects in egg retention time as compared with the respective controls: mun, T48, Mes-4, and Klp67A Surprisingly, none of these genes has a recognized role in ovulation control, but three of the four genes have known effects on fertility or have high expression in the ovaries. We also found that the SNP set associated with egg retention time was enriched for clinal SNPs. The majority of clinal SNPs had alleles associated with longer egg retention present at higher frequencies in higher latitudes. Our results support previous studies that show higher frequency of long retention times at higher latitude, providing evidence for the adaptive value of virgin egg-retention.

QTL mapping of parthenocarpic fruit set in North American processing cucumber

Through a novel phenotyping method, four QTLs were consistently associated with increased parthenocarpic fruit set in North American processing cucumber that accounted for over 75 % of observed phenotypic variation. Parthenocarpy is a desirable trait with potential for increasing yield and quality in processing cucumber production. Although many successful parthenocarpic fresh market cucumber varieties have been developed, the genetic and molecular mechanisms behind parthenocarpic expression in cucumber remain largely unknown. Since parthenocarpy is an important yield component, it is difficult to separate the true parthenocarpic character from other yield related traits. In the present study, we developed a novel phenotypic approach for parthenocarpic fruit set focusing on early fruit development. Two hundred and five F₃ families derived from a cross between the highly parthenocarpic line 2A and low parthenocarpic line Gy8 were phenotypically evaluated in three greenhouse experiments. Seven QTLs associated with parthenocarpic fruit set were detected. Among them, one each on chromosomes 5 and 7 (parth5.1 and parth7.1) and two on chromosome 6 (parth6.1 and parth6.2) were consistently identified in all experiments, but their relative contribution to the total phenotypic variation was dependent on plant growth stages. While each of the four QTLs had almost equal contribution to the expression of the trait at commercial harvest stage, parth7.1 played an important role in early parthenocarpic fruit set. The results suggested that parthenocarpic fruit set can be accurately evaluated with as few as 20 nodes of growth. The QTLs identified in this study for parthenocarpic fruit set are a valuable resource for cucumber breeders interested in developing parthenocarpic cultivars and to researchers interested in the genetic and molecular mechanisms of parthenocarpic fruit set.

Maintenance of Xist Imprinting Depends on Chromatin Condensation State and Rnf12 Dosage in Mice

In female mammals, activation of Xist (X-inactive specific transcript) is essential for establishment of X chromosome inactivation. During early embryonic development in mice, paternal Xist is preferentially expressed whereas maternal Xist (Xm-Xist) is silenced. Unlike autosomal imprinted genes, Xist imprinting for Xm-Xist silencing was erased in cloned or parthenogenetic but not fertilized embryos. However, the molecular mechanism underlying the variable nature of Xm-Xist imprinting is poorly understood. Here, we revealed that Xm-Xist silencing depends on chromatin condensation states at the Xist/Tsix genomic region and on Rnf12 expression levels. In early preimplantation, chromatin decondensation via H3K9me3 loss and histone acetylation gain caused Xm-Xist derepression irrespective of embryo type. Although the presence of the paternal genome during pronuclear formation impeded Xm-Xist derepression, Xm-Xist was robustly derepressed when the maternal genome was decondensed before fertilization. Once Xm-Xist was derepressed by chromatin alterations, the derepression was stably maintained and rescued XmXp^Δ lethality, indicating that loss of Xm-Xist imprinting was irreversible. In late preimplantation, Oct4 served as a chromatin opener to create transcriptional permissive states at Xm-Xist/Tsix genomic loci. In parthenogenetic embryos, Rnf12 overdose caused Xm-Xist derepression via Xm-Tsix repression; physiological Rnf12 levels were essential for Xm-Xist silencing maintenance in fertilized embryos. Thus, chromatin condensation and fine-tuning of Rnf12 dosage were crucial for Xist imprint maintenance by silencing Xm-Xist.

X-inactive specific transcript (Xist) is essential a large non-coding RNA for establishment of X chromosome inactivation in female mammals. The aberrant X chromosome inactivation critically affects cellular viability. Therefore, spatiotemporal regulation of Xist expression is required for proper development. In mice, Xist expression is imprinted in early embryonic development and maternal Xist is never expressed during preimplantation phases irrespective of the presence of Xist activator, maternal Rnf12. Generally, parental origin-specific expression pattern of autosomal imprinted genes is maintained in various types of embryos. However, Xist imprinting for transcriptional silencing of maternal Xist was erased in cloned or parthenogenetic but not fertilized embryos. Here, we dissect the molecular mechanism underlying the variable nature of Xist imprinting. We show that in fertilized embryos, chromatin condensation states are essential maternal Xist repression in early preimplantation phases, whereas at late preimplantation stages, pluripotency factor Oct4 serves as a chromatin opener and the maintenance of Xist silencing depends on Rnf12 expression dosage. Although the Oct4 mediated chromatin decondensation also occurs in parthenogetic embryos, Rnf12 overdose causes maternal Xist derepression at late preimplantation phases. Thus these findings reveal that the chromatin regulation by pluripotency factor and Xist activator dose define Xist imprinting state.

De Novo Transcriptome Assembly and Sex-Biased Gene Expression in the Cyclical Parthenogenetic Daphnia galeata

Daphnia species have become models for ecological genomics and exhibit interesting features, such as high phenotypic plasticity and a densely packed genome with many lineage-specific genes. They are also cyclic parthenogenetic, with alternating asexual and sexual cycles and environmental sex determination. Here, we present a de novo transcriptome assembly of over 32,000 D. galeata genes and use it to investigate gene expression in females and spontaneously produced males of two clonal lines derived from lakes in Germany and the Czech Republic. We find that only a low percentage (18%) of genes shows sex-biased expression and that there are many more female-biased gene (FBG) than male-biased gene (MBG). Furthermore, FBGs tend to be more conserved between species than MBGs in both sequence and expression. These patterns may be a consequence of cyclic parthenogenesis leading to a relaxation of purifying selection on MBGs. The two clonal lines show considerable differences in both number and identity of sex-biased genes, suggesting that they may have reproductive strategies differing in their investment in sexual reproduction. Orthologs of key genes in the sex determination and juvenile hormone pathways, which are thought to be important for the transition from asexual to sexual reproduction, are present in D. galeata and highly conserved among Daphnia species.

Promoter Architecture and Sex-Specific Gene Expression in Daphnia pulex

Large-scale transcription start site (TSS) profiling produces a high-resolution, quantitative picture of transcription initiation and core promoter locations within a genome. However, application of TSS profiling to date has largely been restricted to a small set of prominent model systems. We sought to characterize the cis-regulatory landscape of the water flea Daphnia pulex, an emerging model arthropod that reproduces both asexually (via parthenogenesis) and sexually (via meiosis). We performed Cap Analysis of Gene Expression (CAGE) with RNA isolated from D. pulex within three developmental states: sexual females, asexual females, and males. Identified TSSs were utilized to generate a "Daphnia Promoter Atlas," i.e., a catalog of active promoters across the surveyed states. Analysis of the distribution of promoters revealed evidence for widespread alternative promoter usage in D. pulex, in addition to a prominent fraction of compactly-arranged promoters in divergent orientations. We carried out de novo motif discovery using CAGE-defined TSSs and identified eight candidate core promoter motifs; this collection includes canonical promoter elements (e.g., TATA and Initiator) in addition to others lacking obvious orthologs. A comparison of promoter activities found evidence for considerable state-specific differential gene expression between states. Our work represents the first global definition of transcription initiation and promoter architecture in crustaceans. The Daphnia Promoter Atlas presented here provides a valuable resource for comparative study of cis-regulatory regions in metazoans, as well as for investigations into the circuitries that underpin meiosis and parthenogenesis.

Combined hybridization and mitochondrial capture shape complex phylogeographic patterns in hybridogenetic Cataglyphis desert ants

Some species of Cataglyphis desert ants have evolved a hybridogenetic mode of reproduction at the social scale. In hybridogenetic populations, two distinct genetic lineages coexist. Non-reproductive offspring (workers) are hybrids of the two lineages, whereas sexual offspring (males and new queens) are produced by parthenogenesis and belong to the mother queen lineage. How this unusual reproductive system affects phylogeographic patterns and speciation processes remains completely unknown to date. Using one mitochondrial and four nuclear genes, we examined the phylogenetic relationships between three species of Cataglyphis (C. hispanica, C. humeya and C. velox) where complex DNA inheritance through social hybridogenesis may challenge phylogenetic inference. Our results bring two important insights. First, our data confirm a hybridogenetic mode of reproduction across the whole distribution range of the species C. hispanica. In contrast, they do not provide support for hybridogenesis in the populations sampled of C. humeya and C. velox. This suggests that these populations are not hybridogenetic, or that hybridogenesis is too recent to result in reciprocally monophyletic lineages on nuclear genes. Second, due to mitochondrial introgression between lineages (Darras and Aron, 2015), the faster-evolving COI marker is not lineage specific, hence, unsuitable to further investigate the segregation of lineages in the species studied. Different mitochondrial haplotypes occur in each locality sampled, resulting in strongly structured populations. This micro-allopatric structure leads to over-splitting species delimitation on mitochondrial gene, as every locality could potentially be considered a putative species; haploweb analyses of nuclear markers, however, yield species delimitations that are consistent with morphology. Overall, this study highlights how social hybridogenesis varies across species and shapes complex phylogeographic patterns.

Sexual versus Asexual Reproduction: Distinct Outcomes in Relative Abundance of Parthenogenetic Mealybugs following Recent Colonization

Asexual reproduction, including parthenogenesis in which embryos develop within a female without fertilization, is assumed to confer advantages over sexual reproduction, which includes a “cost of males.” Sexual reproduction largely predominates in animals, however, indicating that this cost is outweighed by the genetic and/or ecological benefits of sexuality, including the acquisition of advantageous mutations occurring in different individuals and the elimination of deleterious mutations. But the evolution of sexual reproduction remains unclear, because we have limited examples that demonstrate the relative success of sexual lineages in the face of competition from asexual lineages in the same environment. Here we investigated a sympatric occurrence of sexual and asexual reproduction in the pineapple mealybug, Dysmicoccus brevipes. This pest invaded southwestern Japan, including Okinawa and Ishigaki Islands, in the 1930s in association with imported pineapple plants. Our recent censuses demonstrated that on Okinawa sexually reproducing individuals can coexist with and even dominate asexual individuals in the presence of habitat and resource competition, which is considered to be severe for this nearly immobile insect. Molecular phylogeny based on partial DNA sequences in the mitochondrial and nuclear genomes, as well as the endosymbiotic bacterial genome, revealed that the asexual lineage diverged from a common sexual ancestor in the relatively recent past. In contrast, only the asexual lineage exhibiting obligate apomictic thelytoky was discovered on Ishigaki. Co-existence of the two lineages cannot be explained by the results of laboratory experiments, which showed that the intrinsic rate of increase in the sexual lineage was not obviously superior to that of the asexual lineage. Differences in biotic and/or abiotic selective forces operating on the two islands might be the cause of this discrepancy. This biological system offers a unique opportunity to assess the relative success of sexual versus asexual lineages with an unusual morphology and life cycle.