Sporadic Genetic Connectivity among Small Insular Populations of the Rare Geoendemic Plant Caulanthus amplexicaulis var. barbarae (Santa Barbara Jewelflower)

Globally, a small number of plants have adapted to terrestrial outcroppings of serpentine geology, which are characterized by soils with low levels of essential mineral nutrients (N, P, K, Ca, Mo) and toxic levels of heavy metals (Ni, Cr, Co). Paradoxically, many of these plants are restricted to this harsh environment. Caulanthus ampexlicaulis var. barbarae (Brassicaceae) is a rare annual plant that is strictly endemic to a small set of isolated serpentine outcrops in the coastal mountains of central California. The goals of the work presented here were to 1) determine the patterns of genetic connectivity among all known populations of C. ampexlicaulis var. barbarae, and 2) estimate contemporary effective population sizes (Ne), to inform ongoing genomic analyses of the evolutionary history of this taxon, and to provide a foundation upon which to model its future evolutionary potential and long-term viability in a changing environment. Eleven populations of this taxon were sampled, and population-genetic parameters were estimated using 11 nuclear microsatellite markers. Contemporary effective population sizes were estimated using multiple methods and found to be strikingly small (typically Ne < 10). Further, our data showed that a substantial component of genetic connectivity of this taxon is not at equilibrium, and instead showed sporadic gene flow. Several lines of evidence indicate that gene flow between isolated populations is maintained through long-distance seed dispersal (e.g., >1 km), possibly via zoochory.

Development of 15 microsatellite loci in the endangered Streptanthus glandulosus subsp. niger (Brassicaceae)

PREMISE OF THE STUDY

The endangered Streptanthus glandulosus subsp. niger (Brassicaceae) is endemic to a single peninsula in California and threatened by fragmentation. We developed microsatellite markers to investigate genetic diversity in the two extant populations and the degree to which they have diverged from one another.

METHODS AND RESULTS

We used Illumina HiSeq high-throughput sequencing to develop 15 microsatellite markers, 14 of which were polymorphic. These di- and trinucleotide repeats yielded one to 11 alleles per locus in 61 plants across the two populations. Levels of observed and expected heterozygosities ranged from 0.108 to 0.946 and 0.257 to 0.839, respectively. We demonstrated cross-amplification in a second rare subspecies, S. glandulosus subsp. secundus, and in the widespread congener S. tortuosus.

CONCLUSIONS

These are the first microsatellites reported for this subspecies, and they will aid in the inclusion of genetic information in conservation planning. Cross-amplification was demonstrated in two related taxa, including one of conservation concern.

Transcriptome Signatures of Selection, Drift, Introgression, and Gene Duplication in the Evolution of an Extremophile Endemic Plant

Plants on serpentine soils provide extreme examples of adaptation to environment, and thus offer excellent models for the study of evolution at the molecular and genomic level. Serpentine outcrops are derived from ultramafic rock and have extremely low levels of essential plant nutrients (e.g., N, P, K, and Ca), as well as toxic levels of heavy metals (e.g., Ni, Cr, and Co) and low moisture availability. These outcrops provide habitat to a number of endemic plant species, including the annual mustard Caulanthus amplexicaulis var. barbarae (Cab) (Brassicaceae). Its sister taxon, C. amplexicaulis var. amplexicaulis (Caa), is intolerant to serpentine soils. Here, we assembled and annotated comprehensive reference transcriptomes of both Caa and Cab for use in protein coding sequence comparisons. A set of 29,443 reciprocal best Blast hit (RBH) orthologs between Caa and Cab was compared with identify coding sequence variants, revealing a high genome-wide dN/dS ratio between the two taxa (mean = 0.346). We show that elevated dN/dS likely results from the composite effects of genetic drift, positive selection, and the relaxation of negative selection. Further, analysis of paralogs within each taxon revealed the signature of a period of elevated gene duplication (∼10 Ma) that is shared with other species of the tribe Thelypodieae, and may have played a role in the striking morphological and ecological diversity of this tribe. In addition, distribution of the synonymous substitution rate, dS, is strongly bimodal, indicating a history of reticulate evolution that may have contributed to serpentine adaptation.

Diversity analysis of cotton (Gossypium hirsutum L.) germplasm using the CottonSNP63K Array

BACKGROUND

Cotton germplasm resources contain beneficial alleles that can be exploited to develop germplasm adapted to emerging environmental and climate conditions. Accessions and lines have traditionally been characterized based on phenotypes, but phenotypic profiles are limited by the cost, time, and space required to make visual observations and measurements. With advances in molecular genetic methods, genotypic profiles are increasingly able to identify differences among accessions due to the larger number of genetic markers that can be measured. A combination of both methods would greatly enhance our ability to characterize germplasm resources. Recent efforts have culminated in the identification of sufficient SNP markers to establish high-throughput genotyping systems, such as the CottonSNP63K array, which enables a researcher to efficiently analyze large numbers of SNP markers and obtain highly repeatable results. In the current investigation, we have utilized the SNP array for analyzing genetic diversity primarily among cotton cultivars, making comparisons to SSR-based phylogenetic analyses, and identifying loci associated with seed nutritional traits.

RESULTS

The SNP markers distinctly separated G. hirsutum from other Gossypium species and distinguished the wild from cultivated types of G. hirsutum. The markers also efficiently discerned differences among cultivars, which was the primary goal when designing the CottonSNP63K array. Population structure within the genus compared favorably with previous results obtained using SSR markers, and an association study identified loci linked to factors that affect cottonseed protein content.

CONCLUSIONS

Our results provide a large genome-wide variation data set for primarily cultivated cotton. Thousands of SNPs in representative cotton genotypes provide an opportunity to finely discriminate among cultivated cotton from around the world. The SNPs will be relevant as dense markers of genome variation for association mapping approaches aimed at correlating molecular polymorphisms with variation in phenotypic traits, as well as for molecular breeding approaches in cotton.

QTL mapping for flowering-time and photoperiod insensitivity of cotton Gossypium darwinii Watt

