Small Extracellular Vesicles Promote Stiffness-mediated Metastasis.

Tissue stiffness is a critical prognostic factor in breast cancer and is associated with metastatic progression. Here we show an alternative and complementary hypothesis of tumor progression whereby physiologic matrix stiffness affects the quantity and protein cargo of small extracellular vesicles (EV) produced by cancer cells, which in turn aid cancer cell dissemination. Primary patient breast tissue released by cancer cells on matrices that model human breast tumors (25 kPa; stiff EVs) feature increased adhesion molecule presentation (ITGα2ß1, ITGα6ß4, ITGα6ß1, CD44) compared with EVs from softer normal tissue (0.5 kPa; soft EVs), which facilitates their binding to extracellular matrix proteins including collagen IV, and a 3-fold increase in homing ability to distant organs in mice. In a zebrafish xenograft model, stiff EVs aid cancer cell dissemination. Moreover, normal, resident lung fibroblasts treated with stiff and soft EVs change their gene expression profiles to adopt a cancer-associated fibroblast phenotype. These findings show that EV quantity, cargo, and function depend heavily on the mechanical properties of the extracellular microenvironment. SIGNIFICANCE: Here we show that the quantity, cargo, and function of breast cancer-derived EVs vary with mechanical properties of the extracellular microenvironment.

Over-expression of Medicago Acyl-CoA-binding 2 genes enhance salt and drought tolerance in Arabidopsis.

Acyl-CoA-binding proteins (ACBPs) are mainly involved in acyl-CoA ester binding and trafficking in eukaryotic cells, and they function in lipid metabolism, membrane biosynthesis, cellular signaling, stress response, disease resistance, and other biological activities in plants. However, the roles of ACBP family members in Medicago remain unclear. In this study, a total of eight ACBP genes were identified in the genome of Medicago truncatula and Medicago sativa, and they were clustered into four sub-families (Class I-IV). Many cis-acting elements related to abiotic response were identified in the promoter region of these ACBP genes, in particular light-responsive elements. These ACBP genes exhibited distinct expression pattern in various tissues, and the expression level of MtACBP1/MsACBP1 and MtACBP2/MsACBP2 gene pairs were significantly increased under NaCl treatment. Subcellular localization analysis showed that MtACBP1/MsACBP1 and MtACBP2/MsACBP2 were localized in the endoplasmic reticulum of tobacco epidermal cells. Arabidopsis seedlings over-expressing MtACBP2/MsACBP2 displayed increased root length than the wild type under short light, Cu2+, ABA, PEG, and NaCl treatments. Over-expression of MtACBP2/MsACBP2 also significantly enhanced Arabidopsis tolerance under NaCl and PEG treatments in mature plants. Collectively, our study identified salt and drought responsive ACBP genes in Medicago and verified their functions in increasing resistance against salt and drought stresses.

A 3D in vitro assay to study combined immune cell infiltration and cytotoxicity.

Immune cell-mediated killing of cancer cells in a solid tumor is prefaced by a multi-step infiltration cascade of invasion, directed migration, and cytotoxic activities. In particular, immune cells must invade and migrate through a series of different extracellular matrix (ECM) boundaries and domains before reaching and killing their target tumor cells. These infiltration events are a central challenge to the clinical success of CAR T cells against solid tumors. The current standard in vitro cell killing assays measure cell cytotoxicity in an obstacle-free, two-dimensional (2D) microenvironment, which precludes the study of 3D immune cell-ECM interactions. Here, we present a 3D combined infiltration/cytotoxicity assay based on an oil-in-water microtechnology. This assay measures stromal invasion following extravasation, migration through the stromal matrix, and invasion of the solid tumor in addition to cell killing. We compare this 3D cytotoxicity assay to the benchmark 2D assay through tumor assembloid cocultures with immune cells and engineered immune cells. This assay is amenable to an array of imaging techniques, which allows direct observation and quantification of each stage of infiltration in different immune and oncological contexts. We establish the 3D infiltration/cytotoxicity assay as an important tool for the mechanistic study of immune cell interactions with the tumor microenvironment.

