Genetic Modification of Recombination Rate in TRIBOLIUM CASTANEUM

Asymmetrical responses were obtained in a replicated study of 15 generations of two-way selection for recombination rate between the ruby (rb) and jet (j) loci in Tribolium castaneum. Recombination rates in the two replicate high lines increased from an average of 0.22 in the base populations to an average of 0.42 at generation 15. Recombination rate pooled over the 15 generations of selection in each low line was significantly less than the control but there was no clear downward trend in response to selection for decreased recombination rate. The realized heritabilities were 0.16 +/- 0.03 and 0.17 +/- 0.02 in the two high lines, and were not significantly different from zero in the two low lines. Selection was based on crossing over in cis females only; however, rates measured in cis males after 12 generations showed the same response patterns as female rates. Similar response patterns were also determined for recombination measured in trans males and females at generation 18 following three generations of relaxed selection. The distribution of recombination rates measured in backcross beetles [(H X L) X H and (H X L) X L] at generation 12 indicated polygenic control with those genes decreasing recombination rate being dominant. Detailed analysis of recombination rates in F1's produced by interline crosses at generation 15 confirmed the directional dominance findings. Under a polygenic model of recombination modifiers in which low recombination is dominant to high, average recombination rates will increase as inbreeding progresses, thus providing a mechanism for the production of new gene combinations in small populations.     

Genetic analysis of extended life span in Drosophila melanogaster II. Replication of the backcross test and molecular characterization of the N14 Locus

We are interested in localizing chromosomal regions that extend life span in Drosophila. Using stocks artificially selected for long life by Luckinbill and his colleagues, we have identified marker loci that are highly divergent in allelic frequencies between replicated long-lived lines and controls (Curtsinger et al., 1998). Several of the most divergent loci have been found to be associated with effects on life span in segregating backcross populations. Here we report an independent replication of the backcross test for the N14 marker locus, previously reported to extend male life spans by 12 days. The life span effect successfully replicates in males. N14 accounts for 30% of the total selection response in males. Life span extension occurs by a decrease in age-specific mortality rates at all ages, and is not attributable to modification of the slope of the age-specific mortality curve. The effect in females is small or nonexistent. Sequencing of the N14 locus shows that it is non-coding and not obviously regulatory, suggesting that the phenotypic effect arises from linkage disequilibrium with another locus or loci that directly affect life span. N14 DNA hybridizes to 63F/64A on the left arm of chromosome 3. The location is consistent with previous whole-chromosome substitution studies, and suggests new candidate genes for life span extension in Drosophila, including ras2. This revised version was published online in July 2006 with corrections to the Cover Date.     

Inheritance of resistance to the Cry1Ab Bacillus thuringiensis toxin in Ostrinia nubilalis (Lepidoptera: Crambidae)

Laboratory selection with Cry1Ab, the predominant Bacillus thuringiensis (Bt) toxin in transgenic corn, Zea mays L., produced >1000-fold resistance in two laboratory strains of European corn borer, Ostrinia nubilalis (Hübner). We tested the offspring of various crosses to determine the mode of inheritance of resistance to Cry1Ab. Patterns of inheritance of resistance were similar in the two resistant strains. The progeny of reciprocal F1 crosses (resistant male x susceptible female and vice versa) responded alike in bioassays, indicating autosomal inheritance. The median lethal concentrations (LC50 values) of F1 were intermediate between the resistant and susceptible parents, indicating approximately additive inheritance. However, the dominance of resistance increased as the concentration of Cry1Ab decreased. Analysis of progeny from backcrosses (F1 x susceptible strain) suggests that resistance was controlled by more than one locus. In particular, the fit of observed to expected mortality improved as the number of putative loci increased from 1 to 10. The polygenic nature of resistance in these two laboratory strains suggests that major genes for resistance to Cry1Ab were not common in the founding populations of O. nubilalis. A low initial frequency of major genes for Cry1Ab resistance might be an important factor in delaying evolution of resistance to Bt corn in this pest.     

Genetic and biochemical approach for characterization of resistance to Bacillus thuringiensis toxin Cry1Ac in a field population of the diamondback moth, Plutella xylostella

