First-principles study of orthorhombic CdTiO3 perovskite

In this work we perform an ab-initio study of CdTiO3 perovskite in its orthorhombic phase using FLAPW method. Our calculations help to decide between the different cristallographic structures proposed for this perovskite from X-Ray measurements. We compute the electric field gradient tensor (EFG) at Cd site and obtain excellent agreement with available experimental information from a perturbed angular correlation (PAC) experiment. We study EFG under an isotropic change of volume and show that in this case the widely used "point charge model approximation" to determine EFG works quite well.Comment: 4 pages, 1 figure. Accepted in Physical Review

Angular momentum sharing in dissipative collisions

Light charged particles emitted by the projectile-like fragment were measured in the direct and reverse collision of $^{93}$Nb and $^{116}$Sn at 25 AMeV. The experimental multiplicities of Hydrogen and Helium particles as a function of the primary mass of the emitting fragment show evidence for a correlation with net mass transfer. The ratio of Hydrogen and Helium multiplicities points to a dependence of the angular momentum sharing on the net mass transfer.Comment: 8 pages, 2 figure

Transport Properties, Thermodynamic Properties, and Electronic Structure of SrRuO3

SrRuO$_3$ is a metallic ferromagnet. Its electrical resistivity is reported for temperatures up to 1000K; its Hall coefficient for temperatures up to 300K; its specific heat for temperatures up to 230K. The energy bands have been calculated by self-consistent spin-density functional theory, which finds a ferromagnetic ordered moment of 1.45$\mu_{{\rm B}}$ per Ru atom. The measured linear specific heat coefficient $\gamma$ is 30mJ/mole, which exceeds the theoretical value by a factor of 3.7. A transport mean free path at room temperature of $\approx 10 \AA$ is found. The resistivity increases nearly linearly with temperature to 1000K in spite of such a short mean free path that resistivity saturation would be expected. The Hall coefficient is small and positive above the Curie temperature, and exhibits both a low-field and a high-field anomalous behavior below the Curie temperature.Comment: 6 pages (latex) and 6 figures (postscript, uuencoded.) This paper will appear in Phys. Rev. B, Feb. 15, 199

Half-metallic antiferromagnets in double perovskites: LaAVRuO6_6 (A=Ca, Sr, and Ba)

Based on the theoretical exploration of electronic structures, we propose that the ordered double perovskites LaAVRuO$_6$ and LaVO$_3$/ARuO$_3$ (001) superlattice (A = Ca, Sr and Ba) are strong candidates for half-metallic (HM) antiferromagnets (AFMs). %LaAVRuO$_6$ and LaVO$_3$/ARuO$_3$ have the %100% spin polarizations at the Fermi level but with zero %total magnetic moments. We have shown that the HM-AFM nature in LaAVRuO$_6$ is very robust regardless of (i) divalent ion replacement at A-sites, (ii) oxygen site relaxation, (iii) the inclusion of the Coulomb correlation, and (iv) cation disorder. A type of the double exchange interaction is expected to be responsible for the half-metallicity and the antiferromagnetism in these systems.Comment: 4 pages, 4 figure

Comparison of the Electronic Structures and Energetics of Ferroelectric LiNbO3 and LiTaO3

This paper explains the origin of the ferroelectric instability in LiNbO3 and LiTaO3 and compares the electronic structures and energetics of the two materials.Comment: 31 pages, 11 Postscript figure

The genomic basis of the plant island syndrome in Darwin’s giant daisies

The repeated, rapid and often pronounced patterns of evolutionary divergence observed in insular plants, or the ‘plant island syndrome’, include changes in leaf phenotypes, growth, as well as the acquisition of a perennial lifestyle. Here, we sequence and describe the genome of the critically endangered, Galápagos-endemic species Scalesia atractyloides Arnot., obtaining a chromosome-resolved, 3.2-Gbp assembly containing 43,093 candidate gene models. Using a combination of fossil transposable elements, k-mer spectra analyses and orthologue assignment, we identify the two ancestral genomes, and date their divergence and the polyploidization event, concluding that the ancestor of all extant Scalesia species was an allotetraploid. There are a comparable number of genes and transposable elements across the two subgenomes, and while their synteny has been mostly conserved, we find multiple inversions that may have facilitated adaptation. We identify clear signatures of selection across genes associated with vascular development, growth, adaptation to salinity and flowering time, thus finding compelling evidence for a genomic basis of the island syndrome in one of Darwin’s giant daisies

