AP-1/c-Fos supports SIV and HIV-1 latency in CD4 T cells infected in vivo.

Persistent HIV-1 reservoirs of infected CD4 T cells are a major barrier to HIV-1 cure, although the mechanisms by which they are established and maintained in vivo remain poorly characterized. To elucidate host cell gene expression patterns that govern virus gene expression, we analyzed viral RNA+ (vRNA) CD4 T cells of untreated simian immunodeficiency virus (SIV)-infected macaques by single-cell RNA sequencing. A subset of vRNA+ cells distinguished by spliced and high total vRNA (7-10% of reads) expressed diminished FOS, a component of the Activator protein 1 (AP-1) transcription factor, relative to vRNA-low and -negative cells. Conversely, FOS and JUN, another AP-1 component, were upregulated in HIV DNA+ infected cells compared to uninfected cells from people with HIV-1 on suppressive therapy. Inhibiting c-Fos in latently infected primary cells augmented reactivatable HIV-1 infection. These findings implicate AP-1 in latency establishment and maintenance and as a potential therapeutic target to limit HIV-1 reservoirs.

Single-cell transcriptomics identifies prothymosin α restriction of HIV-1 in vivo.

Host restriction factors play key roles in innate antiviral defense, but it remains poorly understood which of them restricts HIV-1 in vivo. Here, we used single-cell transcriptomic analysis to identify host factors associated with HIV-1 control during acute infection by correlating host gene expression with viral RNA abundance within individual cells. Wide sequencing of cells from one participant with the highest plasma viral load revealed that intracellular viral RNA transcription correlates inversely with expression of the gene PTMA, which encodes prothymosin α. This association was genome-wide significant (Padjusted < 0.05) and was validated in 28 additional participants from Thailand and the Americas with HIV-1 CRF01_AE and subtype B infections, respectively. Overexpression of prothymosin α in vitro confirmed that this cellular factor inhibits HIV-1 transcription and infectious virus production. Our results identify prothymosin α as a host factor that restricts HIV-1 infection in vivo, which has implications for viral transmission and cure strategies.

HLA-B*57 and B*58 Associate with Predictors of Reservoir Size in an Acutely Treated HIV Cohort.

Much has been learnt about the role of human leukocyte antigen (HLA) alleles during natural infection of HIV-1, but far less is known about their role in people living with HIV (PLWH) on suppressive antiretroviral therapy (ART). In this study we used variable selection to identify predictors of HIV reservoir size, as measured by total HIV DNA in 192 participants in an acute HIV infection (AHI) cohort. Baseline clinical data including pre-ART CD4 T cell counts and plasma viral load (VL) were available from all participants along with longitudinal measurements after ART initiation during AHI. Time to VL suppression, time to CD4 reconstitution, and pre-ART viremia were the strongest predictors of undetectable total HIV DNA at 24 weeks after ART initiation. We next performed HLA typing in 526 participants from the same cohort and investigated associations with the three predictors of reservoir size. HLA-B*57 and B*58 both associated significantly with time to VL suppression, which was one of the predictors of the size of the HIV reservoir. These findings are significant in PLWH and have to be considered in the context of therapeutic intervention when conducting analytic treatment interruption studies as participants with these alleles could impact clinical findings given the small sizes of these studies.

Long-term antiretroviral therapy initiated in acute HIV infection prevents residual dysfunction of HIV-specific CD8+ T cells.

BACKGROUND: Harnessing CD8+ T cell responses is being explored to achieve HIV remission. Although HIV-specific CD8+ T cells become dysfunctional without treatment, antiretroviral therapy (ART) partially restores their function. However, the extent of this recovery under long-term ART is less understood. METHODS: We analyzed the differentiation status and function of HIV-specific CD8+ T cells after long-term ART initiated in acute or chronic HIV infection ex vivo and upon in vitro recall. FINDINGS: ART initiation in any stage of acute HIV infection promoted the persistence of long-lived HIV-specific CD8+ T cells with high expansion (P<0·0008) and cytotoxic capacity (P=0·02) after in vitro recall, albeit at low cell number (P=0·003). This superior expansion capacity correlated with stemness (r=0·90, P=0·006), measured by TCF-1 expression, similar to functional HIV-specific CD8+ T cells found in spontaneous controllers. Importanly, TCF-1 expression in these cells was associated with longer time to viral rebound ranging from 13 to 48 days after ART interruption (r =0·71, P=0·03). In contrast, ART initiation in chronic HIV infection led to more differentiated HIV-specific CD8+ T cells lacking stemness properties and exhibiting residual dysfunction upon recall, with reduced proliferation and cytolytic activity. INTERPRETATION: ART initiation in acute HIV infection preserves functional HIV-specific CD8+ T cells, albeit at numbers too low to control viral rebound post-ART. HIV remission strategies may need to boost HIV-specific CD8+ T cell numbers and induce stem cell-like properties to reverse the residual dysfunction persisting on ART in people treated after acute infection prior to ART release. FUNDING: U.S. National Institutes of Health and U.S. Department of Defense.

