Phylogenetic position of the acariform mites: sensitivity to homology assessment under total evidence

Background: Mites (Acari) have traditionally been treated as monophyletic, albeit composed of two major lineages: Acariformes and Parasitiformes. Yet recent studies based on morphology, molecular data, or combinations thereof, have increasingly drawn their monophyly into question. Furthermore, the usually basal (molecular) position of one or both mite lineages among the chelicerates is in conflict to their morphology, and to the widely accepted view that mites are close relatives of Ricinulei. Results: The phylogenetic position of the acariform mites is examined through employing SSU, partial LSU sequences, and morphology from 91 chelicerate extant terminals (forty Acariformes). In a static homology framework, molecular sequences were aligned using their secondary structure as guide, whereby regions of ambiguous alignment were discarded, and pre-aligned sequences analyzed under parsimony and different mixed models in a Bayesian inference. Parsimony and Bayesian analyses led to trees largely congruent concerning infraordinal, well-supported branches, but with low support for inter-ordinal relationships. An exception is Solifugae + Acariformes (P. P = 100%, J. = 0.91). In a dynamic homology framework, two analyses were run: a standard POY analysis and an analysis constrained by secondary structure. Both analyses led to largely congruent trees; supporting a (Palpigradi (Solifugae Acariformes)) clade and Ricinulei as sister group of Tetrapulmonata with the topology (Ricinulei (Amblypygi (Uropygi Araneae))). Combined analysis with two different morphological data matrices were run in order to evaluate the impact of constraining the analysis on the recovered topology when employing secondary structure as a guide for homology establishment. The constrained combined analysis yielded two topologies similar to the exclusively molecular analysis for both morphological matrices, except for the recovery of Pedipalpi instead of the (Uropygi Araneae) clade. The standard (direct optimization) POY analysis, however, led to the recovery of trees differing in the absence of the otherwise well-supported group Solifugae + Acariformes. Conclusions: Previous studies combining ribosomal sequences and morphology often recovered topologies similar to purely morphological analyses of Chelicerata. The apparent stability of certain clades not recovered here, like Haplocnemata and Acari, is regarded as a byproduct of the way the molecular homology was previously established using the instrumentalist approach implemented in POY. Constraining the analysis by a priori homology assessment is defended here as a way of maintaining the severity of the test when adding new data to the analysis. Although the strength of the method advocated here is keeping phylogenetic information from regions usually discarded in an exclusively static homology framework; it still has the inconvenience of being uninformative on the effect of alignment ambiguity on resampling methods of clade support estimation. Finally, putative morphological apomorphies of Solifugae + Acariformes are the reduction of the proximal cheliceral podomere, medial abutting of the leg coxae, loss of sperm nuclear membrane, and presence of differentiated germinative and secretory regions in the testis delivering their products into a common lumen.

Desenvolupament computacional de la semblança molecular quàntica

