The isolation, characterization, and application of WRKY genes as useful molecular markers in tropical trees

The improvement of tropical tree crops using conventional breeding methods faces challenges due to the length of time involved. Thus, like most crops, there is an effort to utilize molecular genetic markers in breeding programs to select for desirable agronomic traits. Known as marker assisted breeding or marker assisted selection, genetic markers associated with a phenotype of interest are used to screen and select material reducing the time necessary to evaluate candidates. As the focus of this research was improving disease resistance in tropical trees, the usefulness of the WRKY gene superfamily was investigated as candidates for generating useful molecular genetic markers. WRKY genes encode plant-specific transcriptional factors associated with regulating plants' responses to both biotic and abiotic stress. ^ One pair of degenerate primers amplified 48 WRKY gene fragments from three taxonomically distinct, economically important, tropical tree crop species: 18 from Theobroma cacao L., 21 from Cocos nucifera L. and 9 from Persea americana Mill. Several loci from each species were polymorphic because of single nucleotide substitutions present within a putative non-coding region of the loci. Capillary array electrophoresis-single strand conformational polymorphism (CAE-SSCP) mapped four WRKY loci onto a genetic linkage map of a T. cacao F2 population segregating for resistance to witches' broom disease. Additionally, PCR primers specific for four T. cacao loci successfully amplified WRKY loci from 15 members of the Byttneriae tribe. A method was devised to allow the reliable discrimination of alleles by CAE-SSCP using only the mobility assigned to the sample peaks. Once this method was validated, the diversity of three WRKY loci was evaluated in a germplasm collection of T. cacao . One locus displayed high diversity in the collection, with at least 18 alleles detected from mobility differences of the product peaks. The number of WRKY loci available within the genome, ease of isolation by degenerate PCR, codominant segregation demonstrated in the F2 population, and usefulness for screening germplasm collections and closely related wild species demonstrates that the WRKY superfamily of genes are excellent candidates for developing a number of genetic molecular markers for breeding purposes in tropical trees. ^

Investigating the origin of PM2.5 in Baltimore using highly time-resolved organic molecular markers measured during the Baltimore PM-supersite

The detailed organic composition of atmospheric fine particles with an aerodynamic diameter smaller than or equal to 2.5 micrometers (PM2.5) is an integral part of the knowledge needed in order to fully characterize its sources and transformation in the environment. For the study presented here, samples were collected at 3-hour intervals. This high time resolution allows gaining unique insights on the influence of short- and long-range transport phenomena, and dynamic atmospheric processes. A specially designed sequential sampler was deployed at the 2002-2003 Baltimore PM-Supersite to collect PM2.5 samples at a 3-hourly resolution for extended periods of consecutive days, during both summer and winter seasons. Established solvent-extraction and GC-MS techniques were used to extract and analyze the organic compounds in 119 samples from each season. Over 100 individual compounds were quantified in each sample. For primary organics, averaging the diurnal ambient concentrations over the sampled periods revealed ambient patterns that relate to diurnal emission patterns of major source classes. Several short-term releases of pollutants from local sources were detected, and local meteorological data was used to pinpoint possible source regions. Biogenic secondary organic compounds were detected as well, and possible mechanisms of formation were evaluated. The relationships between the observed continuous variations of the concentrations of selected organic markers and both the on-site meteorological measurements conducted parallel to the PM2.5 sampling, and the synoptic patterns of weather and wind conditions were also examined. Several one-to-two days episodes were identified from the sequential variation of the concentration observed for specific marker compounds and markers ratios. The influence of the meteorological events on the concentrations of the organic compounds during selected episodes was discussed. It was observed that during the summer, under conditions of pervasive influence of air masses originated from the west/northwest, some organic species displayed characteristics consistent with the measured PM2.5 being strongly influenced by the aged nature of these long-traveling background parcels. During the winter, intrusions from more regional air masses originating from the south and the southwest were more important.

Using soil profiles of seeds and molecular markers as proxies for sawgrass and wet prairie slough vegetation in Shark Slough, Everglades National Park