Destruction of vascular endothelial glycocalyx during formation of pre-metastatic niches.

A special microenvironment called the "pre-metastatic niche" is thought to help primary tumor cells migrate to new tissues and invade them, in part because the normal barrier function of the vascular endothelium is compromised. While the primary tumor itself can promote the creation of such niches by secreting pro-metastatic factors, the underlying molecular mechanisms are still poorly understood. Here, we show that the injection of primary tumor-secreted pro-metastatic factors from B16F10 melanoma or 4T1 breast cancer cells into healthy mice can induce the destruction of the vascular endothelial glycocalyx, which is a polysaccharide coating on the vascular endothelial lumen that normally inhibits tumor cell passage into and out of the circulation. However, when human umbilical vein endothelial cultures were treated in vitro with these secreted pro-metastatic factors, no significant destruction of the glycocalyx was observed, implying that this destruction requires a complex in vivo microenvironment. The tissue section analysis revealed that secreted pro-metastatic factors could clearly upregulate macrophage-related molecules such as CD11b and tumor necrosis factor-α (TNF-α) in the heart, liver, spleen, lung, and kidney, which is associated with the upregulation and activation of heparanase. In addition, macrophage depletion significantly attenuated the degradation of the vascular endothelial glycocalyx induced by secreted pro-metastatic factors. This indicates that the secreted pro-metastatic factors that destroy the vascular endothelial glycocalyx rely primarily on macrophages. Our findings suggest that the formation of pre-metastatic niches involves degradation of the vascular endothelial glycocalyx, which may hence be a useful target for developing therapies to inhibit cancer metastasis.

Medicago truncatula ß-glucosidase 17 contributes to drought and salt tolerance through antioxidant flavonoid accumulation.

Flavonoids are usually present in forms of glucosides in plants, which could be catabolized by ß-glucosidase (BGLU) to form their corresponding flavonoid aglycones. In this study, we isolated three abiotic-responsive BGLU genes (MtBGLU17, MtBGLU21 and MtBGLU22) from Medicago truncatula, and found only the recombinant MtBGLU17 protein could catalyse the hydrolysis of flavonoid glycosides. The recombinant MtBGLU17 protein is active towards a variety of flavonoid glucosides, including glucosides of flavones (apigenin and luteolin), flavonols (kaempferol and quercetin), isoflavones (genistein and daidzein) and flavanone (naringenin). In particular, the recombinant MtBGLU17 protein preferentially hydrolyses flavonoid-7-O-glucosides over their corresponding 3-O-glucosides. The content of luteoin-7-O-glucoside was reduced in the MtBGLU17 overexpression plants but increased in the Tnt-1 insertional mutant lines, whereas luteoin content was increased in the MtBGLU17 overexpression plants but reduced in the Tnt-1 insertional mutant lines. Under drought and salt (NaCl) treatment, the MtBGLU17 overexpression lines showed relatively higher DPPH content, and higher CAT and SOD activity than the wild type control. These results indicated that overexpression lines of MtBGLU17 possess higher antioxidant activity and thus confer drought and salt tolerance, implying MtBGLU17 could be potentially used as a candidate gene to improve plant abiotic stress tolerance.

Overexpression of MsNIP2 improves salinity tolerance in Medicago sativa.