Most wild and semi-wild species of the genus Gossypium are exhibit photoperiod-sensitive flowering. The wild germplasm cotton is a valuable source of genes for genetic improvement of modern cotton cultivars. A bi-parental cotton population segregating for photoperiodic flowering was developed by crossing a photoperiod insensitive irradiation mutant line with its pre-mutagenesis photoperiodic wild-type G. darwinii Watt genotype. Individuals from the F2 and F3 generations were grown with their parental lines and F1 hybrid progeny in the long day and short night summer condition (natural day-length) of Uzbekistan to evaluate photoperiod sensitivity, i.e., flowering-time during the seasons 2008-2009. Through genotyping the individuals of this bi-parental population segregating for flowering-time, linkage maps were constructed using 212 simple-sequence repeat (SSR) and three cleaved amplified polymorphic sequence (CAPS) markers. Six QTLs directly associated with flowering-time and photoperiodic flowering were discovered in the F2 population, whereas eight QTLs were identified in the F3 population. Two QTLs controlling photoperiodic flowering and duration of flowering were common in both populations. In silico annotations of the flanking DNA sequences of mapped SSRs from sequenced cotton (G. hirsutum L.) genome database has identified several potential 'candidate' genes that are known to be associated with regulation of flowering characteristics of plants. The outcome of this research will expand our understanding of the genetic and molecular mechanisms of photoperiodic flowering. Identified markers should be useful for marker-assisted selection in cotton breeding to improve early flowering characteristics.

Development, genetic mapping and QTL association of cotton PHYA, PHYB, and HY5-specific CAPS and dCAPS markers

BACKGROUND

Among SNP markers that become increasingly valuable in molecular breeding of crop plants are the CAPS and dCAPS markers derived from the genes of interest. To date, the number of such gene-based markers is small in polyploid crop plants such as allotetraploid cotton that has A- and D-sub-genomes. The objective of this study was to develop and map new CAPS and dCAPS markers for cotton developmental-regulatory genes that are important in plant breeding programs.

RESULTS

Gossypium hirsutum and G. barbadense, are the two cultivated allotetraploid cotton species. These have distinct fiber quality and other agronomic traits. Using comparative sequence analysis of characterized GSTs of the PHYA1, PHYB, and HY5 genes of G. hirsutum and G. barbadense one PHYA1-specific Mbo I/Dpn II CAPS, one PHYB-specific Alu I dCAPS, and one HY5-specific Hinf I dCAPS cotton markers were developed. These markers have successfully differentiated the two allotetraploid genomes (AD₁ and AD₂) when tested in parental genotypes of 'Texas Marker-1' ('TM-1'), 'Pima 3-79' and their F₁ hybrids. The genetic mapping and chromosome substitution line-based deletion analyses revealed that PHYA1 gene is located in A-sub-genome chromosome 11, PHYB gene is in A-sub-genome chromosome 10, and HY5 gene is in D-sub-genome chromosome 24, on the reference 'TM-1' x 'Pima 3-79' RIL genetic map. Further, it was found that genetic linkage map regions containing phytochrome and HY5-specific markers were associated with major fiber quality and flowering time traits in previously published QTL mapping studies.

CONCLUSION

This study detailed the genome mapping of three cotton phytochrome genes with newly developed CAPS and dCAPS markers. The proximity of these loci to fiber quality and other cotton QTL was demonstrated in two A-subgenome and one D-subgenome chromosomes. These candidate gene markers will be valuable for marker-assisted selection (MAS) programs to rapidly introgress G. barbadense phytochromes and/or HY5 gene (s) into G. hirsutum cotton genotypes or vice versa.

RNA Interference for Functional Genomics and Improvement of Cotton (Gossypium sp.)

RNA interference (RNAi), is a powerful new technology in the discovery of genetic sequence functions, and has become a valuable tool for functional genomics of cotton (Gossypium sp.). The rapid adoption of RNAi has replaced previous antisense technology. RNAi has aided in the discovery of function and biological roles of many key cotton genes involved in fiber development, fertility and somatic embryogenesis, resistance to important biotic and abiotic stresses, and oil and seed quality improvements as well as the key agronomic traits including yield and maturity. Here, we have comparatively reviewed seminal research efforts in previously used antisense approaches and currently applied breakthrough RNAi studies in cotton, analyzing developed RNAi methodologies, achievements, limitations, and future needs in functional characterizations of cotton genes. We also highlighted needed efforts in the development of RNAi-based cotton cultivars, and their safety and risk assessment, small and large-scale field trials, and commercialization.

Development of a 63K SNP Array for Cotton and High-Density Mapping of Intraspecific and Interspecific Populations of Gossypium spp

High-throughput genotyping arrays provide a standardized resource for plant breeding communities that are useful for a breadth of applications including high-density genetic mapping, genome-wide association studies (GWAS), genomic selection (GS), complex trait dissection, and studying patterns of genomic diversity among cultivars and wild accessions. We have developed the CottonSNP63K, an Illumina Infinium array containing assays for 45,104 putative intraspecific single nucleotide polymorphism (SNP) markers for use within the cultivated cotton species Gossypium hirsutum L. and 17,954 putative interspecific SNP markers for use with crosses of other cotton species with G. hirsutum. The SNPs on the array were developed from 13 different discovery sets that represent a diverse range of G. hirsutum germplasm and five other species: G. barbadense L., G. tomentosum Nuttal × Seemann, G. mustelinum Miers × Watt, G. armourianum Kearny, and G. longicalyx J.B. Hutchinson and Lee. The array was validated with 1,156 samples to generate cluster positions to facilitate automated analysis of 38,822 polymorphic markers. Two high-density genetic maps containing a total of 22,829 SNPs were generated for two F2 mapping populations, one intraspecific and one interspecific, and 3,533 SNP markers were co-occurring in both maps. The produced intraspecific genetic map is the first saturated map that associates into 26 linkage groups corresponding to the number of cotton chromosomes for a cross between two G. hirsutum lines. The linkage maps were shown to have high levels of collinearity to the JGI G. raimondii Ulbrich reference genome sequence. The CottonSNP63K array, cluster file and associated marker sequences constitute a major new resource for the global cotton research community.

SNP discovery in complex allotetraploid genomes (Gossypium spp., Malvaceae) using genotyping by sequencing

PREMISE OF THE STUDY

Single-nucleotide polymorphism (SNP) marker discovery in plants with complex allotetraploid genomes is often confounded by the presence of homeologous loci (along with paralogous and orthologous loci). Here we present a strategy to filter for SNPs representing orthologous loci.

METHODS AND RESULTS

