Genome-wide Analysis Research Articles

Genome wide identification and expression profiling of ATP binding cassette (ABC) transporters gene family in lentil (Lens culinaris Medikus) under aluminium stress condition

Adenosine triphosphate-binding cassette transporters (ABC transporters) are involved in regulating plant growth, development and tolerance to environmental stresses. In this study, a total of 138 ABC transporter genes were identified in the lentil genome that were classified into eight subfamilies. Four lentil ABC transporters from subfamily B and I were clustered together with the previously characterized ABC transporter proteins related to aluminium (Al) detoxification. Lentil ABC transporter genes were distributed across the chromosomes. Tandem duplication was the main driving force for expansion of the ABC gene family. Collinearity of lentil with soybean indicated that ABC gene family is closely linked to Glycine max. ABC genes in the same subfamily showed similar gene structure and conserved motifs. The ABC promoter regions harboured a large number of plant hormones and multiple stress responsive cis-regulatory elements. The qRT-PCR showed that ABC genes had varied expression in roots of lentil at different time points under Al stress. This is the first report on genome wide identification and expression analyses of genes encoding ABC transporter genes in lentil which has provided in-depth insight for future research on evolution and elucidation of molecular mechanisms for aluminium tolerance.

Exploring the relationship between admixture and genetic susceptibility to attention deficit hyperactivity disorder in two Latin American cohorts.

Contemporary research on the genomics of Attention Deficit Hyperactivity Disorder (ADHD) often underrepresents admixed populations of diverse genomic ancestries, such as Latin Americans. This study explores the relationship between admixture and genetic associations for ADHD in Colombian and Mexican cohorts. Some 546 participants in two groups, ADHD and Control, were genotyped with Infinium PsychArray®. Global ancestry levels were estimated using overall admixture proportions and principal component analysis, while local ancestry was determined using a method to estimate ancestral components along the genome. Genome-wide association analysis (GWAS) was conducted to identify significant associations. Differences between Colombia and Mexico were evaluated using appropriate statistical tests. 354 Single-nucleotide polymorphisms (SNPs) and Single-nucleotide variants (SNVs) related to some genes and intergenic regions exhibited suggestive significance (p-value < 5*10e-5) in the GWAS. None of the variants revealed genome-wide significance (p-value < 5*10e-8). The study identified a significant relationship between risk SNPs and the European component of admixture, notably observed in the LOC105379109 gene. Despite differences in risk association loci, such as FOXP2, our findings suggest a possible homogeneity in genetic variation's impact on ADHD between Colombian and Mexican populations. Current reference datasets for ADHD predominantly consist of samples with high European ancestry, underscoring the need for further research to enhance the representation of reference populations and improve the identification of ADHD risk traits in Latin Americans.

Genome-Wide Analysis of the SAUR Gene Family and Its Expression Profiles in Response to Salt Stress in Santalum album.

The SAUR (small auxin-up RNA) family constitutes a category of genes that promptly respond to the hormone auxin and play a pivotal role in diverse biological processes encompassing plant growth and the response to abiotic stress. Santalum album L., a semi-parasitic evergreen tree, is renowned for its economically valuable essential oils, positioning it among the most prized tree species. In this study, a meticulous identification and comprehensive analysis of 43 SAUR genes was conducted within S. album. Based on phylogenetic relationships, the SaSAUR genes were systematically categorized into five groups. A collinearity analysis revealed intriguing insights, disclosing 14 segmental duplications and 9 tandem duplications within the SaSAUR genes, emphasizing the pivotal role of duplication in the expansion of this gene family. Noteworthy variations in the expression levels of SaSAUR genes were observed by delving into the SaSAUR transcriptome data from various tissues, including leaves, roots, and heartwood, as well as under salt-stress conditions. Notably, SaSAUR08 and SaSAUR13 were significantly upregulated in heartwood compared with roots and leaves, while SaSAUR18 was markedly more expressed in roots compared with heartwood and leaves. Furthermore, SaSAUR27 and SaSAUR28 were found to respond closely to salt stress, hinting at their potential involvement in the salt-stress response mechanism. This research offers a comprehensive investigation of SAUR genes in S. album and establishes a foundation for future exploration of the SAUR gene family, particularly its relation to growth and salt-stress responses.