We measured the abundance of Cladium jamaicense (Crantz) seeds and three biomarkers in freshwater marsh soils in Shark River Slough (SRS), Everglades National Park (ENP) to determine the degree to which these paleoecological proxies reflect spatial and temporal variation in vegetation. We found that C. jamaicense seeds and the biomarkers Paq, total lignin phenols (TLP) and kaurenes analyzed from surface soils were all significantly correlated with extant aboveground C. jamaicense biomass quantified along a vegetation gradient from a C. jamaicense to a wet prairie/slough (WPS) community. Our results also suggest that these individual proxies may reflect vegetation over different spatial scales: Paq and kaurenes correlated most strongly (R 2 = 0.88 and 0.99, respectively) with vegetation within 1 m of a soil sample, while seeds and TLP reflected vegetation 0–20 m upstream of soil samples. These differences in the spatial scale depicted by the different proxies may be complementary in understanding aspects of historic landscape patterning. Soil profiles of short (25 cm) cores showed that downcore variation in C. jamaicense seeds was highly correlated with two of the three biomarkers (Paq, R 2 = 0.84, p<0.005; TLP, R 2 = 0.97, p<0.0001), and all four of the proxies indicated a recent increase in C. jamaicense biomass at the site. Using a preliminary depth-to-age relationship based on matching charcoal peaks with available ENP fire records (1980-present) specific to our coring site, we found that peak-depths in C. jamaicense seed concentration appeared to correspond to recent minimum water levels (e.g., 1989 and 2001), and low seed abundance corresponded to high water levels (e.g., 1995), consistent with the known autecology of C. jamaicense. In summary, the combination of C. jamaicense seeds and biomarkers may be useful for paleoecological reconstruction of vegetation change and ultimately in guaging the success of ongoing efforts to restore historic hydrologic conditions in the South Florida Everglades.

Investigating the origin of PM2.5 in Baltimore using highly time-resolved organic molecular markers measured during the Baltimore PM-supersite

The detailed organic composition of atmospheric fine particles with an aerodynamic diameter smaller than or equal to 2.5 micrometers (PM 2.5) is an integral part of the knowledge needed in order to fully characterize its sources and transformation in the environment. For the study presented here, samples were collected at 3-hour intervals. This high time resolution allows gaining unique insights on the influence of short- and long-range transport phenomena, and dynamic atmospheric processes. A specially designed sequential sampler was deployed at the 2002-2003 Baltimore PM Supersite to collect PM2.5 samples at a 3-hourly resolution for extended periods of consecutive days, during both summer and winter seasons. Established solvent-extraction and GC-MS techniques were used to extract and analyze the organic compounds in 119 samples from each season. Over 100 individual compounds were quantified in each sample. For primary organics, averaging the diurnal ambient concentrations over the sampled periods revealed ambient patterns that relate to diurnal emission patterns of major source classes. Several short-term releases of pollutants from local sources were detected, and local meteorological data was used to pinpoint possible source regions. Biogenic secondary organic compounds were detected as well, and possible mechanisms of formation were evaluated. The relationships between the observed continuous variations of the concentrations of selected organic markers and both the on-site meteorological measurements conducted parallel to the PM2.5 sampling, and the synoptic patterns of weather and wind conditions were also examined. Several one-to-two days episodes were identified from the sequential variation of the concentration observed for specific marker compounds and markers ratios. The influence of the meteorological events on the concentrations of the organic compounds during selected episodes was discussed. It was observed that during the summer, under conditions of pervasive influence of air masses originated from the west/northwest, some organic species displayed characteristics consistent with the measured PM2.5 being strongly influenced by the aged nature of these long-traveling background parcels. During the winter, intrusions from more regional air masses originating from the south and the southwest were more important.

Organic matter dynamics in the Florida coastal Everglades: A molecular marker and isotopic approach

Everglades National Park (ENP) is about to undergo the world's largest wetland restoration with the aim of improving the quality, timing and distribution of water flow. The changes in water flow are hypothesized to alter the nutrient fluxes and organic matter (OM) dynamics within ENP, especially in the estuarine areas. This study used a multi-proxy approach of molecular markers and stable δ 13C isotope measurements, to determine the present day OM dynamics in ENP. ^ OM dynamics in wetland soils/sediments have proved to be difficult to understand using traditional geochemical approaches. These are often inadequate to describe the multitude of OM sources (e.g. higher land plant, emergent vegetation, submerged vegetation) to the soils/sediments and the complex diagenetic processes that can alter the OM characteristics. A multi-proxy approach, however, that incorporates both molecular level and bulk parameter information is ideal to comprehend complex OM dynamics in aquatic environments. Therefore, biomass-specific molecular markers or proxies can be useful in tracing the sources and processing of OM. This approach was used to examine the OM dynamics in the two major drainage basins, Shark River Slough and Taylor River Slough, of ENP. Freshwater to marine transects were sampled in both systems for soils/sediments and suspended particulate organic matter (SPOM) to be characterized through bulk OM analyses, lipid biomarker determinations (e.g. sterols, fatty acids, hydrocarbons and triterpenoids) and compound-specific stable carbon isotope (δ 13C) determinations. ^ One key accomplishment of the research was the assessment of a molecular marker proxy (Paq) to distinguish between emergent/higher plant vegetation from submerged vegetation within ENP. This proxy proved to be quite useful at tracing OM inputs to the soils/sediments of ENP. A second key accomplishment was the development of a 3-way model using vegetation specific molecular markers. This novel, descriptive model was successfully applied to the estuarine areas of Taylor and Shark River sloughs, providing clear evidence of mixing of freshwater, estuarine and marine derived OM in these areas. In addition, diagenetic transformations of OM in these estuaries were found to be quite different between Taylor and Shark Rivers, and are likely a result of OM quality and hydrological differences. ^