HLA-B∗46 associates with rapid HIV disease progression in Asian cohorts and prominent differences in NK cell phenotype.

Human leukocyte antigen (HLA) alleles have been linked to HIV disease progression and attributed to differences in cytotoxic T lymphocyte (CTL) epitope representation. These findings are largely based on treatment-naive individuals of European and African ancestry. We assessed HLA associations with HIV-1 outcomes in 1,318 individuals from Thailand and found HLA-B∗46:01 (B∗46) associated with accelerated disease in three independent cohorts. B∗46 had no detectable effect on HIV-specific T cell responses, but this allele is unusual in containing an HLA-C epitope that binds inhibitory receptors on natural killer (NK) cells. Unbiased transcriptomic screens showed increased NK cell activation in people with HIV, without B∗46, and simultaneous single-cell profiling of surface proteins and transcriptomes revealed a NK cell subset primed for increased responses in the absence of B∗46. These findings support a role for NK cells in HIV pathogenesis, revealed by the unique properties of the B∗46 allele common only in Asia.

Associations of human leukocyte antigen with neutralizing antibody titers in a tetravalent dengue vaccine phase 2 efficacy trial in Thailand.

The recombinant, live, attenuated, tetravalent dengue vaccine CYD-TDV has shown efficacy against all four dengue serotypes. In this exploratory study (CYD59, NCT02827162), we evaluated potential associations of host human leukocyte antigen (HLA) alleles with dengue antibody responses, CYD-TDV vaccine efficacy, and virologically-confirmed dengue (VCD) cases. Children 4-11 years old, who previously completed a phase 2b efficacy study of CYD-TDV in a single center in Thailand, were included in the study. Genotyping of HLA class I and II loci was performed by next-generation sequencing from DNA obtained from 335 saliva samples. Dengue neutralizing antibody titers (NAb) were assessed as a correlate of risk and protection. Regression analyses were used to assess associations between HLA alleles and NAb responses, vaccine efficacy, and dengue outcomes. Month 13 NAb log geometric mean titers (GMTs) were associated with decreased risk of VCD. In the vaccine group, HLA-DRB1*11 was significantly associated with higher NAb log GMT levels (beta: 0.76; p = 0.002, q = 0.13). Additionally, in the absence of vaccination, HLA associations were observed between the presence of DPB1*03:01 and increased NAb log GMT levels (beta: 1.24; p = 0.005, q = 0.17), and between DPB1*05:01 and reduced NAb log GMT levels (beta: -1.1; p = 0.001, q = 0.07). This study suggests associations of HLA alleles with NAb titers in the context of dengue outcomes. This study was registered with clinicaltrials.gov: NCT02827162.

Monocyte-derived transcriptome signature indicates antibody-dependent cellular phagocytosis as a potential mechanism of vaccine-induced protection against HIV-1.

A gene signature was previously found to be correlated with mosaic adenovirus 26 vaccine protection in simian immunodeficiency virus and simian-human immunodeficiency virus challenge models in non-human primates. In this report, we investigated the presence of this signature as a correlate of reduced risk in human clinical trials and potential mechanisms of protection. The absence of this gene signature in the DNA/rAd5 human vaccine trial, which did not show efficacy, strengthens our hypothesis that this signature is only enriched in studies that demonstrated protection. This gene signature was enriched in the partially effective RV144 human trial that administered the ALVAC/protein vaccine, and we find that the signature associates with both decreased risk of HIV-1 acquisition and increased vaccine efficacy (VE). Total RNA-seq in a clinical trial that used the same vaccine regimen as the RV144 HIV vaccine implicated antibody-dependent cellular phagocytosis (ADCP) as a potential mechanism of vaccine protection. CITE-seq profiling of 53 surface markers and transcriptomes of 53,777 single cells from the same trial showed that genes in this signature were primarily expressed in cells belonging to the myeloid lineage, including monocytes, which are major effector cells for ADCP. The consistent association of this transcriptome signature with VE represents a tool both to identify potential mechanisms, as with ADCP here, and to screen novel approaches to accelerate the development of new vaccine candidates.