Using Illumina next-generation sequencing, 54 million reads were collected from restriction enzyme-digested DNA libraries of a diversity of Gossypium taxa. Loci with one to three SNPs were discovered using the Stacks software package, yielding 25,529 new cotton SNP combinations, including those that are polymorphic at both interspecific and intraspecific levels. Frequencies of predicted dual-homozygous (aa/bb) marker polymorphisms ranged from 6.7-11.6% of total shared fragments in intraspecific comparisons and from 15.0-16.4% in interspecific comparisons.

CONCLUSIONS

This resource provides dual-homozygous (aa/bb) marker polymorphisms. Both in silico and experimental validation efforts demonstrated that these markers are enriched for single orthologous loci that are homozygous for alternative alleles.

Mapping genomic loci for cotton plant architecture, yield components, and fiber properties in an interspecific (Gossypium hirsutum L. × G. barbadense L.) RIL population

A quantitative trait locus (QTL) mapping was conducted to better understand the genetic control of plant architecture (PA), yield components (YC), and fiber properties (FP) in the two cultivated tetraploid species of cotton (Gossypium hirsutum L. and G. barbadense L.). One hundred and fifty-nine genomic regions were identified on a saturated genetic map of more than 2,500 SSR and SNP markers, constructed with an interspecific recombinant inbred line (RIL) population derived from the genetic standards of the respective cotton species (G. hirsutum acc. TM-1 × G. barbadense acc. 3-79). Using the single nonparametric and MQM QTL model mapping procedures, we detected 428 putative loci in the 159 genomic regions that confer 24 cotton traits in three diverse production environments [College Station F&B Road (FB), TX; Brazos Bottom (BB), TX; and Shafter (SH), CA]. These putative QTL loci included 25 loci for PA, 60 for YC, and 343 for FP, of which 3, 12, and 60, respectively, were strongly associated with the traits (LOD score ≥ 3.0). Approximately 17.7 % of the PA putative QTL, 32.9 % of the YC QTL, and 48.3 % of the FP QTL had trait associations under multiple environments. The At subgenome (chromosomes 1-13) contributed 72.7 % of loci for PA, 46.2 % for YC, and 50.4 % for FP while the Dt subgenome (chromosomes 14-26) contributed 27.3 % of loci for PA, 53.8 % for YC, and 49.6 % for FP. The data obtained from this study augment prior evidence of QTL clusters or gene islands for specific traits or biological functions existing in several non-homoeologous cotton chromosomes. DNA markers identified in the 159 genomic regions will facilitate further dissection of genetic factors underlying these important traits and marker-assisted selection in cotton.

De novo assembly and characterization of the transcriptome of the toxic dinoflagellate Karenia brevis

Background

Karenia brevis is a harmful algal species that blooms in the Gulf of Mexico and produces brevetoxins that cause neurotoxic shellfish poisoning. Elevated brevetoxin levels in K. brevis cells have been measured during laboratory hypo-osmotic stress treatments. To investigate mechanisms underlying K. brevis osmoacclimation and osmoregulation and establish a valuable resource for gene discovery, we assembled reference transcriptomes for three clones: Wilson-CCFWC268, SP3, and SP1 (a low-toxin producing variant). K. brevis transcriptomes were annotated with gene ontology terms and searched for putative transmembrane proteins that may elucidate cellular responses to hypo-osmotic stress. An analysis of single nucleotide polymorphisms among clones was used to characterize genetic divergence.

Results

K. brevis reference transcriptomes were assembled with 58.5 (Wilson), 78.0 (SP1), and 51.4 million (SP3) paired reads. Transcriptomes contained 86,580 (Wilson), 93,668 (SP1), and 84,309 (SP3) predicted transcripts. Approximately 40% of the transcripts were homologous to proteins in the BLAST nr database with an E value ≤ 1.00E-6. Greater than 80% of the highly conserved CEGMA core eukaryotic genes were identified in each transcriptome, which supports assembly completeness. Seven putative voltage-gated Na⁺ or Ca²⁺ channels, two aquaporin-like proteins, and twelve putative VATPase subunits were discovered in all clones using multiple bioinformatics approaches. Furthermore, 45% (Wilson) and 43% (SP1 and SP3) of the K. brevis putative peptides > 100 amino acids long produced significant hits to a sequence in the NCBI nr protein database. Of these, 77% (Wilson and SP1) and 73% (SP3) were successfully assigned gene ontology functional terms. The predicted single nucleotide polymorphism (SNP) frequencies between clones were 0.0028 (Wilson to SP1), 0.0030 (Wilson to SP3), and 0.0028 (SP1 to SP3).

Conclusions

The K. brevis transcriptomes assembled here provide a foundational resource for gene discovery and future RNA-seq experiments. The identification of ion channels, VATPases, and aquaporins in all three transcriptomes indicates that K. brevis regulates cellular ion and water concentrations via transmembrane proteins. Additionally, > 40,000 unannotated loci may include potentially novel K. brevis genes. Ultimately, the SNPs identified among the three ecologically diverse clones with different toxin profiles may help to elucidate variations in K. brevis brevetoxin production.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-888) contains supplementary material, which is available to authorized users.

Phytochrome RNAi enhances major fibre quality and agronomic traits of the cotton Gossypium hirsutum L

Simultaneous improvement of fibre quality, early-flowering, early-maturity and productivity in Upland cotton (G. hirsutum) is a challenging task for conventional breeding. The influence of red/far-red light ratio on the fibre length prompted us to examine the phenotypic effects of RNA interference (RNAi) of the cotton PHYA1 gene. Here we show a suppression of up to ~70% for the PHYA1 transcript, and compensatory overexpression of up to ~20-fold in the remaining phytochromes in somatically regenerated PHYA1 RNAi cotton plants. Two independent transformants of three generations exhibited vigorous root and vegetative growth, early-flowering, significantly improved upper half mean fibre length and an improvement in other major fibre characteristics. Small decreases in lint traits were observed but seed cotton yield was increased an average 10-17% compared with controls. RNAi-associated phenotypes were heritable and transferable via sexual hybridization. These results should aid in the development of early-maturing and productive Upland cultivars with superior fibre quality.

Clinical findings and survival in cats naturally infected with feline immunodeficiency virus

BACKGROUND

The clinical course and outcome of natural feline immunodeficiency virus (FIV) infection are variable and incompletely understood. Assigning clinical relevance to FIV infection in individual cats represents a considerable clinical challenge.