La present Tesi Doctoral, titulada desenvolupament computacional de la semblança molecular quàntica, tracta, fonamentalment, els aspectes de càlcul de mesures de semblança basades en la comparació de funcions de densitat electrònica.El primer capítol, Semblança quàntica, és introductori. S'hi descriuen les funcions de densitat de probabilitat electrònica i llur significança en el marc de la mecànica quàntica. Se n'expliciten els aspectes essencials i les condicions matemàtiques a satisfer, cara a una millor comprensió dels models de densitat electrònica que es proposen. Hom presenta les densitats electròniques, mencionant els teoremes de Hohenberg i Kohn i esquematitzant la teoria de Bader, com magnituds fonamentals en la descripció de les molècules i en la comprensió de llurs propietats.En el capítol Models de densitats electròniques moleculars es presenten procediments computacionals originals per l'ajust de funcions densitat a models expandits en termes de gaussianes 1s centrades en els nuclis. Les restriccions físico-matemàtiques associades a les distribucions de probabilitat s'introdueixen de manera rigorosa, en el procediment anomenat Atomic Shell Approximation (ASA). Aquest procediment, implementat en el programa ASAC, parteix d'un espai funcional quasi complert, d'on se seleccionen variacionalment les funcions o capes de l'expansió, d'acord als requisits de no negativitat. La qualitat d'aquestes densitats i de les mesures de semblança derivades es verifica abastament. Aquest model ASA s'estén a representacions dinàmiques, físicament més acurades, en quant que afectades per les vibracions nuclears, cara a una exploració de l'efecte de l'esmorteïment dels pics nuclears en les mesures de semblança molecular. La comparació de les densitats dinàmiques respecte les estàtiques evidencia un reordenament en les densitats dinàmiques, d'acord al que constituiria una manifestació del Principi quàntic de Le Chatelier. El procediment ASA, explícitament consistent amb les condicions de N-representabilitat, s'aplica també a la determinació directe de densitats electròniques hidrogenoides, en un context de teoria del funcional de la densitat.El capítol Maximització global de la funció de semblança presenta algorismes originals per la determinació de la màxima sobreposició de les densitats electròniques moleculars. Les mesures de semblança molecular quàntica s'identifiquen amb el màxim solapament, de manera es mesuri la distància entre les molècules, independentment dels sistemes de referència on es defineixen les densitats electròniques. Partint de la solució global en el límit de densitats infinitament compactades en els nuclis, es proposen tres nivells de aproximació per l'exploració sistemàtica, no estocàstica, de la funció de semblança, possibilitant la identificació eficient del màxim global, així com també dels diferents màxims locals. Es proposa també una parametrització original de les integrals de recobriment a través d'ajustos a funcions lorentzianes, en quant que tècnica d'acceleració computacional. En la pràctica de les relacions estructura-activitat, aquests avenços possibiliten la implementació eficient de mesures de semblança quantitatives, i, paral·lelament, proporcionen una metodologia totalment automàtica d'alineació molecular. El capítol Semblances d'àtoms en molècules descriu un algorisme de comparació dels àtoms de Bader, o regions tridimensionals delimitades per superfícies de flux zero de la funció de densitat electrònica. El caràcter quantitatiu d'aquestes semblances possibilita la mesura rigorosa de la noció química de transferibilitat d'àtoms i grups funcionals. Les superfícies de flux zero i els algorismes d'integració usats han estat publicats recentment i constitueixen l'aproximació més acurada pel càlcul de les propietats atòmiques. Finalment, en el capítol Semblances en estructures cristal·lines hom proposa una definició original de semblança, específica per la comparació dels conceptes de suavitat o softness en la distribució de fonons associats a l'estructura cristal·lina. Aquests conceptes apareixen en estudis de superconductivitat a causa de la influència de les interaccions electró-fonó en les temperatures de transició a l'estat superconductor. En aplicar-se aquesta metodologia a l'anàlisi de sals de BEDT-TTF, s'evidencien correlacions estructurals entre sals superconductores i no superconductores, en consonància amb les hipòtesis apuntades a la literatura sobre la rellevància de determinades interaccions.Conclouen aquesta tesi un apèndix que conté el programa ASAC, implementació de l'algorisme ASA, i un capítol final amb referències bibliogràfiques.

Molecular characterization of Iranian wheat stripe virus shows its taxonomic position as a distinct species in the genus Tenuivirus

The full lengths of three genome segments of Iranian wheat stripe virus (IWSV) were amplified by reverse transcription (RT) followed by polymerase chain reaction (PCR) using a primer complementary to tenuivirus conserved terminal sequences. The segments were sequenced and found to comprise 3469, 2337, and 1831 nt, respectively. The gene organization of these segments is similar to that of other known tenuiviruses, each displaying an ambisense coding strategy. IWSV segments, however, are different from those of other viruses with respect to the number of nucleotides and deduced amino acid sequence for each ORF. Depending on the segment, the first 16-22 nt at the 5' end and the first 16 nt at the 3' end are highly conserved among IWSV and rice hoja blanca virus (RHBV), rice stripe virus (RSV) and maize stripe virus ( MStV). In addition, the first 15-18 nt at the 5' end are complementary to the first 16-18 nt at the 3' end. Phylogenetic analyses showed close similarity and a common ancestor for IWSV, RHBV, and Echinochloa hoja blanca virus (EHBV). These findings confirm the position of IWSV as a distinct species in the genus Tenuivirus.

Molecular characterization of Iranian wheat stripe virus shows its taxonomic position as a distinct species in the genus Tenuivirus