Four subpopulations of a Plutella xylostella (L.) strain from Malaysia (F(4) to F(8)) were selected with Bacillus thuringiensis subsp. kurstaki HD-1, Bacillus thuringiensis subsp. aizawai, Cry1Ab, and Cry1Ac, respectively, while a fifth subpopulation was left as unselected (UNSEL-MEL). Bioassays at F(9) found that selection with Cry1Ac, Cry1Ab, B. thuringiensis subsp. kurstaki, and B. thuringiensis subsp. aizawai gave resistance ratios of >95, 10, 7, and 3, respectively, compared with UNSEL-MEL (>10,500, 500, >100, and 26, respectively, compared with a susceptible population, ROTH). Resistance to Cry1Ac, Cry1Ab, B. thuringiensis subsp. kurstaki, and B. thuringiensis subsp. aizawai in UNSEL-MEL declined significantly by F(9). The Cry1Ac-selected population showed very little cross-resistance to Cry1Ab, B. thuringiensis subsp. kurstaki, and B. thuringiensis subsp. aizawai (5-, 1-, and 4-fold compared with UNSEL-MEL), whereas the Cry1Ab-, B. thuringiensis subsp. kurstaki-, and B. thuringiensis subsp. aizawai-selected populations showed high cross-resistance to Cry1Ac (60-, 100-, and 70-fold). The Cry1Ac-selected population was reselected (F(9) to F(13)) to give a resistance ratio of >2,400 compared with UNSEL-MEL. Binding studies with (125)I-labeled Cry1Ab and Cry1Ac revealed complete lack of binding to brush border membrane vesicles prepared from Cry1Ac-selected larvae (F(15)). Binding was also reduced, although less drastically, in the revertant population, which indicates that a modification in the common binding site of these two toxins was involved in the resistance mechanism in the original population. Reciprocal genetic crosses between Cry1Ac-reselected and ROTH insects indicated that resistance was autosomal and showed incomplete dominance. At the highest dose of Cry1Ac tested, resistance was recessive while at the lowest dose it was almost completely dominant. The F(2) progeny from a backcross of F(1) progeny with ROTH was tested with a concentration of Cry1Ac which would kill 100% of ROTH moths. Eight of the 12 families tested had 60 to 90% mortality, which indicated that more than one allele on separate loci was responsible for resistance to Cry1Ac.     

Genetic analysis of extended life span in Drosophila melanogaster I. RAPD screen for genetic divergence between selected and control lines

Using lines selected for long life by Luckinbill and his co-workers, we screened two selected and two control lines for allelic frequency differences at 1200 randomly chosen RAPD marker loci. Twenty-three marker loci showed frequency differences in excess of 80%, and five were greater than 90%. Age-specific effects of the five most differentiated loci were estimated by collecting complete survival data in segregating backcross populations. Alleles at four of the five marker loci were associated with significant extension of life span in males, while two marker loci had significant effects in females. Eighty percent of the total selection response in males can be explained by the identified QTL's, under the assumption of additivity. The N14+ marker allele accounted for a 12-day life span extension in males, but had little effect in females. Both sex-limited and sex-shared effects were observed. Analysis of age-specific mortality rates suggests that life span extension occurs by a combination of genetic factors that moderate both the level of mortality and the rate at which mortality increases with age. This revised version was published online in July 2006 with corrections to the Cover Date.     

Interaction ofH-2 and nonH-2 linked genes in the regulation of antibody response to a threshold dose of sheep erythrocytes

The antibody response to a threshold dose (¹⁰) of SE was studied in the High responder line (H) and the Low responder line (L) of mice obtained by bidirectional selective breeding for the character quantitative agglutinin response to an optimal dose of SE, and in interline hybrids: F₁, F₁ and both backcrosses. Whereas the interline difference in agglutinin responses at the optimal dose is due to the additive effect of about ten independently segregating loci, one of which isH-2 linked, the responsiveness to the threshold dose is determined by the effect of two loci. The direction of the dominance effect also varies with the antigen dose: high responsiveness is partially dominant at the optimal dose while at the threshold dose nonresponder character is partially dominant. The role of theH-2 linked locus was investigated. It has been demonstrated that on an identical background (equivalent to that of F₁ hybrids) this locus is responsible for 12% of the interline difference at the optimal antigen dose, and for 61% at the threshold antigen dose. For the two antigen doses, the quantitative effect of theH-2 locus is in agreement with the estimate of the number of loci obtained by variance analysis. The intervention of a second gene, non-H-2 linked, in the regulation of responsiveness to 10⁶ SE is demonstrated by appropriate assortative matings. The interaction between the two genes is discussed.     