OBJECTIVE

To compare signalment, hematologic and biochemical data, major clinical problem, and survival among client-owned, FIV-infected, and uninfected domestic cats.

ANIMALS

Client-owned, domestic cats tested for FIV (n = 520).

METHODS

Retrospective, case control study. Logistic regression analyses were conducted to identify risk factors for FIV infection and to compare hematologic and biochemical data between cases and controls, after adjusting for potential confounders. Survival times were compared using Kaplan-Meier curves.

RESULTS

The prevalence of FIV infection was 14.6%. Mixed breed, male sex, and older age were risk factors for FIV infection. Hematologic abnormalities, biochemical abnormalities or both were common in both FIV-infected and uninfected cats. Lymphoid malignancies were slightly more common in FIV-infected than uninfected cats. Survival of FIV-infected cats was not significantly different from that of uninfected cats.

CONCLUSIONS AND CLINICAL IMPORTANCE

Multiple hematologic and biochemical abnormalities are common in old, sick cats regardless of their FIV status. Their presence should not be assumed to indicate clinical progression of FIV infection. A negative effect of FIV on survival was not apparent in this study.

Genetic analyses of nickel tolerance in a North American serpentine endemic plant, Caulanthus amplexicaulis var. barbarae (Brassicaceae)

PREMISE OF THE STUDY

The evolution of metal tolerance in plants is an important model for studies of adaptation to environment, population genetics, and speciation. Here, we investigated nickel tolerance in the North American serpentine endemic Caulanthus amplexicaulis var. barbarae in comparison with its nonserpentine sister taxon C. amplexicaulis var. amplexicaulis. We hypothesized that the serpentine endemic would have a heritable growth advantage on nickel-containing substrates.

METHODS

We employed an artificial growth assay to quantify biomass accumulation. Study plants were crossed to create an F(2:3) population that was used to determine the heritability of nickel tolerance and to map quantitative trait loci (QTL). Nickel accumulation in both laboratory populations and native specimens was examined using energy-dispersive x-ray fluorescence (EDXRF).

KEY RESULTS

The serpentine endemic had a dramatic growth advantage at concentrations of nickel >30 µmol/L. Caulanthus amplexicaulis var. barbarae and its nonserpentine sister taxon both accumulated nickel to substantial levels. Nickel tolerance was highly heritable (h(2) = 0.59) and not associated with accumulation. The QTL analyses identified two major loci for nickel tolerance, on linkage group 2 (LG2) and linkage group 9 (LG9).

CONCLUSIONS

In our study, nickel tolerance was determined by two major loci with large effects. At both loci, alleles from the serpentine parent conferred positive effects on nickel tolerance, suggesting that they are adaptive in the natural serpentine environment. The mechanism of nickel tolerance in the serpentine plant was not exclusion of nickel. Nickel tolerance may have an inducible component in C. amplexicaulis var. barbarae.

Isolation and characterization of microsatellite markers for Bothriochloa ischaemum (Poaceae)

PREMISE OF THE STUDY

Microsatellite primers were developed for Bothriochloa ischaemum to investigate the structure of invasive populations within Texas and determine the origin of introduction from within the native range.

METHODS AND RESULTS

We adapted the biotinylated nucleotide method of developing an enriched genomic library to isolate and characterize 10 polymorphic microsatellite markers. The number of alleles per locus (A) ranged from five to 18 (mean A = 10.45), expected heterozygosity (H(E)) ranged from 0.620 to 0.895 (mean H(E) = 0.785), and expected heterozygosity corrected for sample size (H(EC)) ranged from 0.635 to 0.909 (mean H(EC) = 0.799). The primers were also tested for amplification in Schizachyrium scoparium var. scoparium, Andropogon gerardii, Bothriochloa saccharoides, and Dichanthium annulatum.

CONCLUSIONS

The use of microsatellite markers may assist in understanding the pattern of spread, determining the source of invasive populations, and developing biological control agents for invasive populations of Bothriochloa ischaemum.

Duplication, divergence and persistence in the Phytochrome photoreceptor gene family of cottons (Gossypium spp.)

BACKGROUND

Phytochromes are a family of red/far-red photoreceptors that regulate a number of important developmental traits in cotton (Gossypium spp.), including plant architecture, fiber development, and photoperiodic flowering. Little is known about the composition and evolution of the phytochrome gene family in diploid (G. herbaceum, G. raimondii) or allotetraploid (G. hirsutum, G. barbadense) cotton species. The objective of this study was to obtain a preliminary inventory and molecular-evolutionary characterization of the phytochrome gene family in cotton.

RESULTS

We used comparative sequence resources to design low-degeneracy PCR primers that amplify genomic sequence tags (GSTs) for members of the PHYA, PHYB/D, PHYC and PHYE gene sub-families from A- and D-genome diploid and AD-genome allotetraploid Gossypium species. We identified two paralogous PHYA genes (designated PHYA1 and PHYA2) in diploid cottons, the result of a Malvaceae-specific PHYA gene duplication that occurred approximately 14 million years ago (MYA), before the divergence of the A- and D-genome ancestors. We identified a single gene copy of PHYB, PHYC, and PHYE in diploid cottons. The allotetraploid genomes have largely retained the complete gene complements inherited from both of the diploid genome ancestors, with at least four PHYA genes and two genes encoding PHYB, PHYC and PHYE in the AD-genomes. We did not identify a PHYD gene in any cotton genomes examined.

CONCLUSIONS

Detailed sequence analysis suggests that phytochrome genes retained after duplication by segmental duplication and allopolyploidy appear to be evolving independently under a birth-and-death-process with strong purifying selection. Our study provides a preliminary phytochrome gene inventory that is necessary and sufficient for further characterization of the biological functions of each of the cotton phytochrome genes, and for the development of 'candidate gene' markers that are potentially useful for cotton improvement via modern marker-assisted selection strategies.

Linkage disequilibrium based association mapping of fiber quality traits in G. hirsutum L. variety germplasm