Restriction associated DNA-genotyping at multiple spatial scales in Arabidopsis lyrata reveals signatures of pathogen-mediated selection

Background: Genome scans based on outlier analyses have revolutionized detection of genes involved in adaptive processes, but reports of some forms of selection, such as balancing selection, are still limited. It is unclear whether high throughput genotyping approaches for identification of single nucleotide polymorphisms have sufficient power to detect modes of selection expected to result in reduced genetic differentiation among populations. In this study, we used Arabidopsis lyrata to investigate whether signatures of balancing selection can be detected based on genomic smoothing of Restriction Associated DNA sequencing (RAD-seq) data. We compared how different sampling approaches (both within and between subspecies) and different background levels of polymorphism (inbreeding or outcrossing populations) affected the ability to detect genomic regions showing key signatures of balancing selection, specifically elevated polymorphism, reduced differentiation and shifts towards intermediate allele frequencies. We then tested whether candidate genes associated with disease resistance (R-gene analogs) were detected more frequently in these regions compared to other regions of the genome. Results: We found that genomic regions showing elevated polymorphism contained a significantly higher density of R-gene analogs predicted to be under pathogen-mediated selection than regions of non-elevated polymorphism, and that many of these also showed evidence for an intermediate site-frequency spectrum based on Tajima’s D. However, we found few genomic regions that showed both elevated polymorphism and reduced FST among populations, despite strong background levels of genetic differentiation among populations. This suggests either insufficient power to detect the reduced population structure predicted for genes under balancing selection using sparsely distributed RAD markers, or that other forms of diversifying selection are more common for the R-gene analogs tested. Conclusions: Genome scans based on a small number of individuals sampled from a wide range of populations were sufficient to confirm the relative scarcity of signatures of balancing selection across the genome, but also identified new potential disease resistance candidates within genomic regions showing signatures of balancing selection that would be strong candidates for further sequencing efforts

Screening synteny blocks in pairwise genome comparisons through integer programming

<p>Abstract</p> <p>Background</p> <p>It is difficult to accurately interpret chromosomal correspondences such as true orthology and paralogy due to significant divergence of genomes from a common ancestor. Analyses are particularly problematic among lineages that have repeatedly experienced whole genome duplication (WGD) events. To compare multiple "subgenomes" derived from genome duplications, we need to relax the traditional requirements of "one-to-one" syntenic matchings of genomic regions in order to reflect "one-to-many" or more generally "many-to-many" matchings. However this relaxation may result in the identification of synteny blocks that are derived from ancient shared WGDs that are not of interest. For many downstream analyses, we need to eliminate weak, low scoring alignments from pairwise genome comparisons. Our goal is to objectively select subset of synteny blocks whose total scores are maximized while respecting the duplication history of the genomes in comparison. We call this "quota-based" screening of synteny blocks in order to appropriately fill a quota of syntenic relationships within one genome or between two genomes having WGD events.</p> <p>Results</p> <p>We have formulated the synteny block screening as an optimization problem known as "Binary Integer Programming" (BIP), which is solved using existing linear programming solvers. The computer program QUOTA-ALIGN performs this task by creating a clear objective function that maximizes the compatible set of synteny blocks under given constraints on overlaps and depths (corresponding to the duplication history in respective genomes). Such a procedure is useful for any pairwise synteny alignments, but is most useful in lineages affected by multiple WGDs, like plants or fish lineages. For example, there should be a 1:2 ploidy relationship between genome A and B if genome B had an independent WGD subsequent to the divergence of the two genomes. We show through simulations and real examples using plant genomes in the rosid superorder that the quota-based screening can eliminate ambiguous synteny blocks and focus on specific genomic evolutionary events, like the divergence of lineages (in cross-species comparisons) and the most recent WGD (in self comparisons).</p> <p>Conclusions</p> <p>The QUOTA-ALIGN algorithm screens a set of synteny blocks to retain only those compatible with a user specified ploidy relationship between two genomes. These blocks, in turn, may be used for additional downstream analyses such as identifying true orthologous regions in interspecific comparisons. There are two major contributions of QUOTA-ALIGN: 1) reducing the block screening task to a BIP problem, which is novel; 2) providing an efficient software pipeline starting from all-against-all BLAST to the screened synteny blocks with dot plot visualizations. Python codes and full documentations are publicly available <url>http://github.com/tanghaibao/quota-alignment</url>. QUOTA-ALIGN program is also integrated as a major component in SynMap <url>http://genomevolution.com/CoGe/SynMap.pl</url>, offering easier access to thousands of genomes for non-programmers.</p

