Zika virus, an emerging flavivirus, as a cause of fever and rash in a traveller returning from Central America.

Fever in a returning traveller is a common clinical scenario for physicians in primary and acute care. Differential diagnoses for these patients are generated based on presenting clinical symptoms, travel destinations, potential exposure activities as well as the incubation period of common etiologic agents. In a case of fever and rash in a woman returning to Canada from El Salvador in November 2015, measles, dengue and chikungunya viral infections were queried as possible causes. Subsequent molecular testing using amplification of conserved regions of the flavivirus genome from nasopharyngeal and urine samples was positive, suggesting an active flavivirus infection. Sequencing was significant for the identification of Zika virus, a flavivirus that has only recently become endemic to Brazil and is now emerging throughout Central America. Zika virus should now be included in the differential diagnosis for travelers returning from Central and South America with a febrile illness and rash. To our knowledge this is the first reported case of Zika virus in Canada related to the most recent outbreak in Central America, South America and the Caribbean.

Atomic electron affinities and the role of symmetry between electron addition and subtraction in a corrected Koopmans approach.

The essential aspects of zero-temperature grand-canonical ensemble density-functional theory are reviewed in the context of spin-density-functional theory and are used to highlight the assumption of symmetry between electron addition and subtraction that underlies the corrected Koopmans approach of Tozer and De Proft (TDP) for computing electron affinities. The issue of symmetry is then investigated in a systematic study of atomic electron affinities, comparing TDP affinities with those from a conventional Koopmans evaluation and electronic energy differences. Although it cannot compete with affinities determined from energy differences, the TDP expression yields results that are a significant improvement over those from the conventional Koopmans expression. Key insight into the results from both expressions is provided by an analysis of plots of the electronic energy as a function of the number of electrons, which highlight the extent of symmetry between addition and subtraction. The accuracy of the TDP affinities is closely related to the nature of the orbitals involved in the electron addition and subtraction, being particularly poor in cases where there is a change in principal quantum number, but relatively accurate within a single manifold of orbitals. The analysis is then extended to a consideration of the ground state Mulliken electronegativity and chemical hardness. The findings further emphasize the key role of symmetry in determining the quality of the results.

Non-perturbative calculation of molecular magnetic properties within current-density functional theory.

We present a novel implementation of Kohn-Sham density-functional theory utilizing London atomic orbitals as basis functions. External magnetic fields are treated non-perturbatively, which enable the study of both magnetic response properties and the effects of strong fields, using either standard density functionals or current-density functionals-the implementation is the first fully self-consistent implementation of the latter for molecules. Pilot applications are presented for the finite-field calculation of molecular magnetizabilities, hypermagnetizabilities, and nuclear magnetic resonance shielding constants, focusing on the impact of current-density functionals on the accuracy of the results. Existing current-density functionals based on the gauge-invariant vorticity are tested and found to be sensitive to numerical details of their implementation. Furthermore, when appropriately regularized, the resulting magnetic properties show no improvement over standard density-functional results. An advantage of the present implementation is the ability to apply density-functional theory to molecules in very strong magnetic fields, where the perturbative approach breaks down. Comparison with high accuracy full-configuration-interaction results show that the inadequacies of current-density approximations are exacerbated with increasing magnetic field strength. Standard density-functionals remain well behaved but fail to deliver high accuracy. The need for improved current-dependent density-functionals, and how they may be tested using the presented implementation, is discussed in light of our findings.

High accuracy ab initio studies of electron-densities for the ground state of Be-like atomic systems.

Benchmark results for electron densities in the ground states of Li(-), Be, C(2+), Ne(6+), and Ar(14+) have been generated from very accurate variational wave functions represented in terms of extensive basis sets of exponentially correlated Gaussian functions. For Ne(6+), and Ar(14+), the upper bounds to the energies improve over previous results known from the literature. For the remaining systems our bounds are from 0.1 to 1.1 µhartree higher than the most accurate ones. We present in graphical and, partially, numerical form results both for the radial electron densities and for the difference radial density distributions (DRD) (defined with respect to the Hartree-Fock radial density) that highlight the impact of correlation effects on electron densities. Next, we have employed these DRD distributions in studies of the performance of several broadly used orbital-based quantum-chemical methods in accounting for correlation effects on the density. Our computed benchmark densities for Be have been also applied for testing the possibility of using the mathematically strict result concerning exact atomic electron densities, obtained by Ahlrichs et al. [Phys. Rev. A 23, 2106 (1981)], for the determination of the reliability range of computed densities in the long-range asymptotic region. The results obtained for Be are encouraging.