Alfalfa (Medicago sativa) is one of the most widely cultivated forage crops in the world. However, alfalfa yield and quality are adversely affected by salinity stress. Nodulin 26-like intrinsic proteins (NIPs) play essential roles in water and small molecules transport and response to salt stress. Here, we isolated a salt stress responsive MsNIP2 gene and demonstrated its functions by overexpression in alfalfa. The open reading frame of MsNIP2 is 816 bp in length, and it encodes 272 amino acids. It has six transmembrane domains and two NPA motifs. MsNIP2 showed high identity to other known NIP proteins, and its tertiary model was similar to the crystal structure of OsNIP2-1 (7cjs) tetramer. Subcellular localization analysis showed that MsNIP2 protein fused with green fluorescent protein (GFP) was localized to the plasma membrane. Transgenic alfalfa lines overexpressing MsNIP2 showed significantly higher height and branch number compared with the non-transgenic control. The POD and CAT activity of the transgenic alfalfa lines was significantly increased and their MDA content was notably reduced compared with the control group under the treatment of NaCl. The transgenic lines showed higher capability in scavenging oxygen radicals with lighter NBT staining than the control under salt stress. The transgenic lines showed relative lower water loss rate and electrolyte leakage, but relatively higher Na+ content than the control line under salt stress. The relative expression levels of abiotic-stress-related genes (MsHSP23, MsCOR47, MsATPase, and MsRD2) in three transgenic lines were compared with the control, among them, only the expression of MsCOR47 was up-regulated. Consequently, this study offers a novel perspective for exploring the function of MsNIP2 in improving salt tolerance of alfalfa.

Engineering self-propelled tumor-infiltrating CAR T cells using synthetic velocity receptors.

Chimeric antigen receptor (CAR) T cells express antigen-specific synthetic receptors, which upon binding to cancer cells, elicit T cell anti-tumor responses. CAR T cell therapy has enjoyed success in the clinic for hematological cancer indications, giving rise to decade-long remissions in some cases. However, CAR T therapy for patients with solid tumors has not seen similar success. Solid tumors constitute 90% of adult human cancers, representing an enormous unmet clinical need. Current approaches do not solve the central problem of limited ability of therapeutic cells to migrate through the stromal matrix. We discover that T cells at low and high density display low- and high-migration phenotypes, respectively. The highly migratory phenotype is mediated by a paracrine pathway from a group of self-produced cytokines that include IL5, TNFα, IFNγ, and IL8. We exploit this finding to "lock-in" a highly migratory phenotype by developing and expressing receptors, which we call velocity receptors (VRs). VRs target these cytokines and signal through these cytokines' cognate receptors to increase T cell motility and infiltrate lung, ovarian, and pancreatic tumors in large numbers and at doses for which control CAR T cells remain confined to the tumor periphery. In contrast to CAR therapy alone, VR-CAR T cells significantly attenuate tumor growth and extend overall survival. This work suggests that approaches to the design of immune cell receptors that focus on migration signaling will help current and future CAR cellular therapies to infiltrate deep into solid tumors.

Small extracellular vesicles promote stiffness-mediated metastasis.

Tissue stiffness is a critical prognostic factor in breast cancer and is associated with metastatic progression. Here we show an alternative and complementary hypothesis of tumor progression whereby physiological matrix stiffness affects the quantity and protein cargo of small EVs produced by cancer cells, which in turn drive their metastasis. Primary patient breast tissue produces significantly more EVs from stiff tumor tissue than soft tumor adjacent tissue. EVs released by cancer cells on matrices that model human breast tumors (25 kPa; stiff EVs) feature increased adhesion molecule presentation (ITGα 2 ß 1 , ITGα 6 ß 4 , ITGα 6 ß 1 , CD44) compared to EVs from softer normal tissue (0.5 kPa; soft EVs), which facilitates their binding to extracellular matrix (ECM) protein collagen IV, and a 3-fold increase in homing ability to distant organs in mice. In a zebrafish xenograft model, stiff EVs aid cancer cell dissemination through enhanced chemotaxis. Moreover, normal, resident lung fibroblasts treated with stiff and soft EVs change their gene expression profiles to adopt a cancer associated fibroblast (CAF) phenotype. These findings show that EV quantity, cargo, and function depend heavily on the mechanical properties of the extracellular microenvironment.

Mechano-induced homotypic patterned domain formation by monocytes.