Cotton is the world's leading cash crop, but it lags behind other major crops for marker-assisted breeding due to limited polymorphisms and a genetic bottleneck through historic domestication. This underlies a need for characterization, tagging, and utilization of existing natural polymorphisms in cotton germplasm collections. Here we report genetic diversity, population characteristics, the extent of linkage disequilibrium (LD), and association mapping of fiber quality traits using 202 microsatellite marker primer pairs in 335 G. hirsutum germplasm grown in two diverse environments, Uzbekistan and Mexico. At the significance threshold (r (2) >or= 0.1), a genome-wide average of LD extended up to genetic distance of 25 cM in assayed cotton variety accessions. Genome wide LD at r (2) >or= 0.2 was reduced to approximately 5-6 cM, providing evidence of the potential for association mapping of agronomically important traits in cotton. Results suggest linkage, selection, inbreeding, population stratification, and genetic drift as the potential LD-generating factors in cotton. In two environments, an average of ~20 SSR markers was associated with each main fiber quality traits using a unified mixed liner model (MLM) incorporating population structure and kinship. These MLM-derived significant associations were confirmed in general linear model and structured association test, accounting for population structure and permutation-based multiple testing. Several common markers, showing the significant associations in both Uzbekistan and Mexican environments, were determined. Between 7 and 43% of the MLM-derived significant associations were supported by a minimum Bayes factor at 'moderate to strong' and 'strong to very strong' evidence levels, suggesting their usefulness for marker-assisted breeding programs and overall effectiveness of association mapping using cotton germplasm resources.

Molecular diversity and association mapping of fiber quality traits in exotic G. hirsutum L. germplasm

The narrow genetic base of cultivated cotton germplasm is hindering the cotton productivity worldwide. Although potential genetic diversity exists in Gossypium genus, it is largely 'underutilized' due to photoperiodism and the lack of innovative tools to overcome such challenges. The application of linkage disequilibrium (LD)-based association mapping is an alternative powerful molecular tool to dissect and exploit the natural genetic diversity conserved within cotton germplasm collections, greatly accelerating still 'lagging' cotton marker-assisted selection (MAS) programs. However, the extent of genome-wide linkage disequilibrium (LD) has not been determined in cotton. We report the extent of genome-wide LD and association mapping of fiber quality traits by using a 95 core set of microsatellite markers in a total of 285 exotic Gossypium hirsutum accessions, comprising of 208 landrace stocks and 77 photoperiodic variety accessions. We demonstrated the existence of useful genetic diversity within exotic cotton germplasm. In this germplasm set, 11-12% of SSR loci pairs revealed a significant LD. At the significance threshold (r(2)>/=0.1), a genome-wide average of LD declines within the genetic distance at <10 cM in the landrace stocks germplasm and >30 cM in variety germplasm. Genome wide LD at r(2)>/=0.2 was reduced on average to approximately 1-2 cM in the landrace stock germplasm and 6-8 cM in variety germplasm, providing evidence of the potential for association mapping of agronomically important traits in cotton. We observed significant population structure and relatedness in assayed germplasm. Consequently, the application of the mixed liner model (MLM), considering both kinship (K) and population structure (Q) detected between 6% and 13% of SSR markers associated with the main fiber quality traits in cotton. Our results highlight for the first time the feasibility and potential of association mapping, with consideration of the population structure and stratification existing in cotton germplasm resources. The number of SSR markers associated with fiber quality traits in diverse cotton germplasm, which broadly covered many historical meiotic events, should be useful to effectively exploit potentially new genetic variation by using MAS programs.

Latitudinal variation in cold hardiness in introduced Tamarix and native Populus

To investigate the evolution of clinal variation in an invasive plant, we compared cold hardiness in the introduced saltcedar (Tamarix ramosissima, Tamarix chinensis, and hybrids) and the native plains cottonwood (Populus deltoides subsp. monilifera). In a shadehouse in Colorado (41°N), we grew plants collected along a latitudinal gradient in the central United States (29-48°N). On 17 occasions between September 2005 and June 2006, we determined killing temperatures using freeze-induced electrolyte leakage and direct observation. In midwinter, cottonwood survived cooling to -70°C, while saltcedar was killed at -33 to -47°C. Frost sensitivity, therefore, may limit northward expansion of saltcedar in North America. Both species demonstrated inherited latitudinal variation in cold hardiness. For example, from September through January killing temperatures for saltcedar from 29.18°N were 5-21°C higher than those for saltcedar from 47.60°N, and on September 26 and October 11, killing temperatures for cottonwood from 33.06°N were >43°C higher than those for cottonwood from 47.60°N. Analysis of nine microsatellite loci showed that southern saltcedars are more closely related to T. chinensis while northern plants are more closely related to T. ramosissima. Hybridization may have introduced the genetic variability necessary for rapid evolution of the cline in saltcedar cold hardiness.

The Arabidopsis thaliana carboxyl-terminal domain phosphatase-like 2 regulates plant growth, stress and auxin responses

More than 20 genes in the Arabidopsis genome encode proteins similar to phosphatases that act on the carboxyl-terminal domain (CTD) of RNA polymerase II. One of these CTD-phosphatase-like (CPL) proteins, CPL2, dephosphorylates CTD-Ser5-PO4 in an intact RNA polymerase II complex and contains a double-stranded (ds)-RNA-binding motif (DRM). Although the dsRNA-binding activity of CPL2 DRM has not been shown to date, T-DNA insertion mutants that express CPL2 variants lacking either a part of DRM (cpl2-1) or the entire DRM (cpl2-2) exhibited leaf expansion defects, early flowering, low fertility, and increased salt sensitivity. cpl2 mutant plants produced shorter hypocotyls than wild-type plants in the light, but were indistinguishable from wild type in the dark. CPL2 was expressed in shoot and root meristems and vasculatures, expanding rosette leaves, and floral organs suggesting a focal role for growth. Microarray and RT-PCR analyses revealed that basal levels of several auxin-responsive transcripts were reduced in cpl2. On the other hand, the levels of endogenous auxin and its conjugates were similar in wild type and cpl2. Overexpression of ARF5 but not all activator ARF transcription factors restored the auxin-responsive DR5-GUS reporter gene expression and the leaf expansion of cpl2 mutant plants but not early flowering phenotype. These results establish CPL2 as a multifunctional regulator that modulates plant growth, stress, and auxin responses.

The role of induced mutation in conversion of photoperiod dependence in cotton