B cell engagement with HIV-1 founder virus envelope predicts development of broadly neutralizing antibodies.

Determining which immunological mechanisms contribute to the development of broad neutralizing antibodies (bNAbs) during HIV-1 infection is a major goal to inform vaccine design. Using samples from a longitudinal HIV-1 acute infection cohort, we found key B cell determinants within the first 14-43 days of viremia that predict the development of bNAbs years later. Individuals who develop neutralization breadth had significantly higher B cell engagement with the autologous founder HIV envelope (Env) within 1 month of initial viremia. A higher frequency of founder-Env-specific naive B cells was associated with increased B cell activation and differentiation and predictive of bNAb development. These data demonstrate that the initial B cell interaction with the founder HIV Env is important for the development of broadly neutralizing antibodies and provide evidence that events within HIV acute infection lead to downstream functional outcomes.

Next-generation sequencing of 11 HLA loci in a large dengue vaccine cohort from the Philippines.

HLA genotyping by next-generation sequencing (NGS) has evolved with significant advancements in the last decade. Here we describe full-length HLA genotyping of 11 loci in 612 individuals comprising a dengue vaccine cohort from Cebu province in the Philippines. The multi-locus individual tagging NGS (MIT-NGS) method that we developed initially for genotyping 4-6 loci in one MiSeq run was expanded to 11 loci including HLA-A, B, C, DPA1, DPB1, DQA1, DQB1, DRB1, and DRB3/4/5. This change did not affect the overall coverage or depth of the sequencing reads. HLA alleles with frequencies greater than 10% were A*11:01:01, A*24:02:01, A*24:07:01, A*34:01:01, B*38:02:01, B*15:35, B*35:05:01, C*07:02:01, C*04:01:01, DPA1*02:02:02, DPB1*05:01:01, DPB1*01:01:01, DQA1*01:02:01, DQA1*06:01:01, DQB1*05:02:01, DQB1*03:01:01, DRB1*15:02:01, DRB1*12:02:01, DRB3*03:01:03, DRB4*01:03:01, and DRB5*01:01:01. Improvements in sequencing library preparation provide uniform and even coverage across all exons and introns. This has led to a marked reduction in allele imbalance and dropout. Furthermore, including more loci, such as DRB3/4/5, decreases cross-mapping and incorrect allele assignment at the DRB1 locus. The increased number of loci sequenced for each sample does not reduce the number of samples that can be multiplexed on a single MiSeq run and is therefore more cost-efficient. We believe that such improvements will help HLA genotyping by NGS to gain momentum over other conventional methods by increasing confidence in the calls.

Potent Zika and dengue cross-neutralizing antibodies induced by Zika vaccination in a dengue-experienced donor.

Zika virus (ZIKV) has caused significant disease, with widespread cases of neurological pathology and congenital neurologic defects. Rapid vaccine development has led to a number of candidates capable of eliciting potent ZIKV-neutralizing antibodies (reviewed in refs. 1-3). Despite advances in vaccine development, it remains unclear how ZIKV vaccination affects immune responses in humans with prior flavivirus immunity. Here we show that a single-dose immunization of ZIKV purified inactivated vaccine (ZPIV)4-7 in a dengue virus (DENV)-experienced human elicited potent cross-neutralizing antibodies to both ZIKV and DENV. Using a unique ZIKV virion-based sorting strategy, we isolated and characterized multiple antibodies, including one termed MZ4, which targets a novel site of vulnerability centered on the Envelope (E) domain I/III linker region and protects mice from viremia and viral dissemination following ZIKV or DENV-2 challenge. These data demonstrate that Zika vaccination in a DENV-experienced individual can boost pre-existing flavivirus immunity and elicit protective responses against both ZIKV and DENV. ZPIV vaccination in Puerto Rican individuals with prior flavivirus experience yielded similar cross-neutralizing potency after a single vaccination, highlighting the potential benefit of ZIKV vaccination in flavivirus-endemic areas.

A vaccine-induced gene expression signature correlates with protection against SIV and HIV in multiple trials.