The full lengths of three genome segments of Iranian wheat stripe virus (IWSV) were amplified by reverse transcription (RT) followed by polymerase chain reaction (PCR) using a primer complementary to tenuivirus conserved terminal sequences. The segments were sequenced and found to comprise 3469, 2337, and 1831 nt, respectively. The gene organization of these segments is similar to that of other known tenuiviruses, each displaying an ambisense coding strategy. IWSV segments, however, are different from those of other viruses with respect to the number of nucleotides and deduced amino acid sequence for each ORF. Depending on the segment, the first 16-22 nt at the 5' end and the first 16 nt at the 3' end are highly conserved among IWSV and rice hoja blanca virus (RHBV), rice stripe virus (RSV) and maize stripe virus ( MStV). In addition, the first 15-18 nt at the 5' end are complementary to the first 16-18 nt at the 3' end. Phylogenetic analyses showed close similarity and a common ancestor for IWSV, RHBV, and Echinochloa hoja blanca virus (EHBV). These findings confirm the position of IWSV as a distinct species in the genus Tenuivirus.

Croton laceratoglandulosus (Euphorbiaceae s.s.), a new glandular-stipulate species from Brazil and Bolivia, and its systematic position based on molecular analysis

Croton laceratoglandulosus, a new species from the dry forests of the Brazilian states of Piaui, Ceara, Bahia and Minas Gerais, and the Bolivian department of Santa Cruz, is described and illustrated here. Molecular sequence data demonstrate that it is most closely related to the taxa of Croton section Cascarilla, and not to sections Medea or Barhamia, which also have glandular calyces and laciniate stipules. (C) 2008 The Linnean Society of London, Botanical Journal of the Linnean Society, 2008, 158, 493-498.

Multiple nucleolus organizer regions in Leptodactylus mystacinus (Amphibia, Anura) and comments on its systematic position in the L fuscus group based on cytogenetic and molecular analyses

Specimens of Leptodactylus mystacinus from Brazil were karyotyped with conventional and differential staining. The 2n = 22 karyotype is similar to that found for the majority of the Leptodactylus, the karyotypic conservatism also confirmed by the similarity of the replication banding patterns with those previously described. L. mystacinus has a small amount of C-banded heterochromatin, located mainly at the centromeres, although telomeric or interstitial bands have also been noticed. With DA/CMA(3) some chromosome regions showed slightly bright fluorescence, and with DA/DAPI, no particular AT-rich repetitive region was observed. Silver staining showed an extensive inter- and intraindividual variation in the number and position of Ag-positive regions, in 1p, 4p, 8p, 8q, and 11p. Nevertheless, FISH using rDNA probes confirmed only the signals on the short arms of chromosomes 4 and 8 as true NORs. The remaining silver stained regions are probably due to the heterochromatin with some affinity to the Ag-staining. Phylogenetic analysis based on partial cytochrome b sequence revealed that L. mystacinus forms a basal branch, so that the presence of multiple NORs in pairs 4 and 8 in this species indicates an autapomorphy.

A New Miniature Characid (Ostariophysi: Characiformes: Characidae), with Phylogenetic Position Inferred from Morphological and Molecular Data

Erythrocharax altipinnis is described from the Serra do Cachimbo, Pará, Brazil. The new taxon is distinguished from all of the Characidae genera by having the pelvic bones firmly attached through the isquiatic processes; a nearly triangular hiatus in the musculature covering the anterior chamber of the swim bladder between the first and second pleural ribs (pseudotympanum); the pedunculate, notably expanded and distally compressed teeth in both jaws; circumorbital series represented by antorbital and four infraorbital bones with laterosensory canals not enclosed; a single tooth row in the premaxillary with the teeth perfectly aligned and similar in shape and cusp number; the first three branched dorsal-fin rays distinctly elongate in males; a bright red adipose and caudal fins in life; a conspicuous dark midlateral stripe extending from the opercle to the tip of the median caudal-fin rays; and by the absence of a humeral spot. The phylogenetic position of the new taxon is discussed using morphological and molecular datasets, with conflicting results of both approaches discussed. Additionally, a summarized discussion on the current problems in the Characidae taxonomy is presented and the principal biases in the morphological dataset are also discussed. © 2013 Netto-Ferreira et al.

Molecular phylogenetic inferences on the position of the Mollusca within the Lophotrochozoa