Range-dependent adiabatic connections.

Recently, we have implemented a scheme for the calculation of the adiabatic connection linking the Kohn-Sham system to the physical, interacting system. This scheme uses a generalized Lieb functional, in which the electronic interaction strength is varied in a simple linear fashion, keeping the potential or the density fixed in the process. In the present work, we generalize this scheme further to accommodate arbitrary two-electron operators, allowing the calculation of adiabatic connections following alternative paths as outlined by Yang [J. Chem. Phys. 109, 10107 (1998)]. Specifically, we examine the error-function and Gaussian-attenuated error-function adiabatic connections. It is shown that while the error-function connection displays some promising features, making it amenable to the possible development of new exchange-correlation functionals by modeling the adiabatic connection integrand, the Gaussian-attenuated error-function connection is less promising. We explore the high-density and strong static correlation regimes for two-electron systems. Implications of this work for the utility of range-separated schemes are discussed.

Accurate calculation and modeling of the adiabatic connection in density functional theory.

Using a recently implemented technique for the calculation of the adiabatic connection (AC) of density functional theory (DFT) based on Lieb maximization with respect to the external potential, the AC is studied for atoms and molecules containing up to ten electrons: the helium isoelectronic series, the hydrogen molecule, the beryllium isoelectronic series, the neon atom, and the water molecule. The calculation of AC curves by Lieb maximization at various levels of electronic-structure theory is discussed. For each system, the AC curve is calculated using Hartree-Fock (HF) theory, second-order Moller-Plesset (MP2) theory, coupled-cluster singles-and-doubles (CCSD) theory, and coupled-cluster singles-doubles-perturbative-triples [CCSD(T)] theory, expanding the molecular orbitals and the effective external potential in large Gaussian basis sets. The HF AC curve includes a small correlation-energy contribution in the context of DFT, arising from orbital relaxation as the electron-electron interaction is switched on under the constraint that the wave function is always a single determinant. The MP2 and CCSD AC curves recover the bulk of the dynamical correlation energy and their shapes can be understood in terms of a simple energy model constructed from a consideration of the doubles-energy expression at different interaction strengths. Differentiation of this energy expression with respect to the interaction strength leads to a simple two-parameter doubles model (AC-D) for the AC integrand (and hence the correlation energy of DFT) as a function of the interaction strength. The structure of the triples-energy contribution is considered in a similar fashion, leading to a quadratic model for the triples correction to the AC curve (AC-T). From a consideration of the structure of a two-level configuration-interaction (CI) energy expression of the hydrogen molecule, a simple two-parameter CI model (AC-CI) is proposed to account for the effects of static correlation on the AC. When parametrized in terms of the same input data, the AC-CI model offers improved performance over the corresponding AC-D model, which is shown to be the lowest-order contribution to the AC-CI model. The utility of the accurately calculated AC curves for the analysis of standard density functionals is demonstrated for the BLYP exchange-correlation functional and the interaction-strength-interpolation (ISI) model AC integrand. From the results of this analysis, we investigate the performance of our proposed two-parameter AC-D and AC-CI models when a simple density functional for the AC at infinite interaction strength is employed in place of information at the fully interacting point. The resulting two-parameter correlation functionals offer a qualitatively correct behavior of the AC integrand with much improved accuracy over previous attempts. The AC integrands in the present work are recommended as a basis for further work, generating functionals that avoid spurious error cancellations between exchange and correlation energies and give good accuracy for the range of densities and types of correlation contained in the systems studied here.

The calculation of adiabatic-connection curves from full configuration-interaction densities: two-electron systems.

