Interrogating 1000 insect genomes for NUMTs: A risk assessment for estimates of species richness.

The nuclear genomes of most animal species include NUMTs, segments of the mitogenome incorporated into their chromosomes. Although NUMT counts are known to vary greatly among species, there has been no comprehensive study of their frequency/attributes in the most diverse group of terrestrial organisms, insects. This study examines NUMTs derived from a 658 bp 5' segment of the cytochrome c oxidase I (COI) gene, the barcode region for the animal kingdom. This assessment is important because unrecognized NUMTs can elevate estimates of species richness obtained through DNA barcoding and derived approaches (eDNA, metabarcoding). This investigation detected nearly 10,000 COI NUMTs ≥ 100 bp in the genomes of 1,002 insect species (range = 0-443). Variation in nuclear genome size explained 56% of the mitogenome-wide variation in NUMT counts. Although insect orders with the largest genome sizes possessed the highest NUMT counts, there was considerable variation among their component lineages. Two thirds of COI NUMTs possessed an IPSC (indel and/or premature stop codon) allowing their recognition and exclusion from downstream analyses. The remainder can elevate species richness as they showed 10.1% mean divergence from their mitochondrial homologue. The extent of exposure to "ghost species" is strongly impacted by the target amplicon's length. NUMTs can raise apparent species richness by up to 22% when a 658 bp COI amplicon is examined versus a doubling of apparent richness when 150 bp amplicons are targeted. Given these impacts, metabarcoding and eDNA studies should target the longest possible amplicons while also avoiding use of 12S/16S rDNA as they triple NUMT exposure because IPSC screens cannot be employed.

Geometrical Properties of the Nucleus and Chromosome Intermingling Are Possible Major Parameters of Chromosome Aberration Formation.

Ionizing radiation causes chromosome aberrations, which are possible biomarkers to assess space radiation cancer risks. Using the Monte Carlo codes Relativistic Ion Tracks (RITRACKS) and Radiation-Induced Tracks, Chromosome Aberrations, Repair and Damage (RITCARD), we investigated how geometrical properties of the cell nucleus, irradiated with ion beams of linear energy transfer (LET) ranging from 0.22 keV/µm to 195 keV/µm, influence the yield of simple and complex exchanges. We focused on the effect of (1) nuclear volume by considering spherical nuclei of varying radii; (2) nuclear shape by considering ellipsoidal nuclei of varying thicknesses; (3) beam orientation; and (4) chromosome intermingling by constraining or not constraining chromosomes in non-overlapping domains. In general, small nuclear volumes yield a higher number of complex exchanges, as compared to larger nuclear volumes, and a higher number of simple exchanges for LET < 40 keV/µm. Nuclear flattening reduces complex exchanges for high-LET beams when irradiated along the flattened axis. The beam orientation also affects yields for ellipsoidal nuclei. Reducing chromosome intermingling decreases both simple and complex exchanges. Our results suggest that the beam orientation, the geometry of the cell nucleus, and the organization of the chromosomes within are important parameters for the formation of aberrations that must be considered to model and translate in vitro results to in vivo risks.

A combinatorial indexing strategy for low-cost epigenomic profiling of plant single cells.

Understanding how cis-regulatory elements facilitate gene expression is a key question in biology. Recent advances in single-cell genomics have led to the discovery of cell-specific chromatin landscapes that underlie transcription programs in animal models. However, the high equipment and reagent costs of commercial systems limit their applications for many laboratories. In this study, we developed a combinatorial index and dual PCR barcode strategy to profile the Arabidopsis thaliana root single-cell epigenome without any specialized equipment. We generated chromatin accessibility profiles for 13 576 root nuclei with an average of 12 784 unique Tn5 integrations per cell. Integration of the single-cell assay for transposase-accessible chromatin sequencing and RNA sequencing data sets enabled the identification of 24 cell clusters with unique transcription, chromatin, and cis-regulatory signatures. Comparison with single-cell data generated using the commercial microfluidic platform from 10X Genomics revealed that this low-cost combinatorial index method is capable of unbiased identification of cell-type-specific chromatin accessibility. We anticipate that, by removing cost, instrumentation, and other technical obstacles, this method will be a valuable tool for routine investigation of single-cell epigenomes and provide new insights into plant growth and development and plant interactions with the environment.