Die phylogenetische Position der Mollusken innerhalb der Trochozoa sowie die interne Evolution der Klassen der Mollusca sind weitgehend unbekannt und wurden in meiner Arbeit anhand molekularer Merkmale untersucht. Phylogenomische Analysen zeigten in der Vergangenheit eine gute Auflösung für ursprüngliche Speziationsereignisse. Daher wurden hier drei neue EST Datensätze generiert: für Sipunculus nudus (Sipuncula), Barentsia elongata (Kamptozoa) und Lepidochitona cinerea, (Polyplacophora, Mollusca). Zusätzlich wurden gezielt Gene verschiedener Mollusken mittels RT-PCR amplifiziert. rnSowohl Kamptozoen als auch Sipunculiden wurden aufgrund morphologischer Kriterien bisher als mögliche Schwestergruppe der Mollusken gehandelt, aber die hier erzielten Ergebnisse zur Evolution der Hämerythrine, Gen-Anordnungen der mitochondrialen Genome und phylogenetische Analysen der ribosomalen und der mitochondriellen Proteine stützen diese Hypothese nicht. Die Position der Kamptozoa erwies sich hier generell als unbeständig; phylogenomische Analysen deuten eine Nähe zu den Bryozoen an, aber diese Position wird stark durch die Auswahl der Taxa beeinflusst. Dagegen weisen meine Analysen klar auf eine nähere Beziehung zwischen Annelida und Sipuncula hin. Die ribosomalen Proteine zeigen Sipuncula (und Echiura) sogar als Subtaxa der Anneliden. Wie den Mollusken fehlt den Sipunculiden jegliche Segmentierung und meine Ergebnisse legen hier die Möglichkeit des Verlusts dieses Merkmals innerhalb der Anneliden bei den Sipunculiden nahe. Innerhalb der Mollusken wurden die Solenogastren bereits als Schwestergruppe aller rezenten Mollusken vorgeschlagen. Im Rahmen meiner Arbeit wurden von drei verschiedenen Solenogastren-Arten die ersten zuverlässigen 18S rRNA-Sequenzen ermittelt, und es zeigte sich, dass alle bisher veröffentlichten 18S-Sequenzen dieser Molluskenklasse höchst unvollständig oder fehlerhaft sind. rnRibosomale Proteine sind gute phylogenetische Marker und hier wurden die Auswahl und Anzahl dieser Gene für phylogenetische Analysen optimiert. Über Sonden-basierte Detektion wurde eine sampling-Strategie getestet, die im Vergleich mit standard-phylogenomischen Ansätzen zukünftige molekulare Stammbaumrekonstruktionen mit größerem Taxonsampling ermöglicht.rn

Molecular studies suggest that cartilaginous fishes have a terminal position in the piscine tree

The Chondrichthyes (cartilaginous fishes) are commonly accepted as being sister group to the other extant Gnathostomata (jawed vertebrates). To clarify gnathostome relationships and to aid in resolving and dating the major piscine divergences, we have sequenced the complete mtDNA of the starry skate and have included it in phylogenetic analysis along with three squalomorph chondrichthyans—the common dogfish, the spiny dogfish, and the star spotted dogfish—and a number of bony fishes and amniotes. The direction of evolution within the gnathostome tree was established by rooting it with the most closely related non-gnathostome outgroup, the sea lamprey, as well as with some more distantly related taxa. The analyses placed the chondrichthyans in a terminal position in the piscine tree. These findings, which also suggest that the origin of the amniote lineage is older than the age of the oldest extant bony fishes (the lungfishes), challenge the evolutionary direction of several morphological characters that have been used in reconstructing gnathostome relationships. Applying as a calibration point the age of the oldest lungfish fossils, 400 million years, the molecular estimate placed the squalomorph/batomorph divergence at ≈190 million years before present. This dating is consistent with the occurrence of the earliest batomorph (skates and rays) fossils in the paleontological record. The split between gnathostome fishes and the amniote lineage was dated at ≈420 million years before present.

The four hexamerin genes in the honey bee: structure, molecular evolution and function deduced from expression patterns in queens, workers and drones