Wild cotton germplasm resources are largely underutilized because of photoperiod-dependent flowering of "exotic" cottons. The objectives of this work were to explore the genome-wide effect of induced mutation in photoperiod-converted induced cotton mutants, estimating the genetic change between mutant and wild-type cottons using simple sequence repeats (SSRs) as well as understand the pattern of SSR mutation in induced mutagenesis. Three groups of photoperiod-converted radiomutants ((32)P) including their wild-type parental lines, A- and D-genome diploids, and typically grown cotton cultivars were screened with 250 cotton SSR primer pairs. Forty SSRs revealed the same SSR mutation profile in, at least, 2 independent mutant lines that were different from the original wild types. Induced mutagenesis both increased and decreased the allele sizes of SSRs in mutants with the higher mutation rate in SSRs containing dinucleotide motifs. Genetic distance obtained based on 141 informative SSR alleles ranged from 0.09 to 0.60 in all studied cotton genotypes. Genetic distance within all photoperiod-converted induced mutants was in a 0.09-0.25 range. The genetic distance among photoperiod-converted mutants and their originals ranged from 0.28 to 0.50, revealing significant modification of mutants from their original wild types. Typical Gossypium hirsutum cultivar, Namangan-77, revealed mutational pattern similar to induced radiomutants in 40 mutated SSR loci, implying possible pressure to these SSR loci not only in radiomutagenesis but also during common breeding process. Outcomes of the research should be useful in understanding the photoperiod-related mutations, and markers might help in mapping photoperiodic flowering genes in cotton.

Evolution of genome size in Brassicaceae

BACKGROUND AND AIMS

Brassicaceae, with nearly 340 genera and more than 3350 species, anchors the low range of angiosperm genome sizes. The relatively narrow range of DNA content (0.16 pg < 1C < 1.95 pg) was maintained in spite of extensive chromosomal change. The aim of this study was to erect a cytological and molecular phylogenetic framework for a selected subset of the Brassicacae, and use this as a template to examine genome size evolution in Brassicaceae.

METHODS

DNA contents were determined by flow cytometry and chromosomes were counted for 34 species of the family Brassicaceae and for ten Arabidopsis thaliana ecotypes. The amplified and sequenced ITS region for 23 taxa (plus six other taxa with known ITS sequences) were aligned and used to infer evolutionary relationship by parsimony analysis.

KEY RESULTS

DNA content in the species studied ranged over 8-fold (1C = 0.16-1.31 pg), and 4.4-fold (1C = 0.16-0.71 pg) excluding allotetraploid Brassica species. The 1C DNA contents of ten Arabidopsis thaliana ecotypes showed little variation, ranging from 0.16 pg to 0.17 pg.

CONCLUSIONS

The tree roots at an ancestral genome size of approximately 1x = 0.2 pg. Arabidopsis thaliana (1C = 0.16 pg; approximately 157 Mbp) has the smallest genome size in Brassicaceae studied here and apparently represents an evolutionary decrease in genome size. Two other branches that represent probable evolutionary decreases in genome size terminate in Lepidium virginicum and Brassica rapa. Branches in the phylogenetic tree that represent probable evolutionary increases in genome size terminate in Arabidopsis halleri, A. lyrata, Arabis hirsuta, Capsella rubella, Caulanthus heterophyllus, Crucihimalaya, Lepidium sativum, Sisymbrium and Thlaspi arvense. Branches within one clade containing Brassica were identified that represent two ancient ploidy events (2x to 4x and 4x to 6x) that were predicted from published comparative mapping studies.

shl, a New set of Arabidopsis mutants with exaggerated developmental responses to available red, far-red, and blue light

The interaction of light perception with development is the subject of intensive genetic analysis in the model plant Arabidopsis. We performed genetic screens in low white light-a threshold condition in which photomorphogenetic signaling pathways are only partially active-for ethyl methane sulfonate-generated mutants with altered developmental phenotypes. Recessive mutants with exaggerated developmental responses were obtained in eight complementation groups designated shl for seedlings hyperresponsive to light. shl1, shl2, shl5, and shl3 shl4 (double mutant) seedlings showed limited or no phenotypic effects in darkness, but showed significantly enhanced inhibition of hypocotyl elongation in low-white, red, far-red, blue, and green light across a range of fluences. These results reflect developmental hyper-responsiveness to signals generated by both phytochrome and cryptochrome photoreceptors. The shl11 mutant retained significant phenotypic effects on hypocotyl length in both the phyA mutant and phyB mutant backgrounds but may be dependent on CRY1 for phenotypic expression in blue light. The shl2 phenotype was partially dependent on PHYB, PHYA, and CRY1 in red, far-red, and blue light, respectively. shl2 and, in particular, shl1 were partially dependent on HY5 activity for their light-hyperresponsive phenotypes. The SHL genes act (genetically) as light-dependent negative regulators of photomorphogenesis, possibly in a downstream signaling or developmental pathway that is shared by CRY1, PHYA, and PHYB and other photoreceptors (CRY2, PHYC, PHYD, and PHYE).

Evolution of Caulanthus amplexicaulis var. barbarae (Brassicaceae), a rare serpentine endemic plant: a molecular phylogenetic perspective

Intra- and interspecific phylogenetic relationships of the rare serpentine endemic taxon Caulanthus amplexicaulus var. barbarae and related taxa in the "Streptanthoid Complex" of genera (Streptanthus, Caulanthus, Guillenia) were examined using nuclear ribosomal internal transcribed spacer (ITS) and chloroplast trnL intron sequences. Phylogenetic hypotheses generated from 81 variable ITS nucleotide sites and six variable trnL nucleotide sites indicate that Streptanthus and Caulanthus are nonmonophyletic groups. Caulanthus amplexicaulis var. barbarae and its more widespread nonserpentine sister taxon Caulanthus amplexicaulis var. amplexicaulis formed a distinct monophyletic group. Among the taxa in our study, C. amplexicaulis was most closely related to Streptanthus tortuosus. The ITS sequences supported monophyly of subgenus Euclesia, which includes the bulk of the serpentine endemics in the Streptanthoid Complex. The serpentine taxa were nonmonophyletic, occurring in at least three distinct clades, suggesting that tolerance to serpentine may be gained or lost through relatively few genetic changes. Intraspecific ITS1 and ITS2 sequence divergence within C. amplexicaulis (1.3-1.8%) was higher than in comparable species (0.0-0.3%); implications of this genetic differentiation for the conservation status of C. amplexicaulis var. barbarae are discussed. Evidence is presented that supports a "biotype depletion" model for the origin of this rare endemic taxon.