Dominant low responsiveness in the IgG response of mice to the complex protein antigen type 1 fimbriae from Actinomyces viscosus T14V.

Type 1 fimbriae from Actinomyces viscosus T14V, composed of a complex protein of Mr 65,000, mediate the adherence of A. viscosus T14V to the host, whereas type 1 fimbriae-specific antibodies inhibit adherence. Genetic control of the serum IgG response to type 1 fimbriae was evaluated in a series of inbred, hybrid, recombinant inbred, and back-cross mice. Mice were given i.p. injections of 10(8) A. viscosus T14V cells in saline on days 0 and 14, and IgG anti-type 1 fimbriae in sera obtained on day 26 were measured by ELISA. Segregation analysis of the responses of (BALB/cJ x A/J)F1 x A/J backcross mice suggested polygenic control. Linkage analysis in (BALB/cJ x A/J)F1 x A/J backcross and SWXL recombinant inbred mice suggested control by genes linked with H-2 and with Ly-17 and Akp-1. In several F1 hybrid strains derived from H-2-disparate high and low responder parental strains, low responsiveness was dominant. The F1 derived from the H-2-identical high and low responder strains CBA/J and C3H/HeJ was a low responder, suggesting that dominant low responsiveness was mediated by non-H-2-linked genes. A three-gene model is proposed for regulation of the type 1 fimbriae response, including an MHC-linked gene, a gene linked with Ly-17 and Akp-1 on the telomeric portion of chromosome 1, and a background gene whose linkage is unknown.     

Genetics of life span in mice

Thymic involution that occurs earlier in some individuals than others may be the result of complex interactions between genetic factors and the environment. Such interactions may produce defects of thymus-dependent immune regulation associated with susceptibility to developing autoimmune diseases, malignancy, and an increased number of infections associated with aging.The major histocompatibility complex may be important in determining profiles of cause of death and length of life in mice. Genetic influences on life span involve interactions between loci and allelic interactions during life which may change following viral infections or exposure to other environmental factors. We have used different experimental protocols to study the influence of H-2 on life span and found that interactions between genetic regions, are inconsistent, particularly when comparing mice infected or not infected with Sendai virus.Genes important for life span need to be studied against many genetic backgrounds and under differing environmental conditions because of the complexity of the genetics of life span. Several genetic models were used to demonstrate that the MHC is a marker of life span in backcross and intercross male mice of the H-2^d and H-2^b genotypes in B10 congenic mice. Females lived longer than males in backcross and intercross mice, while males lived longer than females in B10 congenics. H-2^d was at a disadvantage for life span in backcross mice of the dilute brown and brown males exposed to Sendai infection, but intercross mice not exposed to Sendai virus of the same genotype were not at a disadvantage. H-2^d mice were not disadvantaged when compared to H-2^b in B10 congenics that had not been exposed to Sendai virus infection but the reverse was true when they were exposed. Overall, all our studies suggest that genetic influences in life span may involve interactions between loci and many allelic interactions in growing animals or humans. These genetic influences on life span may vary after they are exposed to infections or other environmental conditions. This paper emphasizes the need to use several genetic models, especially animals that have been monitored for infections, to study the genetics of life span.     

Genetic control of lipopolysaccharide induced generation of serum colony stimulating factor and proliferation of splenic granulocyte/macrophage precursor cells.

Administration of bacterial lipopolysaccharides (LPS) to mice causes a rise in tissue and serum colony stimulating factor (CSF) levels and in bone marrow and splenic colony forming cells (CFC). Two inbred strains of mice differing in their response to LPS were used to study the genetic control of LPS induced granulopoietic responses: a high responder strain (C3H/eB) which reacts to LPS by an elevation in serum CSF and by an increase in splenic CFC levels, and a low responder strain (C3H/HeJ) which fails to show these responses. The ability to generate serum CSF after administration of LPS is controlled by a single autosomal dominant gene, while the splenic CFC response to LPS follows the characteristic patterns of a polygenic inheritance control. The associated relationships of CSF and CFC responsiveness have been investigated in backcross (F1 X C3H/Hej) and F2 mice. Most mice which generated high levels of CSF showed a high or intermediate CFC response and most mice which did not generate any detectable levels of serum CSF showed a low splenic CFC response. The results suggest that CSF may play a physiologic role in vivo as a granulopoietin. In addition it was shown that the genetic control mechanisms governing the CSF/CFC responses are determined by the lipid A-KDO portion of the LPS molecule, suggesting that lipid A is the active part of the LPS molecule in stimulating granulopoiesis.     