The Lieb formulation of density-functional theory is briefly reviewed and its straightforward generalization to arbitrary electron-electron interaction strengths discussed, leading to the introduction of density-fixed and potential-fixed adiabatic connections. An iterative scheme for the calculation of the Lieb functionals under the appropriate constraints is outlined following the direct optimization approach of Wu and Yang [J. Chem. Phys. 118, 2498 (2003)]. First- and second-order optimization schemes for the calculation of accurate adiabatic-connection integrands are investigated and compared; the latter is preferred both in terms of computational efficiency and accuracy. The scheme is applicable to systems of any number of electrons. However, to determine the accuracy that may be achieved, the present work focuses on two-electron systems for which a number of simplifications may be exploited. The procedure is applied to the helium isoelectronic series and the H(2) molecule. The resulting adiabatic-connection curves yield the full configuration-interaction exchange-correlation energies extrapolated to the basis-set limit. The relationship between the Kohn-Sham and natural orbitals as functions of the electron-electron interaction strength is explored in detail for H(2). The accuracy with which the exchange-correlation contributions to the modified local potential can be determined is discussed. The new accurate adiabatic-connection curves are then compared with some recently investigated approximate forms calculated using accurate full configuration-interaction input data. This study demonstrates that the adiabatic-connection integrand may be determined accurately and efficiently, providing important insights into the link between the Kohn-Sham and traditional quantum-chemical treatments of the exchange-correlation problem in electronic-structure theory.

Maps of current density using density-functional methods.

The performance of several density-functional theory (DFT) methods for the calculation of current densities induced by a uniform magnetic field is examined. Calculations are performed using the BLYP and KT3 generalized-gradient approximations, together with the B3LYP hybrid functional. For the latter, both conventional and optimized effective potential (OEP) approaches are used. Results are also determined from coupled-cluster singles-and-doubles (CCSD) electron densities by a DFT constrained search procedure using the approach of Wu and Yang (WY). The current densities are calculated within the CTOCD-DZ2 distributed origin approach. Comparisons are made with results from Hartree-Fock (HF) theory. Several small molecules for which correlation is known to be especially important in the calculation of magnetic response properties are considered-namely, O(3), CO, PN, and H(2)CO. As examples of aromatic and antiaromatic systems, benzene and planarized cyclooctatetraene molecules are considered, with specific attention paid to the ring current phenomenon and its Kohn-Sham orbital origin. Finally, the o-benzyne molecule is considered as a computationally challenging case. The HF and DFT induced current maps show qualitative differences, while among the DFT methods the maps show a similar qualitative structure. To assess quantitative differences in the calculated current densities with different methods, the maximal moduli of the induced current densities are compared and integration of the current densities to yield shielding constants is performed. In general, the maximal modulus is reduced in moving from HF to B3LYP and BLYP, and further reduced in moving to KT3, OEP(B3LYP), and WY(CCSD). The latter three methods offer the most accurate shielding constants in comparison with both experimental and ab initio data and hence the more reliable route to DFT calculation of induced current density in molecules.

A description of the origins, design and performance of the TRAITS-SGP Atlantic salmon Salmo salar L. cDNA microarray.

The origins, design, fabrication and performance of an Atlantic salmon microarray are described. The microarray comprises 16 950 Atlantic salmon-derived cDNA features, printed in duplicate and mostly sourced from pre-existing expressed sequence tag (EST) collections [SALGENE and salmon genome project (SGP)] but also supplemented with cDNAs from suppression subtractive hybridization libraries and candidate genes involved in immune response, protein catabolism, lipid metabolism and the parr-smolt transformation. A preliminary analysis of a dietary lipid experiment identified a number of genes known to be involved in lipid metabolism. Significant fold change differences (as low as 1.2x) were apparent from the microarray analysis and were confirmed by quantitative real-time polymerase chain reaction analysis. The study also highlighted the potential for obtaining artefactual expression patterns as a result of cross-hybridization of similar transcripts. Examination of the robustness and sensitivity of the experimental design employed demonstrated the greater importance of biological replication over technical (dye flip) replication for identification of a limited number of key genes in the studied system. The TRAITS (TRanscriptome Analysis of Important Traits of Salmon)-salmon genome project microarray has been proven, in a number of studies, to be a powerful tool for the study of key traits of Atlantic salmon biology. It is now available for use by researchers in the wider scientific community.

Interferon type I and type II responses in an Atlantic salmon (Salmo salar) SHK-1 cell line by the salmon TRAITS/SGP microarray.