Al-Tolerant Barley Mutant hvatr.g Shows the ATR-Regulated DNA Damage Response to Maleic Acid Hydrazide.

ATR, a DNA damage signaling kinase, is required for cell cycle checkpoint regulation and detecting DNA damage caused by genotoxic factors including Al3+ ions. We analyzed the function of the HvATR gene in response to chemical clastogen-maleic acid hydrazide (MH). For this purpose, the Al-tolerant barley TILLING mutant hvatr.g was used. We described the effects of MH on the nuclear genome of hvatr.g mutant and its WT parent cv. "Sebastian", showing that the genotoxic effect measured by TUNEL test and frequency of cells with micronuclei was much stronger in hvatr.g than in WT. MH caused a significant decrease in the mitotic activity of root cells in both genotypes, however this effect was significantly stronger in "Sebastian". The impact of MH on the roots cell cycle, analyzed using flow cytometry, showed no differences between the mutant and WT.

A parameter sensitivity study for simulating DNA damage after proton irradiation using TOPAS-nBio.

Monte Carlo (MC) track structure simulation tools are commonly used for predicting radiation induced DNA damage by modeling the physical and chemical reactions at the nanometer scale. However, the outcome of these MC simulations is particularly sensitive to the adopted parameters which vary significantly across studies. In this study, a previously developed full model of nuclear DNA was used to describe the DNA geometry. The TOPAS-nBio MC toolkit was used to investigate the impact of physics and chemistry models as well as three key parameters (the energy threshold for direct damage, the chemical stage time length, and the probability of damage between hydroxyl radical reactions with DNA) on the induction of DNA damage. Our results show that the difference in physics and chemistry models alone can cause differences up to 34% and 16% in the DNA double strand break (DSB) yield, respectively. Additionally, changing the direct damage threshold, chemical stage length, and hydroxyl damage probability can cause differences of up to 28%, 51%, and 71% in predicted DSB yields, respectively, for the configurations in this study.

Assessment of Nuclear and Mitochondrial DNA, Expression of Mitochondria-Related Genes in Different Brain Regions in Rats after Whole-Body X-ray Irradiation.

Studies of molecular changes occurred in various brain regions after whole-body irradiation showed a significant increase in terms of the importance in gaining insight into how to slow down or prevent the development of long-term side effects such as carcinogenesis, cognitive impairment and other pathologies. We have analyzed nDNA damage and repair, changes in mitochondrial DNA (mtDNA) copy number and in the level of mtDNA heteroplasmy, and also examined changes in the expression of genes involved in the regulation of mitochondrial biogenesis and dynamics in three areas of the rat brain (hippocampus, cortex and cerebellum) after whole-body X-ray irradiation. Long amplicon quantitative polymerase chain reaction (LA-QPCR) was used to detect nDNA and mtDNA damage. The level of mtDNA heteroplasmy was estimated using Surveyor nuclease technology. The mtDNA copy numbers and expression levels of a number of genes were determined by real-time PCR. The results showed that the repair of nDNA damage in the rat brain regions occurs slowly within 24 h; in the hippocampus, this process runs much slower. The number of mtDNA copies in three regions of the rat brain increases with a simultaneous increase in mtDNA heteroplasmy. However, in the hippocampus, the copy number of mutant mtDNAs increases significantly by the time point of 24 h after radiation exposure. Our analysis shows that in the brain regions of irradiated rats, there is a decrease in the expression of genes (ND2, CytB, ATP5O) involved in ATP synthesis, although by the same time point after irradiation, an increase in transcripts of genes regulating mitochondrial biogenesis is observed. On the other hand, analysis of genes that control the dynamics of mitochondria (Mfn1, Fis1) revealed that sharp decrease in gene expression level occurred, only in the hippocampus. Consequently, the structural and functional characteristics of the hippocampus of rats exposed to whole-body radiation can be different, most significantly from those of the other brain regions.

A new open-source GPU-based microscopic Monte Carlo simulation tool for the calculations of DNA damages caused by ionizing radiation --- Part I: Core algorithm and validation.