An Accurate Representation of the Number of bZIP Transcription Factors in the Triticum aestivum (Wheat) Genome and the Regulation of Functional Genes during Salt Stress.

Climate change is dramatically increasing the overall area of saline soils around the world, which is increasing by approximately two million hectares each year. Soil salinity decreases crop yields and, thereby, makes farming less profitable, potentially causing increased poverty and hunger in many areas. A solution to this problem is increasing the salt tolerance of crop plants. Transcription factors (TFs) within crop plants represent a key to understanding salt tolerance, as these proteins play important roles in the regulation of functional genes linked to salt stress. The basic leucine zipper (bZIP) TF has a well-documented role in the regulation of salt tolerance. To better understand how bZIP TFs are linked to salt tolerance, we performed a genome-wide analysis in wheat using the Chinese spring wheat genome, which has been assembled by the International Wheat Genome Sequencing Consortium. We identified 89 additional bZIP gene sequences, which brings the total of bZIP gene sequences in wheat to 237. The majority of these 237 sequences included a single bZIP protein domain; however, different combinations of five other domains also exist. The bZIP proteins are divided into ten subfamily groups. Using an in silico analysis, we identified five bZIP genes (ABF2, ABF4, ABI5, EMBP1, and VIP1) that were involved in regulating salt stress. By scrutinizing the binding properties to the 2000 bp upstream region, we identified putative functional genes under the regulation of these TFs. Expression analyses of plant tissue that had been treated with or without 100 mM NaCl revealed variable patterns between the TFs and functional genes. For example, an increased expression of ABF4 was correlated with an increased expression of the corresponding functional genes in both root and shoot tissues, whereas VIP1 downregulation in root tissues strongly decreased the expression of two functional genes. Identifying strategies to sustain the expression of the functional genes described in this study could enhance wheat's salt tolerance.

Spatially revealed roles for lncRNAs in Drosophila spermatogenesis, Y chromosome function and evolution

Unlike coding genes, the number of lncRNA genes in organism genomes is relatively proportional to organism complexity. From plants to humans, the tissues with highest numbers and levels of lncRNA gene expression are the male reproductive organs. To learn why, we initiated a genome-wide analysis of Drosophila lncRNA spatial expression patterns in these tissues. The numbers of genes and levels of expression observed greatly exceed those previously reported, due largely to a preponderance of non-polyadenylated transcripts. In stark contrast to coding genes, the highest numbers of lncRNAs expressed are in post-meiotic spermatids. Correlations between expression levels, localization and previously performed genetic analyses indicate high levels of function and requirement. More focused analyses indicate that lncRNAs play major roles in evolution by controlling transposable element activities, Y chromosome gene expression and sperm construction. A new type of lncRNA-based particle found in seminal fluid may also contribute to reproductive outcomes.

Genome-wide identification and expression pattern analysis of the GRF transcription factor family in Astragalus mongholicus.

Astragalus membranaceus is a plant of the Astragalus genus, which is used as a traditional Chinese herbal medicine with extremely high medicinal and edible value. Astragalus mongholicus, as one of the representative medicinal materials with the same origin of medicine and food, has a rising market demand for its raw materials, but the quality is different in different production areas. Growth-regulating factors (GRF) are transcription factors unique to plants that play important roles in plant growth and development. Up to now, there is no report about GRF in A. mongholicus. This study conducted a genome-wide analysis of the AmGRF gene family, identifying a total of nine AmGRF genes that were classified into subfamily V based on phylogenetic relationships. In the promoter region of the AmGRF gene, we successfully predicted cis-elements that respond to abiotic stress, growth, development, and hormone production in plants. Based on transcriptomic data and real-time quantitative polymerase chain reaction (qPCR) validation, the results showed that AmGRFs were expressed in the roots, stems, and leaves, with overall higher expression in leaves, higher expression of AmGRF1 and AmGRF8 in roots, and high expression levels of AmGRF1 and AmGRF9 in stems. The results of this study provide a theoretical basis for the further exploration of the functions of AmGRFs in plant growth and development.

Genome-wide somatic mutation analysis of sinonasal adenocarcinoma with and without wood dust exposure