Background: Hexamerins are hemocyanin-derived proteins that have lost the ability to bind copper ions and transport oxygen; instead, they became storage proteins. The current study aimed to broaden our knowledge on the hexamerin genes found in the honey bee genome by exploring their structural characteristics, expression profiles, evolution, and functions in the life cycle of workers, drones and queens. Results: The hexamerin genes of the honey bee (hex 70a, hex 70b, hex 70c and hex 110) diverge considerably in structure, so that the overall amino acid identity shared among their deduced protein subunits varies from 30 to 42%. Bioinformatics search for motifs in the respective upstream control regions (UCRs) revealed six overrepresented motifs including a potential binding site for Ultraspiracle (Usp), a target of juvenile hormone (JH). The expression of these genes was induced by topical application of JH on worker larvae. The four genes are highly transcribed by the larval fat body, although with significant differences in transcript levels, but only hex 110 and hex 70a are re-induced in the adult fat body in a caste-and sex-specific fashion, workers showing the highest expression. Transcripts for hex 110, hex 70a and hex70b were detected in developing ovaries and testes, and hex 110 was highly transcribed in the ovaries of egg-laying queens. A phylogenetic analysis revealed that HEX 110 is located at the most basal position among the holometabola hexamerins, and like HEX 70a and HEX 70c, it shares potential orthology relationship with hexamerins from other hymenopteran species. Conclusions: Striking differences were found in the structure and developmental expression of the four hexamerin genes in the honey bee. The presence of a potential binding site for Usp in the respective 5' UCRs, and the results of experiments on JH level manipulation in vivo support the hypothesis of regulation by JH. Transcript levels and patterns in the fat body and gonads suggest that, in addition to their primary role in supplying amino acids for metamorphosis, hexamerins serve as storage proteins for gonad development, egg production, and to support foraging activity. A phylogenetic analysis including the four deduced hexamerins and related proteins revealed a complex pattern of evolution, with independent radiation in insect orders.

The taxonomic status of the endangered thin-spined porcupine, Chaetomys subspinosus (Olfers, 1818), based on molecular and karyologic data

Background: The thin-spined porcupine, also known as the bristle-spined rat, Chaetomys subspinosus (Olfers, 1818), the only member of its genus, figures among Brazilian endangered species. In addition to being threatened, it is poorly known, and even its taxonomic status at the family level has long been controversial. The genus Chaetomys was originally regarded as a porcupine in the family Erethizontidae, but some authors classified it as a spiny-rat in the family Echimyidae. Although the dispute seems to be settled in favor of the erethizontid advocates, further discussion of its affinities should be based on a phylogenetic framework. In the present study, we used nucleotide-sequence data from the complete mitochondrial cytochrome b gene and karyotypic information to address this issue. Our molecular analyses included one individual of Chaetomys subspinosus from the state of Bahia in northeastern Brazil, and other hystricognaths. Results: All topologies recovered in our molecular phylogenetic analyses strongly supported Chaetomys subspinosus as a sister clade of the erethizontids. Cytogenetically, Chaetomys subspinosus showed 2n = 52 and FN = 76. Although the sexual pair could not be identified, we assumed that the X chromosome is biarmed. The karyotype included 13 large to medium metacentric and submetacentric chromosome pairs, one small subtelocentric pair, and 12 small acrocentric pairs. The subtelocentric pair 14 had a terminal secondary constriction in the short arm, corresponding to the nucleolar organizer region (Ag-NOR), similar to the erethizontid Sphiggurus villosus, 2n = 42 and FN = 76, and different from the echimyids, in which the secondary constriction is interstitial. Conclusion: Both molecular phylogenies and karyotypical evidence indicated that Chaetomys is closely related to the Erethizontidae rather than to the Echimyidae, although in a basal position relative to the rest of the Erethizontidae. The high levels of molecular and morphological divergence suggest that Chaetomys belongs to an early radiation of the Erethizontidae that may have occurred in the Early Miocene, and should be assigned to its own subfamily, the Chaetomyinae.

Molecular cloning reveals that the p160 myb-binding protein is a novel, predominantly nucleolar protein which may play a role in transactivation by Myb

We have previously detected two related murine nuclear proteins, p160 and p67, that can bind to the leucine zipper motif within the negative regulatory domain of the Myb transcription factor. We now describe the molecular cloning of cDNA corresponding to murine p160. The P160 gene is located on mouse chromosome 11, and related sequences are found on chromosomes 1 and 12. The predicted p160 protein is novel, and in agreement with previous studies, we find that the corresponding 4.5-kb mRNA is ubiquitously expressed. We showed that p67 is an N-terminal fragment of p160 which is generated by proteolytic cleavage in certain cell types. The protein encoded by the cloned p160 cDNA and an engineered protein (p67*) comprising the amino-terminal region of p160 exhibit binding specificities for the Myb and Jun leucine zipper regions identical to those of endogenous p160 and p67, respectively. This implies that the Myb-binding site of p160 lies within the N-terminal 580 residues and that the Jun-binding site is C-terminal to this position. Moreover, we show that p67* but not p160 can inhibit transactivation by Myb. Unexpectedly, immunofluorescence studies show that p160 is localized predominantly in the nucleolus. The implications of these results for possible functions of p160 are discussed.