dCAPS, a simple technique for the genetic analysis of single nucleotide polymorphisms: experimental applications in Arabidopsis thaliana genetics

PCR-based detection of single nucleotide polymorphisms is a powerful tool for the plant geneticist. Cleaved amplified polymorphic sequence analysis is the most widely used approach for the detection of single nucleotide polymorphisms. However, this technique is limited to mutations which create or disrupt a restriction enzyme recognition site. This paper presents a modification of this technique where mismatches in a PCR primer are used to create a polymorphism based on the target mutation. This technique is useful for following known mutations in segregating populations and genetic mapping of isolated DNAs used for positional based cloning of new genes. In addition, a computer program has been developed that facilitates the design of these PCR primers.

dCAPS, a simple technique for the genetic analysis of single nucleotide polymorphisms: experimental applications in Arabidopsis thaliana genetics

PCR-based detection of single nucleotide polymorphisms is a powerful tool for the plant geneticist. Cleaved amplified polymorphic sequence analysis is the most widely used approach for the detection of single nucleotide polymorphisms. However, this technique is limited to mutations which create or disrupt a restriction enzyme recognition site. This paper presents a modification of this technique where mismatches in a PCR primer are used to create a polymorphism based on the target mutation. This technique is useful for following known mutations in segregating populations and genetic mapping of isolated DNAs used for positional based cloning of new genes. In addition, a computer program has been developed that facilitates the design of these PCR primers.

Carrier and prenatal diagnosis of X-linked severe combined immunodeficiency: mutation detection methods and utilization

IL2RG, the gene encoding the common gamma chain, gamma c, of the receptor for interleukin-2 and other cytokines, has been identified as the disease gene for severe combined immunodeficiency (SCID) of the X-linked type. Specific mutational diagnosis for X-linked SCID has thus become possible. For many women at risk for carrying an IL2RG mutation, no samples were saved from an affected male relative prior to either death or bone marrow transplantation (BMT). To establish optimal methods for genetic evaluation of such women, we compared mutational screening by single-strand conformational polymorphism, heteroduplex analysis and dideoxy fingerprinting (ddF). Abnormally migrating band patterns were followed up with direct sequencing for identification of specific mutations. The most sensitive method, ddF, detected heterozygous alterations, subsequently confirmed to represent significant mutations, in all of 19 unrelated obligate or suspected carriers studied. Some of these women, as well as others at risk for carrying an X-linked SCID mutation, enrolled in a study of prenatal diagnosis after fetal testing for gender determination. Originally using linkage analysis and, more recently, specific detection of IL2RG mutations, we evaluated pregnancies at risk for X-linked SCID prospectively on a research basis. Of 27 male fetuses tested 14 were predicted to be unaffected and confirmed to have normal immune status at birth. Among pregnancies predicted to be affected, 2 were terminated, while 11 affected males were born at term. Nine of these received neonatal BMT, one had BMT at 3 months of age, and one underwent a successful experimental in utero BMT. In our study cohort accurate prenatal diagnosis assisted decision making and expanded treatment options for families at risk for having infants with a severe, but treatable genetic disorder that presents early in life.

Extragenic suppressors of the Arabidopsis det1 mutant identify elements of flowering-time and light-response regulatory pathways

Light regulation of seeding morphogenesis is mediated by photoreceptors that perceive red, far-red, blue and UV light. Photomorphogenetic mutants of Arabidopsis have identified several of the primary photoreceptors, as well as a set of negative regulators of seedling photomorphogenesis, including DET1, that appear to act downstream of the photoreceptors. To study the regulatory context in which DET1 acts to repress photomorphogenesis, we used a simple morphological screen to isolate extragenic mutations in six loci, designated ted (for reversal of the det phenotype), that partially or fully suppress the seeding morphological phenotype of det1-1. Genetic analyses indicate that mutations in the ted4 and ted5 loci identify new alleles of the previously described photomorphogenetic loci hy1 and hy5, respectively. Molecular analyses indicate that the ted mutations partially suppress the dark-grown gene expression phenotype of det1-1, and that the mechanism of suppression does not involve direct remediation of the splicing defect caused by the det1-1 mutation. The ted mutations also partially suppress the light-grown morphological phenotype of mature det1-1 plants, and ted1 and ted2 suppress a daylength insensitivity phenotype of det1. TED1, TED2 and TED3 are newly described genes, whose function appears closely associated with that of DET1. In addition, alleles of ted1 are associated with a moderate late-flowering phenotype, suggesting that TED1 plays a role in the pathways that regulate both seedling morphogenesis and the initiation of flowering.

Mutation analysis of IL2RG in human X-linked severe combined immunodeficiency

Severe combined immunodeficiency (SCID) is a syndrome of profoundly impaired cellular and humoral immunity. In humans, SCID is most commonly caused by mutations in the X-linked gene IL2RG, which encodes the common gamma chain, gamma c, of the leukocyte receptors for interleukin-2 and multiple other cytokines. To investigate the frequency and variety of IL2RG mutations that cause SCID, we analyzed DNA, RNA, and B-cell lines from a total of 103 unrelated SCID-affected males and their relatives using a combination of molecular and immunologic techniques. Sixty-two different mutations spanning all eight IL2RG exons were found in 87 cases, making possible correlations between mutation type and functional consequences. Although skewed maternal X chromosome inactivation, single-strand conformation polymorphism, mRNA expression, and cell surface staining with anti-gamma c antibodies were all helpful in establishing IL2RG defects as the cause of SCID, only dideoxy fingerprinting and DNA sequence determination each detected 100% of the IL2RG mutations in our series. Abnormal gamma c chains may be expressed in the lymphocytes of as many as two thirds of patients with X-linked SCID. Specific mutation diagnosis thus remains technically challenging, but it is important for genetic counseling and perhaps for helping to select appropriate subjects for retroviral gene therapy trials, This is a US government work. There are no restrictions on its use.

Maternal mosaicism for a novel interleukin-2 receptor gamma-chain mutation causing X-linked severe combined immunodeficiency in a Navajo kindred