BackgroundSinonasal adenocarcinoma is a rare cancer, encompassing two different entities, the intestinal-type sinonasal adenocarcinoma (ITAC) and the non-intestinal-type sinonasal adenocarcinoma (non-ITAC). Occurrence of ITAC is strongly associated with exposure to hardwood dusts. In countries with predominant exposure to softwood dust the occurrence of sinonasal adenocarcinomas is lower and the relative amount of non-ITACs to ITACs is higher. The molecular mechanisms behind the tumorigenic effects of wood dust remain largely unknown.MethodsWe carried out whole-genome sequencing of formalin-fixed paraffin-embedded (FFPE) samples of sinonasal adenocarcinomas from ten wood dust-exposed and six non-exposed individuals, with partial tobacco exposure data. Sequences were analyzed for the presence of mutational signatures matching COSMIC database signatures. Driver mutations and CN variant regions were characterized.ResultsMutation burden was higher in samples of wood dust-exposed patients (p = 0.016). Reactive oxygen species (ROS) damage-related mutational signatures were almost exclusively identified in ITAC subtype samples (p = 0.00055). Tobacco smoke mutational signatures were observed in samples of patients with tobacco exposure or missing information, but not in samples from non-exposed patients. A tetraploidy copy number (CN) signature was enriched in ITAC subtype (p = 0.042). CN variation included recurrent gains in COSMIC Cancer Gene Census genes TERT, SDHA, RAC1, ETV1, PCM1, and MYC. Pathogenic variants were observed most frequently in TP53, NF1, CHD2, BRAF, APC, and LRP1B. Driver mutations and copy number gains did not segregate by subtype.ConclusionsOur analysis identified distinct mutational characteristics in ITAC and non-ITAC. Mutational signature analysis may eventually become useful for documentation of occupation-related cancer, while the exact mechanisms behind wood dust-driven carcinogenesis remain elusive. The presence of homologous recombination deficiency signatures implies a novel opportunity for treatment, but further studies are needed.

Genome-Wide Identification and Expression Profile Analysis of the Phenylalanine Ammonia-Lyase Gene Family in Hevea brasiliensis.

The majority of the world's natural rubber comes from the rubber tree (Hevea brasiliensis). As a key enzyme for synthesizing phenylpropanoid compounds, phenylalanine ammonia-lyase (PAL) has a critical role in plant satisfactory growth and environmental adaptation. To clarify the characteristics of rubber tree PAL family genes, a genome-wide characterization of rubber tree PALs was conducted in this study. Eight PAL genes (HbPAL1-HbPAL8), which spread over chromosomes 3, 7, 8, 10, 12, 13, 14, 16, and 18, were found to be present in the genome of H. brasiliensis. Phylogenetic analysis classified HbPALs into groups I and II, and the group I HbPALs (HbPAL1-HbPAL6) displayed similar conserved motif compositions and gene architectures. Tissue expression patterns of HbPALs quantified by quantitative real-time PCR (qPCR) proved that distinct HbPALs exhibited varying tissue expression patterns. The HbPAL promoters contained a plethora of cis-acting elements that responded to hormones and stress, and the qPCR analysis demonstrated that abiotic stressors like cold, drought, salt, and H2O2-induced oxidative stress, as well as hormones like salicylic acid, abscisic acid, ethylene, and methyl jasmonate, controlled the expression of HbPALs. The majority of HbPALs were also regulated by powdery mildew, anthracnose, and Corynespora leaf fall disease infection. In addition, HbPAL1, HbPAL4, and HbPAL7 were significantly up-regulated in the bark of tapping panel dryness rubber trees relative to that of healthy trees. Our results provide a thorough comprehension of the characteristics of HbPAL genes and set the groundwork for further investigation of the biological functions of HbPALs in rubber trees.

Genome-wide identification and functional analysis of the SiCIN gene family in foxtail millet (Setaria italica L.)