Chromosomal basis of viability differences in Tigriopus californicus interpopulation hybrids

Crosses between populations of Tigriopus californicus result in backcross and F2 hybrid breakdown for a variety of fitness related measures. The magnitude of this hybrid breakdown is correlated with evolutionary divergence. We assessed the chromosomal basis of viability differences in nonrecombinant backcross hybrids using markers mapped to individual chromosomes. To assess effects of evolutionary divergence we crossed one population to three different populations: two distantly related (approximately 18% mitochondrial COI sequence divergence) and one closely related (approximately 1% mitochondrial COI sequence divergence). We found that all three interpopulation crosses resulted in significant deviations from expected Mendelian ratios at a majority of the loci studied. In all but one case, deviations were due to a deficit of parental homozygotes. This pattern implies that populations of T. californicus carry a significant genetic load, and that a combination of beneficial dominance and deleterious homozygote-heterozygote interactions significantly affects hybrid viability. Pairwise tests of linkage disequilibrium detected relatively few significant interactions. For the two divergent crosses, effects of individual chromosomes were highly concordant. These two crosses also showed higher heterozygote excess in females than males across the vast majority of chromosomes.     

Additive and over-dominant effects resulting from epistatic loci are the primary genetic basis of heterosis in rice

A set of 148 F9 recombinant inbred lines (RILs) was developed from the cross of an indica cultivar 93-11 and japonica cultivar DTT13,showing strong F1 heterosis.Subsequently,two backcross F1 (BCF1) populations were constructed by backcrossing these 148 RILs to two parents,93-11 and DT713.These three related populations (281BCF1 lines,148 RILs) were phenotyped for six yield-related traits in two locations.Significant inbreeding depression was detected in the population of RILS and a high level of heterosis was observed in the two BCF1 populations.A total of 42 main-effect quantitative trait loci (M-QTLs) and 109 epistatic effect QTL pairs (E-QTLs) were detected in the three related populations using the mixed model approach.By comparing the genetic effects of these QTLs detected in the RILs,BCF1 performance and mid-parental heterosis (HMp),we found that,in both BCF1 populations,the QTLs detected could be classified into two predominant types:additive and over-domlnant loci,which indicated that the additive and over-dominant effect were more important than complete or partially dominance for M-QTLs and E-QTLs.Further,we found that the E-QTLs detected collectively explained a larger portion of the total phenotypic variation than the M-QTLs in both RILs and BCF1 populations.All of these results suggest that additive and over-dominance resulting from epistatic loci might be the primary genetic basis of heterosis in rice.     

Salmonella typhimurium infection in high and low antibody responder mice: inverse correlation between antibody responsiveness and resistance to infection

Susceptibility to Salmonella typhimurium infection was compared in H (high Ab responder) and L (low Ab responder) mice obtained by several selective breeding experiments (Selections I, II, III, IV and IV A). H mice were always much more susceptible to infection than their L mice counterparts within a continuous LD 50 variation range. In three of the selections (I, II and IV A) the low responsiveness character is known to result mainly from rapid Ag degradation in L mice macrophages. It was hypothesized that resistance to multiplication of intracellular pathogens could be related to an increased catabolic activity towards Ag. This was actually demonstrated, in F2 segregant hybrids of selection IV A, by the significant inverse correlation between capacity for Ab production and resistance to infection.     

Determining the mode of inheritance of pesticide resistance with backcross experiments