Comparative Oncogenomic Analysis of Copy Number Alterations in Human and Zebrafish Tumors Enables Cancer Driver Discovery

The identification of cancer drivers is a major goal of current cancer research. Finding driver genes within large chromosomal events is especially challenging because such alterations encompass many genes. Previously, we demonstrated that zebrafish malignant peripheral nerve sheath tumors (MPNSTs) are highly aneuploid, much like human tumors. In this study, we examined 147 zebrafish MPNSTs by massively parallel sequencing and identified both large and focal copy number alterations (CNAs). Given the low degree of conserved synteny between fish and mammals, we reasoned that comparative analyses of CNAs from fish versus human MPNSTs would enable elimination of a large proportion of passenger mutations, especially on large CNAs. We established a list of orthologous genes between human and zebrafish, which includes approximately two-thirds of human protein-coding genes. For the subset of these genes found in human MPNST CNAs, only one quarter of their orthologues were co-gained or co-lost in zebrafish, dramatically narrowing the list of candidate cancer drivers for both focal and large CNAs. We conclude that zebrafish-human comparative analysis represents a powerful, and broadly applicable, tool to enrich for evolutionarily conserved cancer drivers.Kathy and Curt Marble Cancer Research FundArthur C. MerrillNational Institutes of Health (U.S.) (Grant CA106416)National Institutes of Health (U.S.) (Grant ROI RR020833)National Institutes of Health (U.S.) (Grant 1F32GM095213-01

Double Digest RADseq: An Inexpensive Method for De Novo SNP Discovery and Genotyping in Model and Non-Model Species

The ability to efficiently and accurately determine genotypes is a keystone technology in modern genetics, crucial to studies ranging from clinical diagnostics, to genotype-phenotype association, to reconstruction of ancestry and the detection of selection. To date, high capacity, low cost genotyping has been largely achieved via “SNP chip” microarray-based platforms which require substantial prior knowledge of both genome sequence and variability, and once designed are suitable only for those targeted variable nucleotide sites. This method introduces substantial ascertainment bias and inherently precludes detection of rare or population-specific variants, a major source of information for both population history and genotype-phenotype association. Recent developments in reduced-representation genome sequencing experiments on massively parallel sequencers (commonly referred to as RAD-tag or RADseq) have brought direct sequencing to the problem of population genotyping, but increased cost and procedural and analytical complexity have limited their widespread adoption. Here, we describe a complete laboratory protocol, including a custom combinatorial indexing method, and accompanying software tools to facilitate genotyping across large numbers (hundreds or more) of individuals for a range of markers (hundreds to hundreds of thousands). Our method requires no prior genomic knowledge and achieves per-site and per-individual costs below that of current SNP chip technology, while requiring similar hands-on time investment, comparable amounts of input DNA, and downstream analysis times on the order of hours. Finally, we provide empirical results from the application of this method to both genotyping in a laboratory cross and in wild populations. Because of its flexibility, this modified RADseq approach promises to be applicable to a diversity of biological questions in a wide range of organisms