Matrix stiffness and corresponding mechano-signaling play indispensable roles in cellular phenotypes and functions. How tissue stiffness influences the behavior of monocytes, a major circulating leukocyte of the innate system, and how it may promote the emergence of collective cell behavior is less understood. Here, using tunable collagen-coated hydrogels of physiological stiffness, we show that human primary monocytes undergo a dynamic local phase separation to form highly regular, reversible, multicellular, multi-layered domains on soft matrix. Local activation of the ß2 integrin initiates inter-cellular adhesion, while global soluble inhibitory factors maintain the steady state domain pattern over days. Patterned domain formation generated by monocytes is unique among other key immune cells, including macrophages, B cells, T cells, and NK cells. While inhibiting their phagocytic capability, domain formation promotes monocytes' survival. We develop a computational model based on the Cahn-Hilliard equation of phase separation, combined with a Turing mechanism of local activation and global inhibition suggested by our experiments, and provides experimentally validated predictions of the role of seeding density and both chemotactic and random cell migration on domain pattern formation. This work reveals that, unlike active matters, cells can generate complex cell phases by exploiting their mechanosensing abilities and combined short-range interactions and long-range signals to enhance their survival.

Mechano-induced homotypic patterned domain formation by monocytes.

Matrix stiffness and corresponding mechano-signaling play indispensable roles in cellular phenotypes and functions. How tissue stiffness influences the behavior of monocytes, a major circulating leukocyte of the innate system, and how it may promote the emergence of collective cell behavior is less understood. Here, using tunable collagen-coated hydrogels of physiological stiffness, we show that human primary monocytes undergo a dynamic local phase separation to form highly patterned multicellular multi-layered domains on soft matrix. Local activation of the ß2 integrin initiates inter-cellular adhesion, while global soluble inhibitory factors maintain the steady-state domain pattern over days. Patterned domain formation generated by monocytes is unique among other key immune cells, including macrophages, B cells, T cells, and NK cells. While inhibiting their phagocytic capability, domain formation promotes monocytes' survival. We develop a computational model based on the Cahn-Hilliard equation, which includes combined local activation and global inhibition mechanisms of intercellular adhesion suggested by our experiments, and provides experimentally validated predictions of the role of seeding density and both chemotactic and random cell migration on pattern formation.

Genome-wide characterization and expression analysis of the HAK gene family in response to abiotic stresses in Medicago.

The high-affinity K+ transporter (HAK) family plays a vital role in K+ uptake and transport as well as in salt and drought stress responses. In the present study, we identified 22 HAK genes in each Medicago truncatula and Medicago sativa genome. Phylogenetic analysis suggested that these HAK proteins could be divided into four clades, and the members of the same subgroup share similar gene structure and conserved motifs. Many cis-acting elements related with defense and stress were found in their promoter region. In addition, gene expression profiles analyzed with genechip and transcriptome data showed that these HAK genes exhibited distinct expression pattern in different tissues, and in response to salt and drought treatments. Furthermore, co-expression analysis showed that 6 homologous HAK hub gene pairs involved in direct network interactions. RT-qPCR verified that the expression level of six HAK gene pairs was induced by NaCl and mannitol treatment to different extents. In particular, MtHK2/7/12 from M. truncatula and MsHAK2/6/7 from M. sativa were highly induced. The expression level of MsHAK1/2/11 determined by RT-qPCR showed significantly positive correlation with transcriptome data. In conclusion, our study shows that HAK genes play a key role in response to various abiotic stresses in Medicago, and the highly inducible candidate HAK genes could be used for further functional studies and molecular breeding in Medicago.

Transcriptome and functional analyses reveal ERF053 from Medicago falcata as key regulator in drought resistances.