The most widely used method for evaluating the mode of inheritance of pesticide resistance is based on bioassays of individuals from a backcross between F1 (hybrid of resistant and susceptible strains) and parental resistant or susceptible strains. Monte Carlo simulations of the standard backcross method showed that the probability of incorrectly rejecting the null hypothesis of monogenic inheritance (Type I error) was generally more than double the conventional value of P = 0.05. Conversely, the null hypothesis of monogenic inheritance was likely to be accepted in a relatively large proportion of cases in which resistance is controlled by two or more loci. Expected differences in mortality of backcross offspring between monogenic and additive polygenic models approached zero as dose approached extremely low values, extremely high values, and the LD50 of the backcross generation. Thus, the effectiveness of the backcross method depended strongly on dose. The power of the standard backcross method to correctly reject the null hypothesis of monogenic inheritance increased as number of loci, slope of parental dose-mortality lines, magnitude of resistance, and sample size increased. Guidelines for improving the design and interpretation of backcross experiments are presented.     

The quantitative genetics of fluctuating asymmetry

Fluctuating asymmetry (subtle departures from identical expression of a trait across an axis of symmetry) in many taxa is under stabilizing selection for reduced asymmetry. However, lack of reliable estimates of genetic parameters for asymmetry variation hampers our ability to predict the evolutionary outcome of this selection. Here we report on a study, based on analysis of variation within and between isofemale lines and of generation means (line-cross analysis), designed to dissect in detail the quantitative genetics of positional fluctuating asymmetry (PFA) in bristle number in natural populations of Drosophila falleni. PFA is defined as the difference between the two sides of the body in the placement or position of components of a meristic trait. Heritability (measured at 25 degrees C) of two related measures of PFA were 13% and 21%, both of which differed significantly from zero. In contrast, heritability estimates for fluctuating asymmetry in the total number of anterior (0.7%) and transverse (2.4%) sternopleural bristles were smaller, not significant, and in quantitative agreement with previously published estimates. Heritabilities for bristle number (trait size) were considerably greater than that for any asymmetry measure. The experimental design controlled for the potentially confounding effects of common familial environment, and repeated testing revealed that PFA differences between lines were genetically stable for up to 16 generations in the laboratory at 25 degrees C. We performed line cross analysis between strains at the extremes of the PFA distribution (highest and lowest values); parental strains, F1, F1r (reciprocal), F2, backcross, and backcross reciprocal generations were represented. The inheritance of PFA was described best by additive and dominance effects localized to the X-chromosomes, whereas autosomal dominance effects were also detected. Epistatic, maternal, and cytoplasmic effects were not detected. The inheritance of trait size was notably more complex and involved significant autosomal additive, dominance, and epistatic effects; maternal dominance effects; and additive and dominance effects localized to the X-chromosomes. The additive genetic correlation between PFA and its associated measure of trait size was negative (-0.049), but not statistically significant, indicating that the loci contributing additive genetic effects to these traits are probably different. It is suggested that PFA may be a sensitive measure of developmental instability because PFA taps the ability of an organism to integrate interconnected developmental pathways.     

Microsatellite markers associated with quantitative trait loci controlling antibody response to Escherichia coli and Salmonella enteritidis in young broilers

A unique resource population was produced to facilitate detection of microsatellite markers associated with quantitative trait loci controlling antibody (Ab) response in broiler chickens. Three F1 males were produced by mating two lines divergently selected on Ab response to Escherichia coli vaccination. Each F1 male was mated with females from four genetic backgrounds: F1, high-Ab line (HH), low-Ab line and commercial line, producing three resource families, each with four progeny types. About 1700 chicks were immunized with E. coli and Salmonella enteritidis vaccines. Selective genotyping was conducted on the individuals with highest or lowest average Ab to E. coli and S. enteritidis within each progeny type in each sire family. Twelve markers were significantly associated with Ab to E. coli and six of them were also associated with Ab to S. enteritidis, mostly exhibiting a similar low effect (approximately 0.35 phenotypic SD) in all progeny types. Four markers exhibited a highly significant and much larger effect (approximately 1.7 SD), but only in progeny of females from the HH, suggesting that a backcross to the high parental line should be preferred over the commonly used F2 population. Results from two markers suggested a quantitative trait locus on chromosome 2 around 400 cM. The marker MCW0083, significant in two sire families, is closely linked to the bone morphogenetic protein 2 (BMP2) gene, known to be associated with the control of T-cell transformation in humans.     

Sugarcane borer (Lepidoptera: Crambidae) resistance to transgenic Bacillus thuringiensis maize