Current HIV vaccines are only partially efficacious, presenting an opportunity to identify correlates of protection and, thereby, potential insight into mechanisms that prevent HIV acquisition. Two independent preclinical challenge studies in nonhuman primates (NHPs) previously showed partial efficacy of a mosaic adenovirus 26 (Ad26)-based HIV-1 vaccine candidate. To investigate the basis of this protection, we performed whole transcriptomics profiling by RNA sequencing (RNA-seq) in sorted lymphocytes from peripheral blood samples taken during these studies at different time points after vaccination but before challenge. We observed a transcriptional signature in B cells that associated with protection from acquisition of simian immunodeficiency virus (SIV) or the simian-human immunodeficiency virus (SHIV) in both studies. Strong antibody responses were elicited, and genes from the signature for which expression was enriched specifically associated with higher magnitude of functional antibody responses. The same gene expression signature also associated with protection in RV144 in the only human HIV vaccine trial to date that has shown efficacy and in two additional NHP studies that evaluated similar canarypox-based vaccine regimens. A composite gene expression score derived from the gene signature was one of the top-ranked correlates of protection in the NHP vaccine studies. This study aims to bridge preclinical and clinical data with the identification of a gene signature in B cells that is associated with protection from SIV and HIV infection by providing a new approach for evaluating future vaccine candidates.

Full-length next-generation sequencing of HLA class I and II genes in a cohort from Thailand.

The human leukocyte antigen (HLA) genes are highly variable and are known to play an important role in disease outcomes, including infectious diseases. Prior knowledge of HLA polymorphisms in a population usually forms the basis for an effective case-control study design. As a prelude to future disease association analyses, we report HLA class I and II diversity in 334 unrelated donors from a Dengue vaccine efficacy trial conducted in Thailand. Long-range PCR amplification of six HLA loci was performed on DNA extracted from saliva samples. HLA-A, -B, -C, -DPB1, -DQB1 and -DRB1 were genotyped using a next-generation sequencing method presented at the 17th International HLA and Immunogenetics Workshop. In total, we identified 201 HLA alleles, including 35 HLA-A, 57 HLA-B, 28 HLA-C, 24 HLA-DPB1, 21 HLA-DQB1 and 36 HLA-DRB1 alleles. Very common HLA alleles with frequencies greater than 10 percent were A∗11:01:01, A∗33:03:01, A∗24:02:01, B∗46:01:01, C∗07:02:01, C∗01:02:01, C∗08:01:01, DPB1∗05:01:01, DPB1∗13:01:01, DPB1∗04:01:01, DPB1∗02:01:02, DQB1∗03:01:01, DQB1∗05:02:01, DQB1∗03:03:02, DRB1∗12:02:01, DRB1∗09:01:02, and DRB1∗15:02:01. A novel HLA allele, B∗15:450, had a non-synonymous substitution and occurred in more than one donor. Population-based full-length NGS HLA typing is more conclusive and provides a sound foundation for exploring disease association in a given population.

Variations in HLA-B cell surface expression, half-life and extracellular antigen receptivity.

The highly polymorphic human leukocyte antigen (HLA) class I molecules present peptide antigens to CD8+ T cells, inducing immunity against infections and cancers. Quality control mediated by peptide loading complex (PLC) components is expected to ensure the cell surface expression of stable peptide-HLA class I complexes. This is exemplified by HLA-B*08:01 in primary human lymphocytes, with both expression level and half-life at the high end of the measured HLA-B expression and stability hierarchies. Conversely, low expression on lymphocytes is measured for three HLA-B allotypes that bind peptides with proline at position 2, which are disfavored by the transporter associated with antigen processing. Surprisingly, these lymphocyte-specific expression and stability differences become reversed or altered in monocytes, which display larger intracellular pools of HLA class I than lymphocytes. Together, the findings indicate that allele and cell-dependent variations in antigen acquisition pathways influence HLA-B surface expression levels, half-lives and receptivity to exogenous antigens.

High-Throughput Contiguous Full-Length Next-Generation Sequencing of HLA Class I and II Genes from 96 Donors in a Single MiSeq Run.

The human leukocyte antigen (HLA) genes regulate and drive the immune system, and are among the most polymorphic loci in the human genome. HLA diversity is known to play an important role in transplantation and disease association studies. There are multiple approaches to DNA-based HLA genotyping and recent advances in next-generation sequencing (NGS) technologies have facilitated the development of whole gene sequencing methods. We describe an accurate, efficient, scalable, and cost-effective approach to contiguously amplify and sequence full-length genes of six HLA class I and II loci from 96 individuals on a single Illumina MiSeq run.

High-throughput multiplex HLA genotyping by next-generation sequencing using multi-locus individual tagging.