Interferons (IFNs) are cytokines that have proinflammatory, antiviral, and immunomodulatory effects and play a central role during a host response to pathogens. The IFN family contains both type I and type II molecules. While there are a number of type I IFNs, there is only one type II IFN. Recently both type I and type II IFN genes have been cloned in salmonid fish and recombinant proteins produced showing IFN activity. We have stimulated an Atlantic salmon cell line (SHK-1) with both type I and type II recombinant salmonid IFNs and analyzed the transcriptional response by microarray analysis. Cells were exposed to recombinant IFNs for 6 or 24 h or left unexposed as controls. RNA was hybridized to an Atlantic salmon cDNA microarray (salmon 17K feature TRAITS/SGP array) in order to assess differential gene expression in response to IFN exposure. For IFN I and II, 47 and 72 genes were stimulated, respectively; most genes were stimulated by a single IFN type, but some were affected by both IFNs, indicating coregulation of the IFN response in fish. Real-time PCR analysis was employed to confirm the microarray results for selected differentially expressed genes in both a cell line and primary leukocyte cultures.

Characterization and comparison of fatty acyl Delta6 desaturase cDNAs from freshwater and marine teleost fish species.

Fish are the most important dietary source of the n-3 highly unsaturated fatty acids (HUFA), eicosapentaenoic (EPA) and docosahexaenoic acid (DHA), that have particularly important roles in human nutrition reflecting their roles in critical physiological processes. The objective of the study described here was to clone, functionally characterize and compare expressed fatty acid desaturase genes involved in the production of EPA and DHA in freshwater and marine teleost fish species. Putative fatty acid desaturase cDNAs were isolated and cloned from common carp (Cyprinus carpio) and turbot (Psetta maximus). The enzymic activities of the products of these cDNAs, together with those of cDNAs previously cloned from rainbow trout (Oncorhynchus mykiss) and gilthead sea bream (Sparus aurata), were determined by heterologous expression in the yeast Saccharomyces cerevisiae. The carp and turbot desaturase cDNAs included open reading frames (ORFs) of 1335 and 1338 base pairs, respectively, specifying proteins of 444 and 445 amino acids. The protein sequences possessed all the characteristic features of microsomal fatty acid desaturases, including three histidine boxes, two transmembrane regions, and N-terminal cytochrome b(5) domains containing the haem-binding motif, HPGG. Functional expression showed all four fish cDNAs encode basically unifunctional Delta6 fatty acid desaturase enzymes responsible for the first and rate-limiting step in the biosynthesis of HUFA from 18:3n-3 and 18:2n-6. All the fish desaturases were more active towards the n-3 substrate with 59.5%, 31.5%, 23.1% and 7.0% of 18:3n-3 being converted to 18:4n-3 in the case of turbot, trout, sea bream and carp, respectively. The enzymes also showed very low, probably physiologically insignificant, levels of Delta5 desaturase activity, but none of the products showed Delta4 desaturase activity. The cloning and characterization of desaturases from these fish is an important advance, as they are species in which there is a relative wealth of data on the nutritional regulation of fatty acid desaturation and HUFA synthesis, and between which substantive differences occur.

Outcome of women with inadequate cervical smears followed up for five years.

BACKGROUND: The clinical and prognostic significance of "inadequate" cervical smear is unknown, even though women with repeated inadequate smears are referred for colposcopy in the National Health Service (NHS) Cervical Screening Programme. AIM: To follow up a cohort of women with inadequate cervical smears over the following five years to examine outcomes, including detection of high grade cervical intraepithelial neoplasia (CIN). METHODS: The study comprised 1972 women with an inadequate cervical smear reported at Walsall Hospitals NHS Trust between 1 April 1995 and 31 March 1996. Results of cervical smears and biopsies taken over the following five years were reviewed to confirm the outcome. FINDINGS: Within five years, 2.2% of women with an inadequate cervical smear developed histologically confirmed high grade CIN, which was higher than the 1.3% seen among all women with cervical smear tests reported at the same laboratory over the same period, although the difference was not significant at the 95% level of confidence. Where inadequacy resulted from or was contributed to by "polymorphs obscuring", the risk of subsequent development/detection of high grade CIN was 2.6%. CONCLUSIONS: Women with inadequate cervical smears had an increased risk of detection of high grade CIN in the following five years compared with "all women". This increased risk was not significant, although if a larger number of women had been studied significance may have been reached, so that further studies are needed. The increased risk appeared to be at least partially dependent on the reason for inadequacy.

Mapping of quantitative trait loci controlling trypanotolerance in a cross of tolerant West African N'Dama and susceptible East African Boran cattle.