Cell wall invertase (CIN) is a vital member of plant invertase (INV) and plays a key role in the breakdown of sucrose. This enzyme facilitates the hydrolysis of sucrose into glucose and fructose, which is crucial for various aspects of plant growth and development. However, the function of CIN genes in foxtail millet (Setaria italica) is less studied. In this research, we used the blast-p of NCBI and TBtools for bidirectional comparison, and a total of 13 CIN genes (named SiCINs) were identified from foxtail millet by using Arabidopsis and rice CIN sequences as reference sequences. The phylogenetic tree analysis revealed that the CIN genes can be categorized into three subfamilies: group 1, group 2, and group 3. Furthermore, upon conducting chromosomal localization analysis, it was observed that the 13 SiCINs were distributed unevenly across five chromosomes. Cis-acting elements of SiCIN genes can be classified into three categories: plant growth and development, stress response, and hormone response. The largest number of cis-acting elements were those related to light response (G-box) and the cis-acting elements related to seed-specific regulation (RY-element). qRT-PCR analysis further confirmed that the expression of SiCIN7 and SiCIN8 in the grain was higher than that in any other tissues. The overexpression of SiCIN7 in Arabidopsis improved the grain size and thousand-grain weight, suggesting that SiCIN7 could positively regulate grain development. Our findings will help to further understand the grain-filling mechanism of SiCIN and elucidate the biological mechanism underlying the grain development of SiCIN.

Genome-Wide Analysis of Transcription Factor R2R3-MYB Gene Family and Gene Expression Profiles during Anthocyanin Synthesis in Common Walnut (Juglans regia L.).

The R2R3-MYB gene family, encoding plant transcriptional regulators, participates in many metabolic pathways of plant physiology and development, including flavonoid metabolism and anthocyanin synthesis. This study proceeded as follows: the JrR2R3-MYB gene family was analyzed genome-wide, and the family members were identified and characterized using the high-quality walnut reference genome "Chandler 2.0". All 204 JrR2R3-MYBs were established and categorized into 30 subgroups via phylogenetic analysis. JrR2R3-MYBs were unevenly distributed over 16 chromosomes. Most JrR2R3-MYBs had similar structures and conservative motifs. The cis-acting elements exhibit multiple functions of JrR2R3-MYBs such as light response, metabolite response, and stress response. We found that the expansion of JrR2R3-MYBs was mainly caused by WGD or segmental duplication events. Ka/Ks analysis indicated that these genes were in a state of negative purifying selection. Transcriptome results suggested that JrR2R3-MYBs were widely entangled in the process of walnut organ development and differentially expressed in different colored varieties of walnuts. Subsequently, we identified 17 differentially expressed JrR2R3-MYBs, 9 of which may regulate anthocyanin biosynthesis based on the results of a phylogenetic analysis. These genes were present in greater expression levels in 'Zijing' leaves than in 'Lvling' leaves, as revealed by the results of qRT-PCR experiments. These results contributed to the elucidation of the functions of JrR2R3-MYBs in walnut coloration. Collectively, this work provides a foundation for exploring the functional characteristics of the JrR2R3-MYBs in walnuts and improving the nutritional value and appearance quality of walnuts.

PLA2G12A protects against diet-induced obesity and insulin resistance by enhancing energy expenditure and clearance of circulating triglycerides.

Secreted phospholipase A2s are involved in the development of obesity, type 2 diabetes mellitus (T2DM) and cardiovascular disease, which have become serious and growing health concerns worldwide. Integration of genome-wide association study and gene co-expression networks analysis showed that the secreted phospholipase A2 group XIIA (PLA2G12A) may participate in hepatic lipids metabolism. Nevertheless, the role of PLA2G12A in lipid metabolism and its potential mechanism remain elusive. Here, we used AAV9 vector carrying human PLA2G12A gene to exogenously express hPLA2G12A in the liver of mice. We demonstrated that the overexpression of hPLA2G12A resulted in a significant decrease in serum lipid levels in wild-type mice fed with chow diet or high-fat diet (HFD). Moreover, hPLA2G12A treatment protected against diet-induced obesity and insulin resistance in mice fed a HFD. Notably, we found that hPLA2G12A treatment confers protection against obesity and hyperlipidemia independent of its enzymatic activity, but rather by increasing physical activity and energy expenditure. Furthermore, we demonstrated that hPLA2G12A treatment induced upregulation of ApoC2 and Cd36 and downregulation of Angptl8, which contributed to the increase in clearance of circulating triglycerides and hepatic uptake of fatty acids without affecting hepatic de novo lipogenesis, very low-density lipoprotein secretion, or intestinal lipid absorption. Our study highlights the potential of PLA2G12A gene therapy as a promising approach for treating obesity, insulin resistance and T2DM.

Genome wide association analysis reveals the genetic sites of heat resistance in largemouth bass (Micropterus salmoides)