PURPOSE: Monte Carlo (MC) simulation of radiation interactions with water medium at physical, physicochemical, and chemical stages, as well as the computation of biologically relevant quantities such as DNA damages, are of critical importance for the understanding of microscopic basis of radiation effects. Due to the large problem size and many-body simulation problem in the chemical stage, existing CPU-based computational packages encounter the problem of low computational efficiency. This paper reports our development on a GPU-based microscopic Monte Carlo simulation tool gMicroMC using advanced GPU-acceleration techniques. METHODS: gMicroMC simulated electron transport in the physical stage using an interaction-by-interaction scheme to calculate the initial events generating radicals in water. After the physicochemical stage, initial positions of all radicals were determined. Simulation of radicals' diffusion and reactions in the chemical stage was achieved using a step-by-step model using GPU-accelerated parallelization together with a GPU-enabled box-sorting algorithm to reduce the computations of searching for interaction pairs and therefore improve efficiency. A multi-scale DNA model of the whole lymphocyte cell nucleus containing ~6.2 Gbp of DNA was built. RESULTS: Accuracy of physical stage simulation was demonstrated by computing stopping power and track length. The results agreed with published data and the data produced by GEANT4-DNA (version 10.3.3) simulations with 10 -20% difference in most cases. Difference of yield values of major radiolytic species from GEANT4-DNA results was within 10%. We computed DNA damages caused by monoenergetic 662 keV photons, approximately representing 137 Cs decay. Single-strand break (SSB) and double-strand break (DSB) yields were 196 ± 8 SSB/Gy/Gbp and 7.3 ± 0.7 DSB/Gy/Gbp, respectively, which agreed with the result of 188 SSB/Gy/Gbp and 8.4 DSB/Gy/Gbp computed by Hsiao et al. Compared to computation using a single CPU, gMicroMC achieved a speedup factor of ~540x using an NVidia TITAN Xp GPU card. CONCLUSIONS: The achieved accuracy and efficiency demonstrated that gMicroMC can facilitate research on microscopic radiation transport simulation and DNA damage calculation. gMicroMC is an open-source package available to the research community.

European minnows through time: museum collections aid genetic assessment of species introductions in freshwater fishes (Cyprinidae: Phoxinus species complex).

Massive fish introductions have taken place throughout much of the world, mostly over the last 70 years, and present a major threat to the genetic diversity of native fishes. Introductions have been reported for European Phoxinus, a ubiquitous small cyprinid that populates a wide variety of habitats. Species delineation in European Phoxinus has proven difficult with one reason being ranges of distribution that often traverse drainage boundaries. The present study combines recent samples with museum samples to better understand the current distribution of Phoxinus species and their distributions prior to the massive introductions of fishes in Europe, and to evaluate the use of museum specimens for species distribution studies. For these purposes, genetic lineages from sites collected prior to 1900 (n = 14), and between 1900 and 1950 (n = 8), were analysed using two mitochondrial and nuclear markers. Although possible fish introductions were detected, our results show that the distribution of genetic lineages of museum samples is comparable to that of the extant lineages of European Phoxinus present in those areas. These observations suggest that in the studied ranges the distribution of Phoxinus lineages has been driven by natural processes.

BAT Hi-C maps global chromatin interactions in an efficient and economical way.

Chromosome Conformation Capture (3C)-based technologies, such as Hi-C, have represented a significant breakthrough in investigating the structure and function of higher-order genome architecture. However, the mapping of global chromatin interactions remains challenging across many biological conditions due to high background noise and financial constraints, especially for small laboratories. Here, we describe the Bridge linker-Alul-Tn5 Hi-C (BAT Hi-C) method, which is a simple and efficient method for delineating chromatin conformational features of mouse embryonic stem (mES) cells and uncover DNA loops. This protocol combines Alul fragmentation and biotinylated linker-mediated proximity ligation to obtain kilobase (kb) resolution with a marked increase in the amount of unique read pairs. The protocol also includes chromatin isolation to reduce background noise and Tn5 tagmentation to cut down on preparation time. Importantly, with only one-third sequencing depth, our method revealed the same spectrum of chromatin contacts as in situ Hi-C. BAT Hi-C is an economical (i.e., approximately $40 for library preparation) and straightforward (total hands-on time of 3 days) tool that is ideal for the in-depth analysis of long-range chromatin looping events in a genome-wide fashion.