Small molecular probes for G-protein-coupled C5a receptors: Conformationally constrained antagonists derived from the C terminus of the human plasma protein C5a

Activation of the human complement system of plasma proteins in response to infection or injury produces a 4-helix bundle glycoprotein (74 amino acids) known as C5a. C5a binds to G-protein-coupled receptors on cell surfaces triggering receptor-ligand internalization, signal transduction, and powerful inflammatory responses. Since excessive levels of C5a are associated with autoimmune and chronic inflammatory disorders, inhibitors of receptor activation may have therapeutic potential. We now report solution structures and receptor-binding and antagonist activities for some of the first small molecule antagonists of C5a derived from its hexapeptide C terminus. The antagonist NMe-Phe-Lys-Pro-D-Cha-Trp-D-Arg-CO2H (1) surprisingly shows an unusually well-defined solution structure as determined by H-1 NMR spectroscopy. This is one of the smallest acyclic peptides found to possess a defined solution conformation, which can be explained by the constraining role of intramolecular hydrogen bonding. NOE and coupling constant data, slow deuterium exchange, and a low dependence on temperature for the chemical shift of the D-Cha-NH strongly indicate an inverse gamma turn stabilized by a D-Cha-NH ... OC-Lys hydrogen bond. Smaller conformational populations are associated with a hydrogen bond between Trp-NH ... OC-Lys, defining a type II beta turn distorted by the inverse gamma turn incorporated within it. An excellent correlation between receptor-affinity and antagonist activity is indicated for a limited set of synthetic peptides. Conversion of the C-terminal carboxylate of 1 to an amide decreases antagonist potency 5-fold, but potency is increased up to 10-fold over 1 if the amide bond is made between the C-terminal carboxylate and a Lys/Orn side chain to form a cyclic analogue. The solution structure of cycle 6 also shows gamma and beta turns; however, the latter occurs in a different position, and there are clear conformational changes in 6 vs 1 that result in enhanced activity. These results indicate that potent C5a antagonists can be developed by targeting site 2 alone of the C5a receptor and define a novel pharmacophore for developing powerful receptor probes or drug candidates.

Molecular cloning and chromosomal mapping of the human homologue of MYB binding protein (P160) 1A (MYBBP1A) to 17p13.3

We have previously isolated and characterized murine MYB binding protein (p160) 1a, a protein that specifically interacts with the leucine zipper motif within the negative regulatory domain of the c-Myb proto-oncoprotein, We now describe the molecular cloning of the human MYBBP1A cDNA and chromosomal localization to 17p13.3 by fluorescence in situ hybridization analysis, Given the likely presence of a tumor suppressor gene (or genes) within this region of chromosome 17, the position of MYBBP1A was further mapped by radiation hybrid analysis and was found to lie between markers D17S1828 and D17S938. A P1 artificial chromosome clone containing the 5' region of MYBBP1A was isolated and indicates a physical linkage between MYBBP1A and the 15-lipoxygenase gene (ALOX15), A novel, polymorphic (CA)(25) dinucleotide repeat was also isolated from this PAC and may serve as a useful marker for MYBBP1A and this region of chromosome 17. (C) 1999 Academic Press.

Molecular orbital calculations of two-electron states for P-donor solid-state spin qubits

We theoretically study the Hilbert space structure of two neighboring P-donor electrons in silicon-based quantum computer architectures. To use electron spins as qubits, a crucial condition is the isolation of the electron spins from their environment, including the electronic orbital degrees of freedom. We provide detailed electronic structure calculations of both the single donor electron wave function and the two-electron pair wave function. We adopted a molecular orbital method for the two-electron problem, forming a basis with the calculated single donor electron orbitals. Our two-electron basis contains many singlet and triplet orbital excited states, in addition to the two simple ground state singlet and triplet orbitals usually used in the Heitler-London approximation to describe the two-electron donor pair wave function. We determined the excitation spectrum of the two-donor system, and study its dependence on strain, lattice position, and interdonor separation. This allows us to determine how isolated the ground state singlet and triplet orbitals are from the rest of the excited state Hilbert space. In addition to calculating the energy spectrum, we are also able to evaluate the exchange coupling between the two donor electrons, and the double occupancy probability that both electrons will reside on the same P donor. These two quantities are very important for logical operations in solid-state quantum computing devices, as a large exchange coupling achieves faster gating times, while the magnitude of the double occupancy probability can affect the error rate.