X-linked severe combined immunodeficiency disease (SCID) results from mutations of IL2RG, the gene encoding the interleukin-2 receptor gamma chain, also known as the common gamma chain (gamma c). A distinct form of autosomal recessive SCID occurs at an increased frequency among the Navajo Native American population. The disease gene responsible for autosomal Navajo SCID remains to be determined. We report the occurrence of X-linked SCID in a Navajo Native American kindred with two affected brothers. X-linked SCID was suggested by the presence of circulating B cells and the absence of surface gamma c expression in a cell line derived from an affected male. A C-to-T transition was demonstrated in exon 5 of the IL2RG gene, resulting in the substitution of tryptophan for arginine at position 224. This change was not present in the peripheral blood lymphocytes of the mother, thus proving the occurrence of a new mutation in the maternal germline. This report underscores the importance of establishing a specific genetic diagnosis for SCID cases and illustrates the inherent difficulties in providing genetic counseling in cases involving mosaicism.

Two mutational hotspots in the interleukin-2 receptor gamma chain gene causing human X-linked severe combined immunodeficiency

Human severe combined immunodeficiency (SCID), a syndrome of profoundly impaired cellular and humoral immunity, is most commonly caused by mutations in the X-linked gene for interleukin-2 (IL-2) receptor gamma chain (IL2RG). For mutational analysis of IL2RG in males with SCID, SSCP screening was followed by DNA sequencing. Of 40 IL2RG mutations found in unrelated SCID patients, 6 were point mutations at the CpG dinucleotide at cDNA 690-691, encoding amino acid R226. This residue lies in the extracellular domain of the protein in a region not previously recognized to be significantly conserved in the cytokine receptor gene family, 11 amino acids upstream from the highly conserved WSXWS motif. Three additional instances of mutation at another CpG dinucleotide at cDNA 879 produced a premature termination signal in the intracellular domain of IL2RG, resulting in loss of the SH2-homologous intracellular domain known to be essential for signaling from the IL-2 receptor complex. Mutations at these two hotspots constitute > 20% of the X-linked SCID mutations found by our group and a similar proportion of all reported IL2RG mutations.

The GCR1 requirement for yeast glycolytic gene expression is suppressed by dominant mutations in the SGC1 gene, which encodes a novel basic-helix-loop-helix protein

The GCR1 gene product is required for maximal transcription of yeast glycolytic genes and for growth of yeast strains in media containing glucose as a carbon source. Dominant mutations in two genes, SGC1 and SGC2, as well as recessive mutations in the SGC5 gene were identified as suppressors of the growth and transcriptional defects caused by a gcr1 null mutation. The wild-type and mutant alleles of SGC1 were cloned and sequenced. The predicted amino acid sequence of the SGC1 gene product includes a region with substantial similarity to the basic-helix-loop-helix domain of the Myc family of DNA-binding proteins. The SGC1-1 dominant mutant allele contained a substitution of glutamine for a highly conserved glutamic acid residue within the putative basic DNA binding domain. A second dominant mutant, SGC1-2, contained a valine-for-isoleucine substitution within the putative loop region. The SGC1-1 dominant mutant suppressed the GCR1 requirement for enolase, glyceraldehyde-3-phosphate dehydrogenase, phosphoglycerate kinase, phosphoglycerate mutase, and pyruvate kinase gene expression. Expression of the yeast enolase genes was reduced three- to fivefold in strains carrying an sgc1 null mutation, demonstrating that SGC1 is required for maximal enolase gene expression. Expression of the enolase genes in strains carrying gcr1 and sgc1 double null mutations was substantially less than observed for strains carrying either null mutation alone, suggesting that GCR1 and SGC1 function on parallel pathways to activate yeast glycolytic gene expression.

Female germ line mosaicism as the origin of a unique IL-2 receptor gamma-chain mutation causing X-linked severe combined immunodeficiency

The IL2RG gene encoding the gamma chain of the lymphocyte receptor for IL-2 lies in human Xq13.1 and is mutated in males with X-linked severe combined immunodeficiency (SCID). In a large Canadian pedigree genetic linkage studies demonstrated that the proband's grandmother was the source of an X-linked SCID mutation. However, her T cells did not show the expected skewed X chromosome inactivation pattern of female carriers of SCID, despite her having one affected son and two carrier daughters with skewed X inactivation. Single strand conformation polymorphism analysis of IL2RG in the affected proband was abnormal in exon 5; sequencing revealed a nine nucleotide in-frame duplication insertion. The three duplicated amino acids included the first tryptophan of the "WSXWS" motif found in all members of the cytokine receptor gene superfamily. Mutation detection in the pedigree confirmed that the founder grandmother's somatic cells had only normal IL2RG, and further showed that the SCID-associated X chromosome haplotype was inherited by three daughters, one with a wild type IL2RG gene and two others with the insertional mutation. Female germ line mosaicism is unusual, but its presence in this X-linked SCID family emphasizes the limitations of genetic diagnosis by linkage as compared with direct mutation analysis.

Transcription of the constitutively expressed yeast enolase gene ENO1 is mediated by positive and negative cis-acting regulatory sequences

There are two enolase genes, ENO1 and ENO2, per haploid yeast genome. Expression of the ENO1 gene is quantitatively similar in cells grown on glucose or gluconeogenic carbon sources. In contrast, ENO2 expression is induced more than 20-fold in cells grown on glucose as the carbon source. cis-Acting regulatory sequences were mapped within the 5'-flanking region of the constitutively expressed yeast enolase gene ENO1. A complex positive regulatory region was located 445 base pairs (bp) upstream from the transcriptional initiation site which was required for ENO1 expression in cells grown on glycolytic or gluconeogenic carbon sources. A negative regulatory region was located 160 bp upstream from the transcriptional initiation site. Sequences required for the function of this negative regulatory element were mapped to a 38-bp region. Deletion of all or a portion of these latter sequences permitted glucose-dependent induction of ENO1 expression that was quantitatively similar to that of the glucose-inducible ENO2 gene. The negative regulatory element therefore prevents glucose-dependent induction of the ENO1 gene. Hybrid 5'-flanking regions were constructed which contained the upstream regulatory sequences of one enolase gene fused at a site upstream from the TATAAA box in the other enolase gene. Analysis of the expression of enolase genes containing these hybrid 5'-flanking region showed that the positive regulatory regions of ENO1 and ENO2 were functionally similar, as were the regions extending from the TATAAA boxes to the initiation codons. Based on these studies, we conclude that the negative regulatory element plays the critical role in maintaining the constitutive expression of the ENO1 structural gene in cells grown on glucose or gluconeogenic carbon sources.