Mitochondrial and nuclear genetic analyses of the tropical black-lip rock oyster (Saccostrea echinata) reveals population subdivision and informs sustainable aquaculture development.

BACKGROUND: The black-lip rock oyster (Saccostrea echinata) has considerable potential for aquaculture throughout the tropics. Previous attempts to farm S. echinata failed due to an insufficient supply of wild spat; however, the prospect of hatchery-based aquaculture has stimulated renewed interest, and small-scale farming is underway across northern Australia and in New Caledonia. The absence of knowledge surrounding the population genetic structure of this species has raised concerns about the genetic impacts of this emerging aquaculture industry. This study is the first to examine population genetics of S. echinata and employs both mitochondrial cytochrome c oxidase subunit I gene (COI) and single nucleotide polymorphism (SNP) markers. RESULTS: The mitochondrial COI data set included 273 sequences of 594 base pair length, which comprised 74 haplotypes. The SNP data set included 27,887 filtered SNPs for 272 oysters and of these 31 SNPs were identified as candidate adaptive loci. Data from the mitochondrial COI analyses, supports a broad tropical Indo-Pacific distribution of S. echinata, and showed high haplotype and nucleotide diversities (0.887-1.000 and 0.005-0.008, respectively). Mitochondrial COI analyses also revealed a 'star-like' haplotype network, and significant and negative neutrality tests (Tajima's D = - 2.030, Fu's Fs = - 25.638, P < 0.001) support a recent population expansion after a bottleneck. The SNP analyses showed significant levels of population subdivision and four genetic clusters were identified: (1) the Noumea (New Caledonia) sample location; (2) the Bowen (north Queensland, Australia) sample location, and remaining sample locations in the Northern Territory, Australia (n = 8) were differentiated into two genetic clusters. These occurred at either side of the Wessel Islands and were termed (3) 'west' and (4) 'east' clusters, and two migrant individuals were detected between them. The SNP data showed a significant positive correlation between genetic and geographic distance (Mantel test, P < 0.001, R2 = 0.798) and supported isolation by distance. Three candidate adaptive SNPs were identified as occurring within known genes and gene ontology was well described for the sex peptide receptor gene. CONCLUSIONS: Data supports the existence of genetically distinct populations of S. echinata, suggesting that management of wild and farmed stocks should be based upon multiple management units. This research has made information on population genetic structure and connectivity available for a new aquaculture species.

Identification and assessment of variable single-copy orthologous (SCO) nuclear loci for low-level phylogenomics: a case study in the genus Rosa (Rosaceae).

BACKGROUND: With an ever-growing number of published genomes, many low levels of the Tree of Life now contain several species with enough molecular data to perform shallow-scale phylogenomic studies. Moving away from using just a few universal phylogenetic markers, we can now target thousands of other loci to decipher taxa relationships. Making the best possible selection of informative sequences regarding the taxa studied has emerged as a new issue. Here, we developed a general procedure to mine genomic data, looking for orthologous single-copy loci capable of deciphering phylogenetic relationships below the generic rank. To develop our strategy, we chose the genus Rosa, a rapid-evolving lineage of the Rosaceae family in which several species genomes have recently been sequenced. We also compared our loci to conventional plastid markers, commonly used for phylogenetic inference in this genus. RESULTS: We generated 1856 sequence tags in putative single-copy orthologous nuclear loci. Associated in silico primer pairs can potentially amplify fragments able to resolve a wide range of speciation events within the genus Rosa. Analysis of parsimony-informative site content showed the value of non-coding genomic regions to obtain variable sequences despite the fact that they may be more difficult to target in less related species. Dozens of nuclear loci outperform the conventional plastid phylogenetic markers in terms of phylogenetic informativeness, for both recent and ancient evolutionary divergences. However, conflicting phylogenetic signals were found between nuclear gene tree topologies and the species-tree topology, shedding light on the many patterns of hybridization and/or incomplete lineage sorting that occur in the genus Rosa. CONCLUSIONS: With recently published genome sequence data, we developed a set of single-copy orthologous nuclear loci to resolve species-level phylogenomics in the genus Rosa. This genome-wide scale dataset contains hundreds of highly variable loci which phylogenetic interest was assessed in terms of phylogenetic informativeness and topological conflict. Our target identification procedure can easily be reproduced to identify new highly informative loci for other taxonomic groups and ranks.