Largemouth bass (Micropterus salmoides) is an important commercially cultured freshwater fish. High temperature is a common challenge encountered in the summer cultivation of largemouth bass. In this study, we conducted a genome-wide association study using general linear models and linear mixed models to identify significant single nucleotide polymorphisms (SNPs) and candidate genes associated with loss of equilibrium (LOE) and weight gain (WG) during high temperature challenge. We collected fin tissue from 353 individuals of largemouth bass exposed to high water temperature (37.5 °C) and performed whole-genome resequencing. Through association analysis, we successfully identified 10 SNPs and 38 candidate genes that were significantly associated with LOE, and four SNPs and 19 candidate genes that were significantly associated with WG. Using quantitative PCR and Western blot, we tested the expression of eight genes related to LOE in heat-sensitive and heat-tolerant individuals, and identified significant differences in the expression of the following five genes: Tmub2, Ngfr, Hgf, Hdac6, and Sypl. For WG, we assessed the expression of three candidate genes, among which Mrc2l displayed significant differences in expression between the high WG group and low WG group. We discuss the functions and potential biological processes of candidate genes for LOE and WG. Our results provide basic genetic information for breeding largemouth bass with high temperature tolerance.

Variation in WIDTH OF LEAF AND GRAIN contributes to grain and leaf size by controlling LARGE2 stability in rice.

Grain and flag leaf size are two important agronomic traits that influence grain yield in rice (Oryza sativa). Many QTLs and genes that regulate these traits individually have been identified, however, few QTLs and genes that simultaneously control these two traits have been identified. In this study, we conducted a genome-wide association analysis in rice and detected a major locus, WIDTH OF LEAF AND GRAIN (WLG), that associated with both grain and flag leaf width. WLG encodes a RING-domain E3 ubiquitin ligase. WLGhap.B, which possesses five SNP variations compared to WLGhap.A, encodes a protein with enhanced ubiquitination activity that confers increased rice leaf width and grain size, whereas mutation of WLG leads to narrower leaves and smaller grains. Both WLGhap.A and WLGhap.B interact with LARGE2, a HETC-type E3 ligase, however, WLGhap.B exhibits stronger interaction with LARGE2, thus higher ubiquitination activity towards LARGE2 compared with WLGhap.A. Lysine1021 is crucial for the ubiquitination of LARGE2 by WLG. Loss-of-function of LARGE2 in wlg-1 phenocopies large2-c in grain and leaf width, suggesting that WLG acts upstream of LARGE2. These findings reveal the genetic and molecular mechanism by which the WLG-LARGE2 module mediates grain and leaf size in rice, and suggest the potential of WLGhap.B in improving rice yield.

A comprehensive genome-wide cross-trait analysis of sexual factors and uterine leiomyoma.

Age at first sexual intercourse (AFS) and lifetime number of sexual partners (NSP) may influence the pathogenesis of uterine leiomyoma (UL) through their associations with hormonal concentrations and uterine infections. Leveraging summary statistics from large-scale genome-wide association studies conducted in European ancestry for each trait (NAFS = 214,547; NNSP = 370,711; NUL = 302,979), we observed a significant negative genomic correlation for UL with AFS (rg = -0.11, P = 7.83×10-4), but not with NSP (rg = 0.01, P = 0.62). Four specific genomic regions were identified as contributing significant local genetic correlations to AFS and UL, including one genomic region further identified for NSP and UL. Partitioning SNP-heritability with cell-type-specific annotations, a close clustering of UL with both AFS and NSP was identified in immune and blood-related components. Cross-trait meta-analysis revealed 15 loci shared between AFS/NSP and UL, including 7 novel SNPs. Univariable two-sample Mendelian randomization (MR) analysis suggested no evidence for a causal association between genetically predicted AFS/NSP and risk of UL, nor vice versa. Multivariable MR adjusting for age at menarche or/and age at natural menopause revealed a significant causal effect of genetically predicted higher AFS on a lower risk of UL. Such effect attenuated to null when age at first birth was further included. Utilizing participant-level data from the UK Biobank, one-sample MR based on genetic risk scores yielded consistent null findings among both pre-menopausal and post-menopausal females. From a genetic perspective, our study demonstrates an intrinsic link underlying sexual factors (AFS and NSP) and UL, highlighting shared biological mechanisms rather than direct causal effects. Future studies are needed to elucidate the specific mechanisms involved in the shared genetic influences and their potential impact on UL development.