Molecular markers and human history: A tale of two haplogroup systems

To chronicle demographic movement across African Asian corridors, a variety of molecular (sequence analysis, restriction mapping and denaturing high performance liquid chromatography etc.) and statistical (correspondence analysis, AMOVA, calculation of diversity indices and phylogenetic inference, etc.) techniques were employed to assess the phylogeographic patterns of mtDNA control region and Y chromosomal variation among 14 sub-Saharan, North African and Middle Eastern populations. The patterns of genetic diversity revealed evidence of multiple migrations across several African Asian passageways as well within the African continent itself. The two-part analysis uncovered several interesting results which include the following: (1) a north (Egypt and Middle East Asia) to south (sub-Saharan Africa) partitioning of both mtDNA and Y chromosomal haplogroup diversity, (2) a genetic diversity gradient in sub-Saharan Africa from east to west, (3) evidence in favor of the Levantine Corridor over the Horn of Africa as the major genetic conduit since the Last Glacial Maximum, (4) a substantially higher mtDNA versus Y chromosomal sub-Saharan component in the Middle East collections, (5) a higher representation of East versus West African mtDNA haplotypes in the Arabian Peninsula populations versus no such bias in the Levant groups and lastly, (6) genetic remnants of the Bantu demographic expansion in sub-Saharan Africa. ^

A Case-Only Genome-wide Association Study of Gender- and Age-specific Risk Markers for Childhood Leukemia

Males and age group 1 to 5 years show a much higher risk for childhood acute lymphoblastic leukemia (ALL). We performed a case-only genome-wide association study (GWAS), using the Illumina Infinium HumanCoreExome Chip, to unmask gender- and age-specific risk variants in 240 non-Hispanic white children with ALL recruited at Texas Children’s Cancer Center, Houston, Texas. Besides statistically most significant results, we also considered results that yielded the highest effect sizes. Existing experimental data and bioinformatic predictions were used to complement results, and to examine the biological significance of statistical results. Our study identified novel risk variants for childhood ALL. The SNP, rs4813720 (RASSF2), showed the statistically most significant gender-specific associations (P < 2 x 10-6). Likewise, rs10505918 (SOX5) yielded the lowest P value (P < 1 x 10-5) for age-specific associations, and also showed the statistically most significant association with age-at-onset (P < 1 x 10-4). Two SNPs, rs12722042 and 12722039, from the HLA-DQA1 region yielded the highest effect sizes (odds ratio (OR) = 15.7; P = 0.002) for gender-specific results, and the SNP, rs17109582 (OR = 12.5; P = 0.006), showed the highest effect size for age-specific results. Sex chromosome variants did not appear to be involved in gender-specific associations. The HLA-DQA1 SNPs belong to DQA1*01:07and confirmed previously reported male-specific association with DQA1*01:07. Twenty one of the SNPs identified as risk markers for gender- or age-specific associations were located in the transcription factor binding sites and 56 SNPs were non-synonymous variants, likely to alter protein function. Although bioinformatic analysis did not implicate a particular mechanism for gender- and age-specific associations, RASSF2 has an estrogen receptor-alpha binding site in its promoter. The unknown mechanisms may be due to lack of interest in gender- and age-specificity in associations. These results provide a foundation for further studies to examine the gender- and age-differential in childhood ALL risk. Following replication and mechanistic studies, risk factors for one gender or age group may have a potential to be used as biomarkers for targeted intervention for prevention and maybe also for treatment.

A case-only genome-wide association study of gender- and age-specific risk markers for childhood leukemia