A cell-by-cell Monte Carlo simulation for assessing radiation-induced DNA double strand breaks.

This paper presents a cell-by-cell Monte Carlo simulation study that combines charged particle track structure data with an interphase cell nucleus model to quantify DNA double strand breaks (DSBs), spatial distribution of DSBs in a cell nucleus, and resulting potentially lethal or mutagenic events (PLMEs) between DSBs in close proximity. Cell nucleus is simulated according to the chromosome territory-interchromatin compartment (CT-IC) model in that chromatin content is unevenly distributed in chromatin domains (CDs) and IC with a chromatin compaction ratio of 22:1. A particle track structure coordinate (PTSC) library was first generated for each particle type, energy, and dose based on a large number of particle track data obtained by running the Monte Carlo track structure code Geant4-DNA. To assess the DNA DSBs of a cell for a specific particle type, energy, and dose, the corresponding PTSC was selected and "map overlaid" onto 960 unique cell nucleus data sets containing chromatin fiber (CF) locations. Clustering algorithm DBSCAN was next used to identify the clustered energy deposition events occurring inside the CF. These events were then converted to DNA DSBs using a probabilistic approach. The locations of the DSBs thus obtained were, in turn, used to calculate PLMEs within the cell nucleus that can result from DSB proximity and complexity. The results obtained from this simulation study are correctly correlated to the experimental data of DSB yield and the RBE-LET relationships for various types of charged particles and of various energies. The results show agreement with other published radiobiological models.

Evaluation of early radiation DNA damage in a fractal cell nucleus model using Geant4-DNA.

The advancement of multidisciplinary research fields dealing with ionising radiation induced biological damage - radiobiology, radiation physics, radiation protection and, in particular, medical physics - requires a clear mechanistic understanding of how cellular damage is induced by ionising radiation. Monte Carlo (MC) simulations provide a promising approach for the mechanistic simulation of radiation transport and radiation chemistry, towards the in silico simulation of early biological damage. We have recently developed a fully integrated MC simulation that calculates early single strand breaks (SSBs) and double strand breaks (DSBs) in a fractal chromatin based human cell nucleus model. The results of this simulation are almost equivalent to past MC simulations when considering direct/indirect strand break fraction, DSB yields and fragment distribution. The simulation results agree with experimental data on DSB yields within 13.6% on average and fragment distributions agree within an average of 34.8%.

Killer whale genomes reveal a complex history of recurrent admixture and vicariance.

Reconstruction of the demographic and evolutionary history of populations assuming a consensus tree-like relationship can mask more complex scenarios, which are prevalent in nature. An emerging genomic toolset, which has been most comprehensively harnessed in the reconstruction of human evolutionary history, enables molecular ecologists to elucidate complex population histories. Killer whales have limited extrinsic barriers to dispersal and have radiated globally, and are therefore a good candidate model for the application of such tools. Here, we analyse a global data set of killer whale genomes in a rare attempt to elucidate global population structure in a nonhuman species. We identify a pattern of genetic homogenisation at lower latitudes and the greatest differentiation at high latitudes, even between currently sympatric lineages. The processes underlying the major axis of structure include high drift at the edge of species' range, likely associated with founder effects and allelic surfing during postglacial range expansion. Divergence between Antarctic and non-Antarctic lineages is further driven by ancestry segments with up to four-fold older coalescence time than the genome-wide average; relicts of a previous vicariance during an earlier glacial cycle. Our study further underpins that episodic gene flow is ubiquitous in natural populations, and can occur across great distances and after substantial periods of isolation between populations. Thus, understanding the evolutionary history of a species requires comprehensive geographic sampling and genome-wide data to sample the variation in ancestry within individuals.

tRNA Fragments Show Intertwining with mRNAs of Specific Repeat Content and Have Links to Disparities.