Medicago falcata L. is an important legume forage grass with strong drought resistant, which could be utilized as an important gene pool in molecular breed of forage grass. In this study, M. falcata seedlings were treated with 400 mM mannitol to simulate drought stress, and the morphological and physiological changes were investigated, as well as the transcriptome changes of M. falcata seedlings at different treatment time points (0 h, 2 h, 6 h, 12 h, 24 h, 36 h and 48 h). Transcriptome analyses revealed four modules were closely related with drought response in M. falcata by WGCNA analysis, and four ERF transcription factor genes related with drought stress were identified (MfERF053, MfERF9, MfERF034 and MfRAP2.1). Among them, MfERF053 was highly expressed in roots, and MfERF053 protein showed transcriptional activation activity by transient expression in tobacco leaves. Overexpression of MfERF053 in Arabidopsis improved root growth, number of lateral roots and fresh weight under drought, salt stress and exogenous ABA treatments. Transgenic Arabidopsis over-expressing MfERF053 gene grew significantly better than the wild type under both drought stress and salt stress when grown in soil. Taken together, our strategy with transcriptome combined WGCNA analyses identified key transcription factor genes from M. falcata, and the selected MfERF053 gene was verified to be able to enhance drought and salt resistance when over-expressed in Arabidopsis.

Genome-Wide Identification and Characterization of Growth Regulatory Factor Family Genes in Medicago.

Growth Regulatory Factors (GRF) are plant-specific transcription factors that play critical roles in plant growth and development as well as plant tolerance against stress. In this study, a total of 16 GRF genes were identified from the genomes of Medicago truncatula and Medicago sativa. Multiple sequence alignment analysis showed that all these members contain conserved QLQ and WRC domains. Phylogenetic analysis suggested that these GRF proteins could be classified into five clusters. The GRF genes showed similar exon-intron organizations and similar architectures in their conserved motifs. Many stress-related cis-acting elements were found in their promoter region, and most of them were related to drought and defense response. In addition, analyses on microarray and transcriptome data indicated that these GRF genes exhibited distinct expression patterns in various tissues or in response to drought and salt treatments. In particular, qPCR results showed that the expression levels of gene pairs MtGRF2-MsGRF2 and MtGRF6-MsGRF6 were significantly increased under NaCl and mannitol treatments, indicating that they are most likely involved in salt and drought stress tolerance. Collectively, our study is valuable for further investigation on the function of GRF genes in Medicago and for the exploration of GRF genes in the molecular breeding of highly resistant M. sativa.

Identification and Characterization of Abiotic Stress-Responsive NF-YB Family Genes in Medicago.

Nuclear factor YB (NF-YB) are plant-specific transcription factors that play a critical regulatory role in plant growth and development as well as in plant resistance against various stresses. In this study, a total of 49 NF-YB genes were identified from the genomes of Medicago truncatula and Medicago sativa. Multiple sequence alignment analysis showed that all of these NF-YB members contain DNA binding domain, NF-YA interaction domain and NF-YC interaction domain. Phylogenetic analysis suggested that these NF-YB proteins could be classified into five distinct clusters. We also analyzed the exon-intron organizations and conserved motifs of these NF-YB genes and their deduced proteins. We also found many stress-related cis-acting elements in their promoter region. In addition, analyses on genechip for M. truncatula and transcriptome data for M. sativa indicated that these NF-YB genes exhibited a distinct expression pattern in various tissues; many of these could be induced by drought and/or salt treatments. In particular, RT-qPCR analysis revealed that the expression levels of gene pairs MsNF-YB27/MtNF-YB15 and MsNF-YB28/MtNF-YB16 were significantly up-regulated under NaCl and mannitol treatments, indicating that they are most likely involved in salt and drought stress response. Taken together, our study on NF-YB family genes in Medicago is valuable for their functional characterization, as well as for the application of NF-YB genes in genetic breeding for high-yield and high-resistance alfalfa.

Fish female-biased gene cyp19a1a leads to female antiviral response attenuation between sexes by autophagic degradation of MITA.