Genome-wide identification and analysis of UBP gene family in Quercus robur

Genome-wide identification and analysis of UBP gene family in Quercus robur

Genomic Testing in Patients with Kidney Failure of an Unknown Cause: a National Australian Study

Background The cause of kidney failure is unknown in approximately 10% of patients with stage 5 chronic kidney disease (CKD). For those who first present to nephrology care with kidney failure, standard investigations of serology, imaging, urinalysis and kidney biopsy are limited differentiators of etiology. We aimed to determine the diagnostic utility of whole-genome sequencing (WGS) with analysis of a broad kidney gene panel in patients with kidney failure of unknown cause. Methods We prospectively recruited 100 participants who reached CKD stage 5 at 50 years of age and had an unknown cause of kidney failure after standard investigation. Clinically-accredited WGS was performed in this national cohort after genetic counselling. The primary analysis was targeted to 388 kidney-related genes with second-tier genome-wide and mitochondrial analysis. Results The cohort was 61% male and the average age of participants at stage 5 CKD was 32 years (9 months to 50 years). A genetic diagnosis was made in 25% of participants. Disease-causing variants were identified across autosomal dominant tubulointerstitial kidney disease (6), glomerular disorders (4), ciliopathies (3), tubular disorders (2), Alport syndrome (4) and mitochondrial disease (1). Most diagnoses (80%) were in autosomal dominant, X-linked or mitochondrial conditions (UMOD; COL4A5; INF2; CLCN5; TRPC6; COL4A4; EYA1; HNF1B; WT1; NBEA; m.3243A&gt;G). Patients with a family history of CKD were more likely to have a positive result (OR 3.29, 95% CI 1.10-11.29). Thirteen percent of participants without a CKD family history had a positive result. In those who first presented in stage 5 CKD, WGS with broad analysis of a curated kidney-disease gene panel was diagnostically more informative than kidney biopsy, with biopsy being inconclusive in 24 of 25 participants. Conclusions In this prospectively ascertained Australian cohort, we identified a genetic diagnosis in 25% of patients with kidney failure of unknown cause.

Epigenetically active chromatin in neonatal iWAT reveals GABPα as a potential regulator of beige adipogenesis.

Thermogenic beige adipocytes, which dissipate energy as heat, are found in neonates and adults. Recent studies show that neonatal beige adipocytes are highly plastic and contribute to >50% of beige adipocytes in adults. Neonatal beige adipocytes are distinct from recruited beige adipocytes in that they develop independently of temperature and sympathetic innervation through poorly defined mechanisms. We characterized the neonatal beige adipocytes in the inguinal white adipose tissue (iWAT) of C57BL6 postnatal day 3 and 20 mice (P3 and P20) by imaging, genome-wide RNA-seq analysis, ChIP-seq analysis, qRT-PCR validation, and biochemical assays. We found an increase in acetylated histone 3 lysine 27 (H3K27ac) on the promoter and enhancer regions of beige-specific gene UCP1 in iWAT of P20 mice. Furthermore, H3K27ac ChIP-seq analysis in the iWAT of P3 and P20 mice revealed strong H3K27ac signals at beige adipocyte-associated genes in the iWAT of P20 mice. The integration of H3K27ac ChIP-seq and RNA-seq analysis in the iWAT of P20 mice reveal epigenetically active signatures of beige adipocytes, including oxidative phosphorylation and mitochondrial metabolism. We identify the enrichment of GA-binding protein alpha (GABPα) binding regions in the epigenetically active chromatin regions of the P20 iWAT, particularly on beige genes, and demonstrate that GABPα is required for beige adipocyte differentiation. Moreover, transcriptomic analysis and glucose oxidation assays revealed increased glycolytic activity in the neonatal iWAT from P20. Our findings demonstrate that epigenetic mechanisms regulate the development of peri-weaning beige adipocytes via GABPα. Further studies to better understand the upstream mechanisms that regulate epigenetic activation of GABPα and characterization of the metabolic identity of neonatal beige adipocytes will help us harness their therapeutic potential in metabolic diseases.

Genome-wide identification and expression analysis of the SET domain-containing gene family in potato (Solanum tuberosum L.)