BACKGROUND: Unambiguous human leukocyte antigen (HLA) typing is important in transplant matching and disease association studies. High-resolution HLA typing that is not restricted to the peptide-binding region can decrease HLA allele ambiguities. Cost and technology constraints have hampered high-throughput and efficient high resolution unambiguous HLA typing. We have developed a method for HLA genotyping that preserves the very high-resolution that can be obtained by next-generation sequencing (NGS) but also achieves substantially increased efficiency. Unambiguous HLA-A, B, C and DRB1 genotypes can be determined for 96 individuals in a single run of the Illumina MiSeq. RESULTS: Long-range amplification of full-length HLA genes from four loci was performed in separate polymerase chain reactions (PCR) using primers and PCR conditions that were optimized to reduce co-amplification of other HLA loci. Amplicons from the four HLA loci of each individual were then pooled and subjected to enzymatic library generation. All four loci of an individual were then tagged with one unique index combination. This multi-locus individual tagging (MIT) method combined with NGS enabled the four loci of 96 individuals to be analyzed in a single 500 cycle sequencing paired-end run of the Illumina-MiSeq. The MIT-NGS method generated sequence reads from the four loci were then discriminated using commercially available NGS HLA typing software. Comparison of the MIT-NGS with Sanger sequence-based HLA typing methods showed that all the ambiguities and discordances between the two methods were due to the accuracy of the MIT-NGS method. CONCLUSIONS: The MIT-NGS method enabled accurate, robust and cost effective simultaneous analyses of four HLA loci per sample and produced 6 or 8-digit high-resolution unambiguous phased HLA typing data from 96 individuals in a single NGS run.

Floral Transcriptomes in Woodland Strawberry Uncover Developing Receptacle and Anther Gene Networks.

Flowers are reproductive organs and precursors to fruits and seeds. While the basic tenets of the ABCE model of flower development are conserved in angiosperms, different flowering plants exhibit different and sometimes unique characteristics. A distinct feature of strawberry (Fragaria spp.) flowers is the development of several hundreds of individual apocarpous (unfused) carpels. These individual carpels are arranged in a spiral pattern on the subtending stem tip, the receptacle. Therefore, the receptacle is an integral part of the strawberry flower and is of significant agronomic importance, being the precursor to strawberry fruit. Taking advantage of next-generation sequencing and laser capture microdissection, we generated different tissue- and stage-transcriptomic profiling of woodland strawberry (Fragaria vesca) flower development. Using pairwise comparisons and weighted gene coexpression network analysis, we identified modules of coexpressed genes and hub genes of tissue-specific networks. Of particular importance is the discovery of a developing receptacle-specific module exhibiting similar molecular features to those of young floral meristems. The strawberry homologs of a number of meristem regulators, including LOST MERISTEM and WUSCHEL, are identified as hub genes operating in the developing receptacle network. Furthermore, almost 25% of the F-box genes in the genome are transiently induced in developing anthers at the meiosis stage, indicating active protein degradation. Together, this work provides important insights into the molecular networks underlying strawberry's unique reproductive developmental processes. This extensive floral transcriptome data set is publicly available and can be readily queried at the project Web site, serving as an important genomic resource for the plant biology research community.

HLA class I, KIR, and genome-wide SNP diversity in the RV144 Thai phase 3 HIV vaccine clinical trial.

RV144 is the first phase 3 HIV vaccine clinical trial to demonstrate efficacy. This study consisted of more than 8,000 individuals in each arm of the trial, representing the four major regions of Thailand. Human leukocyte antigen (HLA) class I and killer cell immunoglobulin-like receptor (KIR) genes, as well as 96 genome-wide ancestry informative markers (AIMs) were genotyped in 450 placebo HIV-1-uninfected individuals to identify the immunogenetic diversity and population structure of this cohort. High-resolution genotyping identified the common HLA alleles as A*02:03, A*02:07, A*11:01, A*24:02, A*24:07, A*33:03, B*13:01, B*15:02, B*18:01, B*40:01, B*44:03, B*46:01, B*58:01, C*01:02, C*03:02, C*03:04, C*07:01, C*07:02, C*07:04, and C*08:01. The most frequent three-loci haplotype was B*46:01-C*01:02-A*02:07. Framework genes KIR2DL4, 3DL2, and 3DL3 were present in all samples, and KIR2DL1, 2DL3, 3DL1, 2DS4, and 2DP1 occurred at frequencies greater than 90 %. The combined HLA and KIR profile suggests admixture with neighboring Asian populations. Principal component and correspondence analyses comparing the RV144 samples to the phase 3 International HapMap Project (HapMap3) populations using AIMs corroborated these findings. Structure analyses identified a distinct profile in the Thai population that did not match the Asian or other HapMap3 samples. This shows genetic variability unique to Thais in RV144, making it essential to take into account population stratification while performing genetic association studies. The overall analyses from all three genetic markers indicate that the RV144 samples are representative of the Thai population. This will inform subsequent host genetic analyses in the RV144 cohort and provide insight for future genetic association studies in the Thai population.