Males and age group 1 to 5 years show a much higher risk for childhood acute lymphoblastic leukemia (ALL). We performed a case-only genome-wide association study (GWAS), using the Illumina Infinium HumanCoreExome Chip, to unmask gender- and age-specific risk variants in 240 non-Hispanic white children with ALL recruited at Texas Children’s Cancer Center, Houston, Texas. Besides statistically most significant results, we also considered results that yielded the highest effect sizes. Existing experimental data and bioinformatic predictions were used to complement results, and to examine the biological significance of statistical results. ^ Our study identified novel risk variants for childhood ALL. The SNP, rs4813720 (RASSF2), showed the statistically most significant gender-specific associations (P < 2 x 10-6). Likewise, rs10505918 (SOX5) yielded the lowest P value (P < 1 x 10-5 ) for age-specific associations, and also showed the statistically most significant association with age-at-onset (P < 1 x 10-4). Two SNPs, rs12722042 and 12722039, from the HLA-DQA1 region yielded the highest effect sizes (odds ratio (OR) = 15.7; P = 0.002) for gender-specific results, and the SNP, rs17109582 (OR = 12.5; P = 0.006), showed the highest effect size for age-specific results. Sex chromosome variants did not appear to be involved in gender-specific associations. ^ The HLA-DQA1 SNPs belong to DQA1*01:07and confirmed previously reported male-specific association with DQA1*01:07. Twenty one of the SNPs identified as risk markers for gender- or age-specific associations were located in the transcription factor binding sites and 56 SNPs were non-synonymous variants, likely to alter protein function. Although bioinformatic analysis did not implicate a particular mechanism for gender- and age-specific associations, RASSF2 has an estrogen receptor-alpha binding site in its promoter. The unknown mechanisms may be due to lack of interest in gender- and age-specificity in associations. These results provide a foundation for further studies to examine the gender- and age-differential in childhood ALL risk. Following replication and mechanistic studies, risk factors for one gender or age group may have a potential to be used as biomarkers for targeted intervention for prevention and maybe also for treatment.^

Horizontal genetic transfer in asexual fungi

Four aspects of horizontal genetic transfer during heterokaryon formation were examined in the asexual pathogen Fusarium oxysporum f.sp. cubense (Foc): (1) variability based on method of heterokaryon formation; (2) differences in nuclear and mitochondrial inheritance; (3) the occurrence of recombination without nuclear fusion; (4) the occurrence of horizontal genetic transfer between distantly related isolates. The use of non-pathogenic strains of Fusarium oxysporum as biocontrol agents warrants a closer examination at the reproductive life cycle of this fungus, particularly if drug resistance or pathogenicity genes can be transmitted horizontally. Experiments were divided into three phases. Phase I looked at heterokaryon formation by hyphal anastomosis and protoplast fusion. Phase II was a time course of heterokaryon formation to look at patterns of nuclear and mitochondrial inheritance. Phase III examined the genetic relatedness of the different vegetative compatibility groups using a multilocus analysis approach. Heterokaryon formation was evident within and between vegetative compatibility groups. Observation of non-parental genotypes after heterokaryon formation confirmed that, although a rare event, horizontal genetic transfer occurred during heterokaryon formation. Uniparental mitochondria inheritance was observed in heterokaryons formed either by hyphal anastomosis or protoplast fusion. Drug resistance was expressed during heterokaryon formation, even across greater genetic distances than those distances imposed by vegetative compatibility. Phylogenies inferred from different molecular markers were incongruent at a significant level, challenging the clonal origins of Foc. Mating type genes were identified in this asexual pathogen Polymorphisms were detected within a Vegetative Compatibility Group (VCG) suggesting non-clonal inheritance and/or sexual recombination in Foc. This research was funded in part by a NIH-NIGMS (National Institutes of Health-National Institute of General Medical Sciences) Grant through the MBRS (Minority Biomedical Research Support), the Department of Biological Sciences and the Tropical Biology Program at FIU. ^

Genetic individualization of Cannabis sativa by a short tandem repeat multiplex system

Cannabis sativa is the most frequently used of all illicit drugs in the United States. Cannabis has been used throughout history for its stems in the production of hemp fiber, for its seed for oil and food, and for its buds and leaves as a psychoactive drug. Short tandem repeats (STRs), were chosen as molecular markers because of their distinct advantages over other genetic methods. STRs are co-dominant, can be standardized such that reproducibility between laboratories can be easily achieved, have a high discrimination power and can be multiplexed. ^ In this study, six STR markers previously described for Cannabis were multiplexed into one reaction. The multiplex reaction was able to individualize 98 Cannabis samples (14 hemp and 84 marijuana, authenticated as originating from 33 of the 50 United States) and detect 29 alleles averaging 4.8 alleles per loci. The data did not relate the samples from the same state to each other. This is the first study to report a single reaction six-plex and apply it to the analysis of almost 100 Cannabis samples of known geographic collection site. ^