From insects to mammals, both innate and adaptive immune response are usually higher in females than in males, with the sex chromosome and hormonal differences considered the main reasons. Here, we report that zebrafish cyp19a1a (cytochrome P450, family 19, subfamily A, polypeptide 1a), an autosomal gene with female-biased expression, causes female fish to exhibit a lower antiviral response. First, we successfully constructed an infection model by intraperitoneal injection of spring viremia of carp virus (SVCV) into zebrafish (Danio rerio) and Carassius auratus herpesvirus (CaHV) in gibel carp (Carassius gibelio). Specifically, female fish were more vulnerable to viral infection than males, accompanied by a significantly weaker interferon (IFN) expression. After screening several candidates, cyp19a1a, which was highly expressed in female fish tissues, was selected for further analysis. The IFN expression and antiviral response were significantly higher in cyp19a1a-/- than in cyp19a1a+/+. Further investigation of the molecular mechanism revealed that Cyp19a1a targets mediator of IRF3 activation (MITA) for autophagic degradation. Interestingly, in the absence of MITA, Cyp19a1a alone could not elicit an autophagic response. Furthermore, the autophagy factor ATG14 (autophagy-related 14) was found interacted with Cyp19a1a to either promote or attenuate Cyp19a1a-mediated MITA degradation by either being overexpressed or knocked down, respectively. At the cellular level, both the normal and MITA-enhanced cellular antiviral responses were diminished by Cyp19a1a. These findings demonstrated a sex difference in the antiviral response based on a regulation mechanism controlled by a female-biased gene besides sex chromosome and hormonal differences, supplying the current understanding of sex differences in fish.

Two duplicated gsdf homeologs cooperatively regulate male differentiation by inhibiting cyp19a1a transcription in a hexaploid fish.

Although evolutionary fates and expression patterns of duplicated genes have been extensively investigated, how duplicated genes co-regulate a biological process in polyploids remains largely unknown. Here, we identified two gsdf (gonadal somatic cell-derived factor) homeologous genes (gsdf-A and gsdf-B) in hexaploid gibel carp (Carassius gibelio), wherein each homeolog contained three highly conserved alleles. Interestingly, gsdf-A and gsdf-B transcription were mainly activated by dmrt1-A (dsx- and mab-3-related transcription factor 1) and dmrt1-B, respectively. Loss of either gsdf-A or gsdf-B alone resulted in partial male-to-female sex reversal and loss of both caused complete sex reversal, which could be rescued by a nonsteroidal aromatase inhibitor. Compensatory expression of gsdf-A and gsdf-B was observed in gsdf-B and gsdf-A mutants, respectively. Subsequently, we determined that in tissue culture cells, Gsdf-A and Gsdf-B both interacted with Ncoa5 (nuclear receptor coactivator 5) and blocked Ncoa5 interaction with Rora (retinoic acid-related orphan receptor-alpha) to repress Rora/Ncoa5-induced activation of cyp19a1a (cytochrome P450, family 19, subfamily A, polypeptide 1a). These findings illustrate that Gsdf-A and Gsdf-B can regulate male differentiation by inhibiting cyp19a1a transcription in hexaploid gibel carp and also reveal that Gsdf-A and Gsdf-B can interact with Ncoa5 to suppress cyp19a1a transcription in vitro. This study provides a typical case of cooperative mechanism of duplicated genes in polyploids and also sheds light on the conserved evolution of sex differentiation.

Cell Trafficking at the Intersection of the Tumor-Immune Compartments.

Migration is an essential cellular process that regulates human organ development and homeostasis as well as disease initiation and progression. In cancer, immune and tumor cell migration is strongly associated with immune cell infiltration, immune escape, and tumor cell metastasis, which ultimately account for more than 90% of cancer deaths. The biophysics and molecular regulation of the migration of cancer and immune cells have been extensively studied separately. However, accumulating evidence indicates that, in the tumor microenvironment, the motilities of immune and cancer cells are highly interdependent via secreted factors such as cytokines and chemokines. Tumor and immune cells constantly express these soluble factors, which produce a tightly intertwined regulatory network for these cells' respective migration. A mechanistic understanding of the reciprocal regulation of soluble factor-mediated cell migration can provide critical information for the development of new biomarkers of tumor progression and of tumor response to immuno-oncological treatments. We review the biophysical andbiomolecular basis for the migration of immune and tumor cells and their associated reciprocal regulatory network. We also describe ongoing attempts to translate this knowledge into the clinic.