Trypanosomosis, or sleeping sickness, is a major disease constraint on livestock productivity in sub-Saharan Africa. To identify quantitative trait loci (QTL) controlling resistance to trypanosomosis in cattle, an experimental cross was made between trypanotolerant African N'Dama (Bos taurus) and trypanosusceptible improved Kenya Boran (Bos indicus) cattle. Sixteen phenotypic traits were defined describing anemia, body weight, and parasitemia. One hundred seventy-seven F2 animals and their parents and grandparents were genotyped at 477 molecular marker loci covering all 29 cattle autosomes. Total genome coverage was 82%. Putative QTL were mapped to 18 autosomes at a genomewise false discovery rate of <0.20. The results are consistent with a single QTL on 17 chromosomes and two QTL on BTA16. Individual QTL effects ranged from approximately 6% to 20% of the phenotypic variance of the trait. Excluding chromosomes with ambiguous or nontrypanotolerance effects, the allele for resistance to trypanosomosis originated from the N'Dama parent at nine QTL and from the Kenya Boran at five QTL, and at four QTL there is evidence of an overdominant mode of inheritance. These results suggest that selection for trypanotolerance within an F2 cross between N'Dama and Boran cattle could produce a synthetic breed with higher trypanotolerance levels than currently exist in the parental breeds.

Identification of mechanisms of natural resistance to African trypanosomiasis in cattle.

Natural resistance to African trypanosomiasis in certain Bos taurus cattle in West Africa, called trypanotolerance, may hold solutions for control of this economically crippling disease. Comparison of immune responses between trypanotolerant and trypanosusceptible cattle have shown some differences in antibody response, complement level and cytokine expression, but it is not known whether these differences are the cause of resistance. Two experiments were carried out to assess the contribution of the immune and haemopoietic systems to trypanotolerance. The production of haemopoietic chimaeras from trypanotolerant and susceptible twin calves and comparison of their responses after infection with singleton calves, allowed an assessment of the role of the haemopoietic system in trypanotolerance. An in vivo depletion of CD4 cells in the two breeds allowed an appraisal of the role of T and B lymphocytes in trypanotolerance. The results of the two experiments suggest that natural resistance comprises at least two mechanisms, an innate mechanism that controls parasite growth, and another, involving the haemopoietic system, that is able to limit anaemia. This supports the hypothesis that innate mechanisms in trypanotolerant cattle are more efficient in controlling disease, making them less reliant on antibody responses.

A vertebrate fatty acid desaturase with Delta 5 and Delta 6 activities.

Delta5 and Delta6 fatty acid desaturases are critical enzymes in the pathways for the biosynthesis of the polyunsaturated fatty acids arachidonic, eicosapentaenoic, and docosahexaenoic acids. They are encoded by distinct genes in mammals and Caenorhabditis elegans. This paper describes a cDNA isolated from zebrafish (Danio rerio) with high similarity to mammalian Delta6 desaturase genes. The 1,590-bp sequence specifies a protein that, in common with other fatty acid desaturases, contains an N-terminal cytochrome b(5) domain and three histidine boxes, believed to be involved in catalysis. When the zebrafish cDNA was expressed in Saccharomyces cerevisiae it conferred the ability to convert linoleic acid (18:2n-6) and alpha-linolenic acid (18:3n-3) to their corresponding Delta6 desaturated products, 18:3n-6 and 18:4n-3. However, in addition it conferred on the yeast the ability to convert di-homo-gamma-linoleic acid (20:3n-6) and eicosatetraenoic acid (20:4n-3) to arachidonic acid (20:4n-6) and eicosapentaenoic acid (20:5n-3), respectively, indicating that the zebrafish gene encodes an enzyme having both Delta5 and Delta6 desaturase activity. The zebrafish Delta5/Delta6 desaturase may represent a component of a prototypic vertebrate polyunsaturated fatty acids biosynthesis pathway.

Susceptibility of tumour necrosis factor-alpha genetically deficient mice to Trypanosoma congolense infection.