Fatty acids and stable isotopes as indicators of early-life feeding and potential maternal resource dependency in the bull shark Carcharhinus leucas

The degree of reliance of newborn sharks on energy reserves from maternal resource allocation and the timescales over which these animals develop foraging skills are critical factors towards understanding the ecological role of top predators in marine ecosystems. We used muscle tissue stable carbon isotopic composition and fatty acid analysis of bull sharks Carcharhinus leucas to investigate early-life feeding ecology in conjunction with maternal resource dependency. Values of δ13C of some young-of-the-year sharks were highly enriched, reflecting inputs from the marine-based diet and foraging locations of their mothers. This group of sharks also contained high levels of the 20:3ω9 fatty acid, which accumulates during periods of essential fatty acid deficiency, suggesting inadequate or undeveloped foraging skills and possible reliance on maternal provisioning. A loss of maternal signal in δ13C values occurred at a length of approximately 100 cm, with muscle tissue δ13C values reflecting a transition from more freshwater/estuarine-based diets to marine-based diets with increasing length. Similarly, fatty acids from sharks >100 cm indicated no signs of essential fatty acid deficiency, implying adequate foraging. By combining stable carbon isotopes and fatty acids, our results provided important constraints on the timing of the loss of maternal isotopic signal and the development of foraging skills in relation to shark size and imply that molecular markers such as fatty acids are useful for the determination of maternal resource dependency.

Quantitative and Qualitative Aspects of Dissolved Organic Carbon Leached from Senescent Plants in an Oligotrophic Wetland

We conducted a series of experiments whereby dissolved organic matter (DOM) was leached from various wetland and estuarine plants, namely sawgrass (Cladium jamaicense), spikerush (Eleocharis cellulosa), red mangrove (Rhizophora mangle), cattail (Typha domingensis), periphyton (dry and wet mat), and a seagrass (turtle grass; Thalassia testudinum). All are abundant in the Florida Coastal Everglades (FCE) except for cattail, but this species has a potential to proliferate in this environment. Senescent plant samples were immersed into ultrapure water with and without addition of 0.1% NaN3 (w/ and w/o NaN3, respectively) for 36 days. We replaced the water every 3 days. The amount of dissolved organic carbon (DOC), sugars, and phenols in the leachates were analyzed. The contribution of plant leachates to the ultrafiltered high molecular weight fraction of DOM (>1 kDa; UDOM) in natural waters in the FCE was also investigated. UDOM in plant leachates was obtained by tangential flow ultrafiltration and its carbon and phenolic compound compositions were analyzed using solid state 13C cross-polarization magic angle spinning nuclear magnetic resonance (13C CPMAS NMR) spectroscopy and thermochemolysis in the presence of tetramethylammonium hydroxide (TMAH thermochemolysis), respectively. The maximum yield of DOC leached from plants over the 36-day incubations ranged from 13.0 to 55.2 g C kg−1 dry weight. This amount was lower in w/o NaN3 treatments (more DOC was consumed by microbes than produced) except for periphyton. During the first 2 weeks of the 5 week incubation period, 60–85% of the total amount of DOC was leached, and exponential decay models fit the leaching rates except for periphyton w/o NaN3. Leached DOC (w/ NaN3) contained different concentrations of sugars and phenols depending on the plant types (1.09–7.22 and 0.38–12.4 g C kg−1 dry weight, respectively), and those biomolecules comprised 8–34% and 4–28% of the total DOC, respectively. This result shows that polyphenols that readily leach from senescent plants can be an important source of chromophoric DOM (CDOM) in wetland environments. The O-alkyl C was found to be the major C form (55±9%) of UDOM in plant leachates as determined by 13C CPMAS NMR. The relative abundance of alkyl C and carbonyl C was consistently lower in plant-leached UDOM than that in natural water UDOM in the FCE, which suggests that these constituents increase in relative abundance during diagenetic processing. TMAH thermochemolysis analysis revealed that the phenolic composition was different among the UDOM leached from different plants, and was expected to serve as a source indicator of UDOM in natural water. Polyphenols are, however, very reactive and photosensitive in aquatic environments, and thus may loose their plant-specific molecular characteristics shortly. Our study suggests that variations in vegetative cover across a wetland landscape will affect the quantity and quality of DOM leached into the water, and such differences in DOM characteristics may affect other biogeochemical processes.