tRNA-derived fragments (tRF) are a class of potent regulatory RNAs. We mined the datasets from The Cancer Genome Atlas (TCGA) representing 32 cancer types with a deterministic and exhaustive pipeline for tRNA fragments. We found that mitochondrial tRNAs contribute disproportionally more tRFs than nuclear tRNAs. Through integrative analyses, we uncovered a multitude of statistically significant and context-dependent associations between the identified tRFs and mRNAs. In many of the 32 cancer types, these associations involve mRNAs from developmental processes, receptor tyrosine kinase signaling, the proteasome, and metabolic pathways that include glycolysis, oxidative phosphorylation, and ATP synthesis. Even though the pathways are common to multiple cancers, the association of specific mRNAs with tRFs depends on and differs from cancer to cancer. The associations between tRFs and mRNAs extend to genomic properties as well; specifically, tRFs are positively correlated with shorter genes that have a higher density in repeats, such as ALUs, MIRs, and ERVLs. Conversely, tRFs are negatively correlated with longer genes that have a lower repeat density, suggesting a possible dichotomy between cell proliferation and differentiation. Analyses of bladder, lung, and kidney cancer data indicate that the tRF-mRNA wiring can also depend on a patient's sex. Sex-dependent associations involve cyclin-dependent kinases in bladder cancer, the MAPK signaling pathway in lung cancer, and purine metabolism in kidney cancer. Taken together, these findings suggest diverse and wide-ranging roles for tRFs and highlight the extensive interconnections of tRFs with key cellular processes and human genomic architecture. SIGNIFICANCE: Across 32 TCGA cancer contexts, nuclear and mitochondrial tRNA fragments exhibit associations with mRNAs that belong to concrete pathways, encode proteins with particular destinations, have a biased repeat content, and are sex dependent.

Authentication and Assessment of Contamination in Ancient DNA.

Contamination from both present-day humans and postmortem microbial sources is a common challenge in ancient DNA studies. Here we present a suite of tools to assist in the assessment of contamination in ancient DNA data sets. These tools perform standard tests of authenticity of ancient DNA data including detecting the presence of postmortem damage signatures in sequence alignments and quantifying the amount of present-day human contamination.

Influence of chromatin compaction on simulated early radiation-induced DNA damage using Geant4-DNA.

PURPOSE: In this work, we present simulated double-strand breaks (DSBs) obtained for two human cell nucleus geometries. The first cell nucleus represents fibroblasts, filled with DNA molecules in different compaction forms: heterochromatin or euchromatin only. The second one represents an endothelial cell nucleus, either filled with heterochromatin only or with a uniform distribution of 48% of heterochromatin and 52% of euchromatin, obtained from measurements carried out at IRSN. Protons and alpha particles of different energies were used as projectiles. Each cell nucleus model includes a multi-scale description of the DNA target from the molecular level to the whole human genome representation. METHODS: The cell nucleus models were generated using an extended version of the DnaFabric software in which a new model of euchromatin was implemented in addition to the existing model of heterochromatin. Thus, each nucleus model contains the complete human genome (a total of 6 Gbp) in the G0/G1 phase of the cycle, filled with a continuous chromatin fiber per chromosome that can take into account the heterochromatin and the euchromatin compaction. These geometries were then exported to a simulation chain using the Monte Carlo toolkit Geant4-DNA to perform computations of the physical, physicochemical, and chemical stages, in order to evaluate the influence of chromatin compaction on DSB induction and the contribution of direct and indirect damage, as well as DSB complexity. RESULTS: More direct damage and less indirect damage were observed in the heterochromatin than in the euchromatin. Nevertheless, no difference in terms of DSB complexity was observed between those formed in the heterochromatin or the euchromatin models. Yields of DSB/Gy/Gbp show an increase when both heterochromatin and euchromatin models are taken into account, compared to when only heterochromatin is considered. CONCLUSIONS: The results presented indicate that the chromatin compaction decreases DNA damage generated by ionizing radiation and thus, DNA compaction should be considered for the simulation of DNA repair and other cellular outcomes.

Assessment of developmental potential of human single pronucleated zygotes derived from conventional in vitro fertilization.