Transgenic maize, Zea mays L., expressing the Bacillus thuringiensis (Bt) CrylAb toxin has been planted to extensive areas across the United States and several other countries, but no resistance has been documented in field populations oflepidopteran target pests. This article describes the first report of resistance alleles to commercially available Cry1Ab Bt maize in a Louisiana population of sugarcane borer, Diatraea saccharalis (F.) (Lepidoptera: Crambidae). Two hundred thirteen two-parent isolines of D. saccharalis were screened for Cry1Ab resistance on Bt maize leaf tissue using an F2 screening technique. Larvae representing three isolines survived >15 d on Bt tissue in the F2 generation. The second generation backcross progeny (B1F2) derived from isoline 52 completed larval development on Bt maize in the greenhouse. Segregation and resistance frequency analysis associated with isoline 52 suggested that Bt resistance is probably determined by a nearly completely recessive allele at a single locus. With this assumption, the estimated resistance allele frequency in this population is 0.0023, within a 95% confidence interval of 0.0003-0.0064.     

Susceptibility of bollworm and tobacco budworm (Lepidoptera: Noctuidae) to Cry2Ab2 insecticidal protein

Susceptibilities of 82 bollworm, Helicoverpa zea (Boddie), and 44 tobacco budworm, Heliothis virescens (F.) (Lepidoptera: Noctuidae), populations to Cry2Ab2 protein were measured in diet incorporated assays at the University of Arkansas from 2002 to 2005. Resulting data were used to calculate overall (pooled data) estimates of species susceptibility for future benchmarks of resistance. Variabilities among populations also were studied by comparing regressions for individual populations and calculating mean susceptibilities for different subgroups of the colonies studied. Individual lethal concentration (LC50) estimates for nine laboratory, seven laboratory-cross, and 28 field populations of H. virescens varied up to 48-fold when adjusted for the response of the most susceptible laboratory colony studied. Mean susceptibilities of all laboratory, laboratory-cross, or field colonies varied only two-fold. When grouped by host plants, populations collected on tobacco, Nicotiana tabacum (L.), seemed to be less susceptible than those collected on other host plants. Individual LC50 values for 82 laboratory, laboratory-cross and field populations of H. zea varied up to 37-fold. Mean LC50 values of all laboratory, laboratory-cross, or field populations varied only three-fold. Susceptibilities of populations from Bollgard cotton were up to four-fold less than those from Bacillus thuringiensis corn, Zea mays L. Field populations collected during late season were generally less susceptible than those collected early in the season. Across the two species, H. zea was less sensitive to Cry2Ab2 than H. virescens. Both species seem to be less sensitive to Cry2Ab2 than to CrylAc.     

A nearest-neighboring-end algorithm for genetic mapping

MOTIVATION: High-throughput methods are beginning to make possible the genotyping of thousands of loci in thousands of individuals, which could be useful for tightly associating phenotypes to candidate loci. Current mapping algorithms cannot handle so many data without building hierarchies of framework maps. RESULTS: A version of Kruskal's minimum spanning tree algorithm can solve any genetic mapping problem that can be stated as marker deletion from a set of linkage groups. These include backcross, recombinant inbred, haploid and double-cross recombinational populations, in addition to conventional deletion and radiation hybrid populations. The algorithm progressively joins linkage groups at increasing recombination fractions between terminal markers, and attempts to recognize and correct erroneous joins at peaks in recombination fraction. The algorithm is O (mn3) for m individuals and n markers, but the mean run time scales close to mn2. It is amenable to parallel processing and has recovered true map order in simulations of large backcross, recombinant inbred and deletion populations with up to 37,005 markers. Simulations were used to investigate map accuracy in response to population size, allelic dominance, segregation distortion, missing data and random typing errors. It produced accurate maps when marker distribution was sufficiently uniform, although segregation distortion could induce translocated marker orders. The algorithm was also used to map 1003 loci in the F7 ITMI population of bread wheat, Triticum aestivum L. emend Thell., where it shortened an existing standard map by 16%, but it failed to associate blocks of markers properly across gaps within linkage groups. This was because it depends upon the rankings of recombination fractions at individual markers, and is susceptible to sampling error, typing error and joint selection involving the terminal markers of nearly finished linkage groups. Therefore, the current form of the algorithm is useful mainly to improve local marker ordering in linkage groups obtained in other ways. AVAILABILITY: The source code and supplemental data are http://www.iubio.bio.indiana.edu/soft/molbio/qtl/flipper/ CONTACT: ccrane@purdue.edu.