Genes containing the SET domain can catalyse histone lysine methylation, which in turn has the potential to cause changes to chromatin structure and regulation of the transcription of genes involved in diverse physiological and developmental processes. However, the functions of SET domain-containing (StSET) genes in potato still need to be studied. The objectives of our study can be summarized as in silico analysis to (i) identify StSET genes in the potato genome, (ii) systematically analyse gene structure, chromosomal distribution, gene duplication events, promoter sequences, and protein domains, (iii) perform phylogenetic analyses, (iv) compare the SET domain-containing genes of potato with other plant species with respect to protein domains and orthologous relationships, (v) analyse tissue-specific expression, and (vi) study the expression of StSET genes in response to drought and heat stresses. In this study, we identified 57 StSET genes in the potato genome, and the genes were physically mapped onto eleven chromosomes. The phylogenetic analysis grouped these StSET genes into six clades. We found that tandem duplication through sub-functionalisation has contributed only marginally to the expansion of the StSET gene family. The protein domain TDBD (PFAM ID: PF16135) was detected in StSET genes of potato while it was absent in all other previously studied species. This study described three pollen-specific StSET genes in the potato genome. Expression analysis of four StSET genes under heat and drought in three potato clones revealed that these genes might have non-overlapping roles under different abiotic stress conditions and durations. The present study provides a comprehensive analysis of StSET genes in potatoes, and it serves as a basis for further functional characterisation of StSET genes towards understanding their underpinning biological mechanisms in conferring stress tolerance.

Genome-wide prediction and functional analysis of WOX genes in blueberry.

WOX genes are a class of plant-specific transcription factors. The WUSCHEL-related homeobox (WOX) family is a member of the homeobox transcription factor superfamily. Previous studies have shown that WOX members play important roles in plant growth and development. However, studies of the WOX gene family in blueberry plants have not been reported. In order to understand the biological function of the WOX gene family in blueberries, bioinformatics were used methods to identify WOX gene family members in the blueberry genome, and analyzed the basic physical and chemical properties, gene structure, gene motifs, promoter cis-acting elements, chromosome location, evolutionary relationships, expression pattern of these family members and predicted their functions. Finally, 12 genes containing the WOX domain were identified and found to be distributed on eight chromosomes. Phylogenetic tree analysis showed that the blueberry WOX gene family had three major branches: ancient branch, middle branch, and WUS branch. Blueberry WOX gene family protein sequences differ in amino acid number, molecular weight, isoelectric point and hydrophobicity. Predictive analysis of promoter cis-acting elements showed that the promoters of the VdWOX genes contained abundant light response, hormone, and stress response elements. The VdWOX genes were induced to express in both stems and leaves in response to salt and drought stress. Our results provided comprehensive characteristics of the WOX gene family and important clues for further exploration of its role in the growth, development and resistance to various stress in blueberry plants.

Genome Evolution and Introgression in the New Zealand mud Snails Potamopyrgus estuarinus and Potamopyrgus kaitunuparaoa.

We have sequenced, assembled, and analyzed the nuclear and mitochondrial genomes and transcriptomes of Potamopyrgus estuarinus and Potamopyrgus kaitunuparaoa, two prosobranch snail species native to New Zealand that together span the continuum from estuary to freshwater. These two species are the closest known relatives of the freshwater species Potamopyrgus antipodarum-a model for studying the evolution of sex, host-parasite coevolution, and biological invasiveness-and thus provide key evolutionary context for understanding its unusual biology. The P. estuarinus and P. kaitunuparaoa genomes are very similar in size and overall gene content. Comparative analyses of genome content indicate that these two species harbor a near-identical set of genes involved in meiosis and sperm functions, including seven genes with meiosis-specific functions. These results are consistent with obligate sexual reproduction in these two species and provide a framework for future analyses of P. antipodarum-a species comprising both obligately sexual and obligately asexual lineages, each separately derived from a sexual ancestor. Genome-wide multigene phylogenetic analyses indicate that P. kaitunuparaoa is likely the closest relative to P. antipodarum. We nevertheless show that there has been considerable introgression between P. estuarinus and P. kaitunuparaoa. That introgression does not extend to the mitochondrial genome, which appears to serve as a barrier to hybridization between P. estuarinus and P. kaitunuparaoa. Nuclear-encoded genes whose products function in joint mitochondrial-nuclear enzyme complexes exhibit similar patterns of nonintrogression, indicating that incompatibilities between the mitochondrial and the nuclear genome may have prevented more extensive gene flow between these two species.