PURPOSE: The aim of this study was to non-invasively validate the developmental potential of human single pronucleated (1PN) zygotes derived from conventional in vitro fertilization (c-IVF) at the zygote stage. METHODS: Fifty 1PN zygotes derived from 45 patients undergoing c-IVF were used. Immunohistochemistry and fluorescence live cell imaging were used to confirm normal chromosome segregation during the first mitosis. The usefulness of measuring pronuclear diameter was assessed on the basis of the presence or absence of a proper first cleavage and validated by subsequent development. RESULTS: Although approximately 80% (15/19) of 1PN zygotes contained a diploid genome, immunohistochemistry revealed an unequal distribution of paternal and maternal genomes at the first mitosis. Fluorescence live imaging revealed that 73% (8/11) of 1PN zygotes formed a functional mitotic spindle at the first mitosis resulting from diploid genomes, with 25% (2/8) of these forming a tripolar spindle. 1PN zygotes in which the pronucleus disappeared and that subsequently underwent cleavage had a pronuclear diameter ≥ 32.2 µm. The selection of 1PN zygotes based on pronuclear diameter resulted in zygotes that all formed mitotic spindles with poles during cleavage. Furthermore, 63% (5/8) of these zygotes reached the blastocyst stage. CONCLUSIONS: This study demonstrates the usefulness of a non-invasive assessment of 1PN zygotes derived from c-IVF as an indicator of developmental potential. Furthermore, diploid 1PN zygotes do not always exhibit normal chromosome segregation at the first mitosis. A pronuclear diameter ≥ 32.2 µm just before PN breakdown might be a useful criterion to assess 1PN zygotes that are capable of further development.

Coadaptation of mitochondrial and nuclear genes, and the cost of mother's curse.

Strict maternal inheritance renders the mitochondrial genome susceptible to accumulating mutations that harm males, but are otherwise benign or beneficial for females. This 'mother's curse' effect can degrade male survival and fertility if unopposed by counteracting evolutionary processes. Coadaptation between nuclear and mitochondrial genomes-with nuclear genes evolving to compensate for male-harming mitochondrial substitutions-may ultimately resolve mother's curse. However, males are still expected to incur a transient fitness cost during mito-nuclear coevolution, and it remains unclear how severe such costs should be. We present a population genetic analysis of mito-nuclear coadaptation to resolve mother's curse effects, and show that the magnitude of the 'male mitochondrial load'-the negative impact of mitochondrial substitutions on male fitness components-may be large, even when genetic variation for compensatory evolution is abundant. We also find that the male load is surprisingly sensitive to population size: male fitness costs of mito-nuclear coevolution are particularly pronounced in both small and large populations, and minimized in populations of intermediate size. Our results reveal complex interactions between demography and genetic constraints during the resolution of mother's curse, suggesting potentially widespread species differences in susceptibility to mother's curse effects.

Rampant polyphyly in the Arracacia clade (Apiaceae) and an assessment of the phylogenetic utility of 20 noncoding plastid loci.

The Arracacia clade (Apiaceae, Apioideae) is a heterogeneous assemblage of 12 genera, comprising 111 known species distributed in high montane temperate and sub-alpine habitats of meso- and South America. Previous studies have indicated that the genera Arracacia, Coulterophytum, and Prionosciadium are polyphyletic, but for the most part relationships among the members of the clade are largely unknown. Initially, cladistic analyses of nrDNA ITS sequences were carried out on 212 accessions (122 taxa), representing 92 species of the Arracacia clade and outgroups from the closely-related páramo genera Cotopaxia, Niphogeton, and Perissocoeleum and members of the Perennial Endemic North American clade and its allies. Using the ITS results to inform sampling of a small subset of taxa, a pilot study examining the phylogenetic utility of 20 noncoding chloroplast loci was subsequently performed to identify those regions most useful at resolving relationships. A cost-benefit analysis determined that five loci (trnQ-5'rps16, trnD-trnT, rpl32-trnL, psbD-trnT, ndhA intron) would maximize resolution and branch support in the clade. Cladistic analyses of four of these loci (trnQ-5'rps16, trnD-trnT, rpl32-trnL, ndhA intron) and the ITS region, separately and combined, revealed that Arracacia, Coaxana, Coulterophytum, Prionosciadium, and Rhodosciadium are each polyphyletic and that Donnellsmithia and Myrrhidendron are each monophyletic. Although most relationships in the Arracacia clade and among the closely-related genera Cotopaxia, Niphogeton, and Perissocoeleum are poorly resolved and supported, ten groups are recognized for future revisionary studies. Polyploidy and rapid species radiation have likely confounded generic circumscriptions and interpretation of relationships.