Identification and Characterization of Salt- and Drought-Responsive AQP Family Genes in Medicagosativa L.

Aquaporins (AQP) are distributed ubiquitously in plants, and they play important roles in multiple aspects of plant growth and development, as well as in plant resistance to various environmental stresses. In this study, 43 MsAQP genes were identified in the forage crop Medicago sativa. All the MsAQP proteins were clustered into four subfamilies based on sequence similarity and phylogenetic relationship, including 17 TIPs, 14 NIPs, 9 PIPs and 3 SIPs. Analyses of gene structure and conserved domains indicated that the majority of the deduced MsAQP proteins contained the signature transmembrane domains and the NPA motifs. Analyses on cis-acting elements in the promoter region of MsAQP genes revealed the presence of multiple and diverse stress-responsive and hormone-responsive cis-acting elements. In addition, by analyzing the available and comprehensive gene expression data of M. truncatula, we screened ten representative MtAQP genes that were responsive to NaCl or drought stress. By analyzing the sequence similarity and phylogenetic relationship, we finally identified the corresponding ten salt- or drought-responsive AQP genes in M. sativa, including three MsTIPs, three MsPIPs and four MsNIPs. The qPCRs showed that the relative expression levels of these ten selected MsAQP genes responded differently to NaCl or drought treatment in M. sativa. Gene expression patterns showed that most MsAQP genes were preferentially expressed in roots or in leaves, which may reflect their tissue-specific functions associated with development. Our results lay an important foundation for the future characterization of the functions of MsAQP genes, and provide candidate genes for stress resistance improvement through genetic breeding in M. sativa.

Genomic anatomy of male-specific microchromosomes in a gynogenetic fish.

Unisexual taxa are commonly considered short-lived as the absence of meiotic recombination is supposed to accumulate deleterious mutations and hinder the creation of genetic diversity. However, the gynogenetic gibel carp (Carassius gibelio) with high genetic diversity and wide ecological distribution has outlived its predicted extinction time of a strict unisexual reproduction population. Unlike other unisexual vertebrates, males associated with supernumerary microchromosomes have been observed in gibel carp, which provides a unique system to explore the rationales underlying male occurrence in unisexual lineage and evolution of unisexual reproduction. Here, we identified a massively expanded satellite DNA cluster on microchromosomes of hexaploid gibel carp via comparing with the ancestral tetraploid crucian carp (Carassius auratus). Based on the satellite cluster, we developed a method for single chromosomal fluorescence microdissection and isolated three male-specific microchromosomes in a male metaphase cell. Genomic anatomy revealed that these male-specific microchromosomes contained homologous sequences of autosomes and abundant repetitive elements. Significantly, several potential male-specific genes with transcriptional activity were identified, among which four and five genes displayed male-specific and male-biased expression in gonads, respectively, during the developmental period of sex determination. Therefore, the male-specific microchromosomes resembling common features of sex chromosomes may be the main driving force for male occurrence in gynogenetic gibel carp, which sheds new light on the evolution of unisexual reproduction.

Characterization of the complete chloroplast genome of Alopecurus pratensis L. (Poaceae).

Alopecurus pratensis L. is one of the most important fodder grasses distributed in sub-frigid regions of the world. In this study, the complete chloroplast genome of A. pratensis was deciphered and is 136,157 bp in length. The genome includes a large single-copy region of 80,275 bp, small single-copy region of 12,830 bp, and a pair of inverted repeat regions of 21,526 bp. The GC content of the complete chloroplast genome of A. pratensis is 38.30%. Among the 134 unique genes in the circular genome, 38 tRNA, eight rRNA, and 88 protein-coding genes were annotated. We constructed the Maximum likelihood (ML) tree with 13 species from the Poaceae and found that A. pratensis was phylogenetically related to A. arundinaceus. The published A. pratensis chloroplast genome will provide useful information for phylogenetic and evolutionary study of the genus Alopecurus in the Poaceae.