SGR: an online genomic resource for the woodland strawberry.

BACKGROUND: Fragaria vesca, a diploid strawberry species commonly known as the alpine or woodland strawberry, is a versatile experimental plant system and an emerging model for the Rosaceae family. An ancestral F. vesca genome contributed to the genome of the octoploid dessert strawberry (F. ×ananassa), and the extant genome exhibits synteny with other commercially important members of the Rosaceae family such as apple and peach. To provide a molecular description of floral organ and fruit development at the resolution of specific tissues and cell types, RNAs from flowers and early developmental stage fruit tissues of the inbred F. vesca line YW5AF7 were extracted and the resulting cDNA libraries sequenced using an Illumina HiSeq2000. To enable easy access as well as mining of this two-dimensional (stage and tissue) transcriptome dataset, a web-based database, the Strawberry Genomic Resource (SGR), was developed. DESCRIPTION: SGR is a web accessible database that contains sample description, sample statistics, gene annotation, and gene expression analysis. This information can be accessed publicly from a web-based interface at http://bioinformatics.towson.edu/strawberry/Default.aspx. The SGR website provides user friendly search and browse capabilities for all the data stored in the database. Users are able to search for genes using a gene ID or description or obtain differentially expressed genes by entering different comparison parameters. Search results can be downloaded in a tabular format compatible with Microsoft excel application. Aligned reads to individual genes and exon/intron structures are displayed using the genome browser, facilitating gene re-annotation by individual users. CONCLUSIONS: The SGR database was developed to facilitate dissemination and data mining of extensive floral and fruit transcriptome data in the woodland strawberry. It enables users to mine the data in different ways to study different pathways or biological processes during reproductive development.

Genome-scale transcriptomic insights into early-stage fruit development in woodland strawberry Fragaria vesca.

Fragaria vesca, a diploid woodland strawberry with a small and sequenced genome, is an excellent model for studying fruit development. The strawberry fruit is unique in that the edible flesh is actually enlarged receptacle tissue. The true fruit are the numerous dry achenes dotting the receptacle's surface. Auxin produced from the achene is essential for the receptacle fruit set, a paradigm for studying crosstalk between hormone signaling and development. To investigate the molecular mechanism underlying strawberry fruit set, next-generation sequencing was employed to profile early-stage fruit development with five fruit tissue types and five developmental stages from floral anthesis to enlarged fruits. This two-dimensional data set provides a systems-level view of molecular events with precise spatial and temporal resolution. The data suggest that the endosperm and seed coat may play a more prominent role than the embryo in auxin and gibberellin biosynthesis for fruit set. A model is proposed to illustrate how hormonal signals produced in the endosperm and seed coat coordinate seed, ovary wall, and receptacle fruit development. The comprehensive fruit transcriptome data set provides a wealth of genomic resources for the strawberry and Rosaceae communities as well as unprecedented molecular insight into fruit set and early stage fruit development.

Flower and early fruit development in a diploid strawberry, Fragaria vesca.

The diploid woodland strawberry, Fragaria vesca, is being recognized as a model for the more complex octoploid commercial strawberry, Fragaria × ananassa. F. vesca exhibits a short seed to seed cycle, can be easily transformed by Agrobacteria, and a draft genome sequence has been published. These features, together with its similar flower structure, potentially make F. vesca a good model for studying the flower development of other members of the Rosaceae family, which contains many economically important fruit trees and ornamental plants. To propel F. vesca's role in genetic and genomic research and to facilitate the study of its reproductive development, we have investigated in detail F. vesca flower and early fruit development using a seventh generation inbred diploid line, Yellow Wonder 5AF7. We present here standardized developmental staging and detailed descriptions of morphological changes associated with flower and early fruit development based on images of hand dissected flowers, histological sections, and scanning electron microscopy. In situ hybridization with the F. vesca AGAMOUS homolog, FvAG, showed expression in young stamen and carpel primordia. This work lays the essential groundwork and standardization for future molecular, genetic, and genomic studies of F. vesca.