The TNF-alpha gene on mouse chromosome MMU17 is among the candidates for the trypanosomosis resistance QTL Tir1. Tir1 has the largest effect of those loci so far detected which influence degree of resistance to murine trypanosomosis caused by Trypanosoma congolense infection. We therefore studied the survival to 180 days after challenge with T. congolense of mice that were homozygous and hemizygous with respect to a disruption of the TNF-alpha gene on a > 99% C57BL/6 (resistant) background. We also examined the responses of TNF-alpha hemizygous mice produced by crossing the deletion line with mice of the C57BL/6J strain, and with mice of the susceptible A/J strain. Mice lacking a functional TNF-alpha gene were shown to be highly susceptible to challenge with T. congolense with a median survival time of 37 days. This was comparable to 71 days for control wild-type mice, and 61 and 111 days for mice of the susceptible A/J and resistant C57BL/6J strains, respectively. In mice of the deletion line, the C57BL/6 TNF-alpha allele tended to be dominant to the TNF knockout in terms of resistance. We conclude that TNF-alpha plays an important role in resistance to the effects of T. congolense infection in mice.

Fine mapping of trypanosomiasis resistance loci in murine advanced intercross lines.

We have previously reported the results of genome-wide searches in two murine F(2) populations for QTLs that influence survival following Trypanosoma congolense infection. Three loci, Tir1, Tir2, and Tir3, were identified and mapped to mouse Chromosomes (Chrs) 17, 5, and 1 respectively, with confidence intervals (CIs) in the range 10-40 cM. The size of these CIs is to a large degree the consequence of limited numbers of recombination events in small chromosomal regions in F(2) populations. A number of population designs have been proposed to increase recombination levels in crosses, one of which is the advanced intercross line (AIL). Here we report fine mapping of Tir1, Tir2, and Tir3 in G6 populations of two independent murine AILs created by crossing the C57BL/6J strain with the A/J and BALB/cJ strains, respectively. Data were analyzed by two methods that gave equally informative and similar results. The three QTLs were confirmed in the A/J x C57BL/6J AIL and in the combined data set, but Tir2 was apparently lost from the BALB/cJ x C57BL/6J AIL. The reduction in CIs for the Tir loci ranged from 2.5 to more than ten-fold in G6 populations by comparison with CIs obtained previously in the equivalent F(2) generations. Mapping in the AILs also resolved the Tir3 locus into three trypanosomiasis resistance QTLs, revealing a degree of complexity not evident in extensive studies at the F(2) level.

Evidence for genomic imprinting of the major QTL controlling susceptibility to trypanosomiasis in mice.

Inbred strains of laboratory mice exhibit marked differences in survival time following infection with Trypanosoma congolense, the principal cause of trypanosomiasis in African livestock. The difference in survival time between the relatively resistant C57BL/6 J and more susceptible BALB/c inbred strains has been attributed to three quantitative trait loci (QTLs), Tir1, Tir2 and Tir3. In order to determine whether there was a parent-of-origin effect on this trait, four backcross populations derived from the C57BL/6 J and BALB/c parental strains were bred and inoculated with T. congolense. The two populations with F1 fathers and BALB/c mothers had a significantly greater overall survival rate than the two populations with BALB/c fathers and F1 mothers. This pattern of inheritance suggested the involvement of imprinted genes. Genotyping with markers at the three QTLs controlling susceptibility revealed that the difference in survival time was consistent with genomic imprinting of the QTL of largest effect, Tir1.

Strategies to optimize marker-assisted introgression of multiple unlinked QTL.

To optimize designs to implement marker-assisted introgression programs aiming to introgress three unlinked quantitative trait loci (QTL), the present paper studies different alternatives versus a traditional backcross or intercross phase. Four alternative backcross strategies appear to be more advantageous by having 50% less genotyping load than a traditional backcross strategy tracking all three QTL at a time through a single line. A multiplication phase following the selection of homozygous animals at the three QTL as an intercross alternative allows doubling of the number of homozygous animals in a mouse model compared with the first intercross generation. Within the same model, a second intercross alternative with individuals carrying all three QTL at the first intercross results in a 12-fold increase in the number of homozygous animals obtained in the first intercross generation. The same ranges of decrease are observed in the number of animals to be genotyped and the number of genotypings when targeting a fixed number of homozygous animals. An option, with two lines each carrying two QTL through the backcross phase and coupled with the second intercross alternative, appears to be the best introgression alternative. This option requires 76% fewer genotypings, 68% fewer animals to be genotyped, and costs 75% less than an option in which all three QTL are introgressed through a single line.