Gene Expression Domains Research Articles

Molecular mechanisms and behavioral relevance underlying neural correlates of childhood neglect

BackgroundChildhood neglect is associated with brain changes, yet the molecular mechanisms and behavioral relevance underlying such associations remain elusive. MethodsWe calculated fractional amplitude of low-frequency fluctuations (fALFF) using resting-state functional MRI and tested their correlation with childhood neglect across a large sample of 510 healthy young adults. Then, we investigated the spatial relationships of the identified neural correlates of childhood neglect with gene expression, neurotransmitter, and behavioral domain atlases. ResultsWe found that more severe childhood neglect was correlated with higher fALFF in the bilateral anterior cingulate cortex. Remarkably, the identified neural correlates of childhood neglect were spatially correlated with expression of gene categories primarily involving neuron, synapse, ion channel, cognitive and perceptual processes, and physiological response and regulation. Concurrently, there were significant associations between the neural correlates and specific neurotransmitter systems including serotonin and GABA. Finally, neural correlates of childhood neglect were associated with diverse behavioral domains implicating mental disorders, emotion, cognition, and sensory perception. LimitationsThe cross-sectional study design cannot unequivocally establish causality. ConclusionsOur findings may not only add to the current knowledge regarding the relationship between childhood neglect and mental health, but also have clinical implications for developing preventive strategies for individuals exposed to childhood neglect who are at risk for mental disorders.

Unfolding the ventral nerve center of chaetognaths

BackgroundChaetognaths are a clade of marine worm-like invertebrates with a heavily debated phylogenetic position. Their nervous system superficially resembles the protostome type, however, knowledge regarding the molecular processes involved in neurogenesis is lacking. To better understand these processes, we examined the expression profiles of marker genes involved in bilaterian neurogenesis during post-embryonic stages of Spadella cephaloptera. We also investigated whether the transcription factor encoding genes involved in neural patterning are regionally expressed in a staggered fashion along the mediolateral axis of the nerve cord as it has been previously demonstrated in selected vertebrate, insect, and annelid models.MethodsThe expression patterns of genes involved in neural differentiation (elav), neural patterning (foxA, nkx2.2, pax6, pax3/7, and msx), and neuronal function (ChAT and VAChT) were examined in S. cephaloptera hatchlings and early juveniles using whole-mount fluorescent in situ hybridization and confocal microscopy.ResultsThe Sce-elav+ profile of S. cephaloptera hatchlings reveals that, within 24 h of post-embryonic development, the developing neural territories are not limited to the regions previously ascribed to the cerebral ganglion, the ventral nerve center (VNC), and the sensory organs, but also extend to previously unreported CNS domains that likely contribute to the ventral cephalic ganglia. In general, the neural patterning genes are expressed in distinct neural subpopulations of the cerebral ganglion and the VNC in hatchlings, eventually becoming broadly expressed with reduced intensity throughout the CNS in early juveniles. Neural patterning gene expression domains are also present outside the CNS, including the digestive tract and sensory organs. ChAT and VAChT domains within the CNS are predominantly observed in specific subpopulations of the VNC territory adjacent to the ventral longitudinal muscles in hatchlings.ConclusionsThe observed spatial expression domains of bilaterian neural marker gene homologs in S. cephaloptera suggest evolutionarily conserved roles in neurogenesis for these genes among bilaterians. Patterning genes expressed in distinct regions of the VNC do not show a staggered medial-to-lateral expression profile directly superimposable to other bilaterian models. Only when the VNC is conceptually laterally unfolded from the longitudinal muscle into a flat structure, an expression pattern bearing resemblance to the proposed conserved bilaterian mediolateral regionalization becomes noticeable. This finding supports the idea of an ancestral mediolateral patterning of the trunk nervous system in bilaterians.

Neural Correlates of Early-Life Urbanization and Their Spatial Relationships with Gene Expression, Neurotransmitter, and Behavioral Domain Atlases.

Previous neuroimaging research has established associations between urban exposure during early life and alterations in brain function and structure. However, the molecular mechanisms and behavioral relevance of these associations remain largely unknown. Here, we aimed to address this question using a combined analysis of multimodal data. Initially, we calculated amplitude of low-frequency fluctuations (ALFF) and gray matter volume (GMV) using resting-state functional and structural MRI to investigate their associations with early-life urbanization in a large sample of 511 healthy young adults. Then, we examined the spatial relationships of the identified neural correlates of early-life urbanization with gene expression, neurotransmitter, and behavioral domain atlases. Results showed that higher early-life urbanization scores were correlated with increased ALFF of the right fusiform gyrus and decreased GMV of the left dorsal medial prefrontal cortex and left precuneus. Remarkably, the identified neural correlates of early-life urbanization were spatially correlated with expression of gene categories primarily involving immune system process, signal transduction, and cellular metabolic process. Concurrently, there were significant associations between the neural correlates and specific neurotransmitter systems including dopamine, acetylcholine, and serotonin. Finally, we found that the ALFF correlates were associated with behavioral terms including "perception," "sensory," "cognitive control," and "reasoning." Apart from expanding existing knowledge of early-life urban environmental risk for mental disorders and health in general, our findings may contribute to an emerging framework for integrating social science, neuroscience, genetics, and public policy to respond to the major health challenge of world urbanization.

Evolution of two-pore domain potassium channels and their gene expression in zebrafish embryos.

The two-pore domain potassium (K2P) channels are a major type of potassium channels that maintain the cell membrane potential by conducting passive potassium leak currents independent of voltage change. They play prominent roles in multiple physiological processes, including neuromodulation, perception of pain, breathing and mood control, and response to volatile anesthetics. Mutations in K2P channels have been linked to many human diseases, such as neuronal and cardiovascular disorders and cancers. Significant progress has been made to understand their protein structures, physiological functions, and pharmacological modifiers. However, their expression and function during embryonic development remain largely unknown. We employed the zebrafish model and identified 23 k2p genes using BLAST search and gene cloning. We first analyzed vertebrate K2P channel evolution by phylogenetic and syntenic analyses. Our data revealed that the six subtypes of the K2P genes have already evolved in invertebrates long before the emergence of vertebrates. Moreover, the vertebrate K2P gene number increased, most likely due to two whole-genome duplications. Furthermore, we examined zebrafish k2p gene expression during early embryogenesis by in situ hybridization. Each subgroup's genes showed similar but distinct gene expression domains with some exceptions. Most of them were expressed in neural tissues consistent with their known function of neural excitability regulation. However, a few k2p genes were expressed temporarily in specific tissues or organs, suggesting that these K2P channels may be needed for embryonic development. Our phylogenetic and developmental analyses of K2P channels shed light on their evolutionary history and potential roles during embryogenesis related to their physiological functions and human channelopathies.

A Rhodnius prolixus catalytically inactive Calpain protease patterns the insect embryonic dorsal-ventral axis

The calcium dependent Calpain proteases are modulatory enzymes with important roles in cell cycle control, development and immunity. In the fly model Drosophila melanogaster Calpain A cleaves Cactus/IkappaB and consequently modifies Toll signals during embryonic dorsal-ventral (DV) patterning. Here we explore the role of Calpains in the hemiptera Rhodnius prolixus, an intermediate germband insect where the Bone Morphogenetic Protein (BMP) instead of the Toll pathway plays a major role in DV patterning. Phylogenetic analysis of Calpains in species ranging from Isoptera to Diptera indicates an increase of Calpain sequences in the R. prolixus genome and other hemimetabolous species. One locus encoding each of the CalpC, CalpD and Calp7 families, and seven Calpain A/B loci are present in the R. prolixus genome. Several predicted R. prolixus Calpains display a unique architecture, such as loss of Calcium-binding EF-hand domains and loss of catalytic residues in the active site CysPc domain, yielding catalytically dead Calpains A/B. Knockdown for one of these inactive Calpains results in embryonic DV patterning defects, with expansion of ventral and lateral gene expression domains and consequent failure of germ band elongation. In conclusion, our results reveal that Calpains may exert a conserved function in insect DV patterning, despite the changing role of the Toll and BMP pathways in defining gene expression territories along the insect DV axis.

Ultrastructure of the lamprey head mesoderm reveals evolution of the vertebrate head

The cranial muscle is a critical component in the vertebrate head for a predatory lifestyle. However, its evolutionary origin and possible segmental nature during embryogenesis have been controversial. In jawed vertebrates, the presence of pre-otic segments similar to trunk somites has been claimed based on developmental observations. However, evaluating such arguments has been hampered by the paucity of research on jawless vertebrates. Here, we discovered different cellular arrangements in the head mesoderm in lamprey embryos (Lethenteron camtschaticum) using serial block-face scanning electron and laser scanning microscopies. These cell populations were morphologically and molecularly different from somites. Furthermore, genetic comparison among deuterostomes revealed that mesodermal gene expression domains were segregated antero-posteriorly in vertebrates, whereas such segregation was not recognized in invertebrate deuterostome embryos. These findings indicate that the vertebrate head mesoderm evolved from the anteroposterior repatterning of an ancient mesoderm and developmentally diversified before the split of jawless and jawed vertebrates.

The sulfotransferase XB5850668.L is required to apportion embryonic ectodermal domains.

Members of the sulfotransferase superfamily (SULT) influence the activity of a wide range of hormones, neurotransmitters, metabolites and xenobiotics. However, their roles in developmental processes are not well characterized even though they are expressed during embryogenesis. We previously found in a microarray screen that Six1 up-regulates LOC100037047, which encodes XB5850668.L, an uncharacterized sulfotransferase. Since Six1 is required for patterning the embryonic ectoderm into its neural plate, neural crest, preplacodal and epidermal domains, we used loss- and gain-of function assays to characterize the role of XB5850668.L during this process. Knockdown of endogenous XB5850668.L resulted in the reduction of epidermal, neural crest, cranial placode and otic vesicle gene expression domains, concomitant with neural plate expansion. Increased levels had minimal effects, but infrequently expanded neural plate and neural crest gene domains, and infrequently reduced cranial placode and otic vesicle gene domains. Mutation of two key amino acids in the sulfotransferase catalytic domain required for PAPS binding and enzymatic activity tended to reduce the effects of overexpressing the wild-type protein. Our analyses indicates that XB5850668.L is a member of the SULT2 family that plays important roles in patterning the embryonic ectoderm. Some aspects of its influence likely depend on sulfotransferase activity.

How enhancers regulate wavelike gene expression patterns.

A key problem in development is to understand how genes turn on or off at the right place and right time during embryogenesis. Such decisions are made by non-coding sequences called 'enhancers.' Much of our models of how enhancers work rely on the assumption that genes are activated de novo as stable domains across embryonic tissues. Such a view has been strengthened by the intensive landmark studies of the early patterning of the anterior-posterior (AP) axis of the Drosophila embryo, where indeed gene expression domains seem to arise more or less stably. However, careful analysis of gene expression patterns in other model systems (including the AP patterning in vertebrates and short-germ insects like the beetle Tribolium castaneum) painted a different, very dynamic view of gene regulation, where genes are oftentimes expressed in a wavelike fashion. How such gene expression waves are mediated at the enhancer level is so far unclear. Here, we establish the AP patterning of the short-germ beetle Tribolium as a model system to study dynamic and temporal pattern formation at the enhancer level. To that end, we established an enhancer prediction system in Tribolium based on time- and tissue-specific ATAC-seq and an enhancer live reporter system based on MS2 tagging. Using this experimental framework, we discovered several Tribolium enhancers, and assessed the spatiotemporal activities of some of them in live embryos. We found our data consistent with a model in which the timing of gene expression during embryonic pattern formation is mediated by a balancing act between enhancers that induce rapid changes in gene expression patterns (that we call 'dynamic enhancers') and enhancers that stabilize gene expression patterns (that we call 'static enhancers'). However, more data is needed for a strong support for this or any other alternative models.

Levels of APE1 Repair Gene and CLEC4M in Lung Cancer Patients Receiving Cisplatin Chemotherapy.

This study aimed to detect the levels of apurinic/apyrimidinic endonuclease 1 (APE1) gene expression and C-type lectin domain family 4 member M (CLEC4M) and their association with cisplatin chemotherapy in lung cancer patients. Overall, 105 individuals who attended the Al-Amal National Hospital for Cancer Management, Baghdad, Iraq, were enrolled in the study and divided into three equal groups. The groups included the patients newly diagnosed with lung cancer, cancer patients who received cisplatin, and the healthy control group. All study groups were subjected to the sampling of the venous blood for molecular analysis by real-time polymerase chain reaction (RT-PCR) to detect the APE1 gene and enzyme-linked immunosorbent assay (ELISA) for serological testing to measure the concentration of CLEC4M protein. Significantly, the values of both cancer groups were higher than those reported in the control group. The relative index revealed a significant difference in the mean fold change level of APE1 in the newly diagnosed group (3 fold) and cisplatin therapy patients group (2 fold), compared to the control group (P=0.005). No significant differences were detected between the two cancer groups in terms of fold change mean of expression, demographic characteristics, and cancer histological type. Regarding human CLEC4M protein level, cases receiving cisplatin (139.2±25.9) and newly diagnosed patients (331.0±38.1) had a highly significant difference with the control group (100.3±47.5, P<0.001). There was no significant difference between the concentration level of CLEC4M and all parameters in demographic characteristics and cancer histological type. This was the first study to demonstrate that higher expression levels of new APE1, CLEC4M, and glutathione, especially after chemotherapy, are beneficial as diagnostic and prognostic markers for resistance to platinum chemotherapy in Iraqi lung cancer patients.

Functional regionalization of the differentiating cerebellar Purkinje cell population occurs in an activity-dependent manner.

The cerebellum is organized into functional regions each dedicated to process different motor or sensory inputs for controlling different locomotor behaviors. This functional regionalization is prominent in the evolutionary conserved single-cell layered Purkinje cell (PC) population. Fragmented gene expression domains suggest a genetic organization of PC layer regionalization during cerebellum development. However, the establishment of such functionally specific domains during PC differentiation remained elusive. We show the progressive emergence of functional regionalization of PCs from broad responses to spatially restricted regions in zebrafish by means of in vivo Ca2+-imaging during stereotypic locomotive behavior. Moreover, we reveal that formation of new dendritic spines during cerebellar development using in vivo imaging parallels the time course of functional domain development. Pharmacological as well as cell-type specific optogenetic inhibition of PC neuronal activity results in reduced PC dendritic spine density and an altered stagnant pattern of functional domain formation in the PC layer. Hence, our study suggests that functional regionalization of the PC layer is driven by physiological activity of maturing PCs themselves.

Further refining the boundaries of the hippocampus CA2 with gene expression and connectivity: Potential subregions and heterogeneous cell types.

Over the last two decades, the definition of hippocampal area CA2 has evolved from Lorente de Nó's original Golgi-based morphological description with the discovery of specific CA2 gene expression markers. Exploiting the specificity of these molecules has allowed for the genetic dissection of CA2 structure and function in transgenic mice. With this change in criteria, the anatomical boundaries of the CA2 have expanded across the hippocampal axis but the CA2's full rostrocaudal extent is not consistently delineated across atlases. The Hippocampus Gene Expression Atlas (HGEA) provides a comprehensive map of 20 gene expression domains across the entire mouse hippocampus including the CA2. In this commentary, I will review the consensus gene expression patterns that demarcate the expanded CA2 boundaries in the HGEA. Using DropViz single-cell transcriptomics and Mouse Connectome Project connectomics data, I will then suggest potential differences in CA2 cell type heterogeneity and connectivity that may identify and characterize further CA2 subregions.

Differential carbon utilization enables co-existence of recently speciated Campylobacteraceae in the cow rumen epithelial microbiome

The activities of different microbes in the cow rumen have been shown to modulate the host’s ability to utilize plant biomass, while the host–rumen interface has received little attention. As datasets collected worldwide have pointed to Campylobacteraceae as particularly abundant members of the rumen epithelial microbiome, we targeted this group in a subset of seven cows with meta- and isolate genome analysis. We show that the dominant Campylobacteraceae lineage has recently speciated into two populations that were structured by genome-wide selective sweeps followed by population-specific gene import and recombination. These processes led to differences in gene expression and enzyme domain composition that correspond to the ability to utilize acetate, the main carbon source for the host, at the cost of inhibition by propionate. This trade-off in competitive ability further manifests itself in differential dynamics of the two populations in vivo. By exploring population-level adaptations that otherwise remain cryptic in culture-independent analyses, our results highlight how recent evolutionary dynamics can shape key functional roles in the rumen microbiome.

Quantifying Gene Expression Domains in Plant Shoot Apical Meristems.

The shoot apical meristem is the plant tissue that produces the plant aerial organs such as flowers and leaves. To better understand how the shoot apical meristem develops and adapts to the environment, imaging developing shoot meristems expressing fluorescence reporters through laser confocal microscopy is becoming increasingly important. Yet, there are not many computational pipelines enabling a systematic and high-throughput characterization of the produced microscopy images. This chapter provides a simple method to analyze 3D images obtained through laser scanning microscopy and quantitatively characterize radially or axially symmetric 3D fluorescence domains expressed in a tissue or organ by a reporter. Then, it presents different computational pipelines aiming at performing high-throughput quantitative image analysis of gene expression in plant inflorescence and floral meristems. This methodology has notably enabled the quantitative characterization of how stem cells respond to environmental perturbations in the Arabidopsisthaliana inflorescence meristem and will open new avenues in the use of quantitative analysis of gene expression in shoot apical meristems. Overall, the presented methodology provides a simple framework to analyze quantitatively gene expression domains from 3D confocal images at the tissue and organ level, which can be applied to shoot meristems and other organs and tissues.

Seeking Sense in the Hox Gene Cluster.

The Hox gene cluster, responsible for patterning of the head-tail axis, is an ancestral feature of all bilaterally symmetrical animals (the Bilateria) that remains intact in a wide range of species. We can say that the Hox cluster evolved successfully only once since it is commonly the same in all groups, with labial-like genes at one end of the cluster expressed in the anterior embryo, and Abd-B-like genes at the other end of the cluster expressed posteriorly. This review attempts to make sense of the Hox gene cluster and to address the following questions. How did the Hox cluster form in the protostome-deuterostome last common ancestor, and why was this with a particular head-tail polarity? Why is gene clustering usually maintained? Why is there collinearity between the order of genes along the cluster and the positions of their expressions along the embryo? Why do the Hox gene expression domains overlap along the embryo? Why have vertebrates duplicated the Hox cluster? Why do Hox gene knockouts typically result in anterior homeotic transformations? How do animals adapt their Hox clusters to evolve new structural patterns along the head-tail axis?

The SmMYB36-SmERF6/SmERF115 module regulates the biosynthesis of tanshinones and phenolic acids in salvia miltiorrhiza hairy roots.

Tanshinone and phenolic acids are the most important active substances of Salvia miltiorrhiza, and the insight into their transcriptional regulatory mechanisms is an essential process to increase their content in vivo. SmMYB36 has been found to have important regulatory functions in the synthesis of tanshinone and phenolic acid; paradoxically, its mechanism of action in S. miltiorrhiza is not clear. Here, we demonstrated that SmMYB36 functions as a promoter of tanshinones accumulation and a suppressor of phenolic acids through the generation of SmMYB36 overexpressed and chimeric SmMYB36-SRDX (EAR repressive domain) repressor hairy roots in combination with transcriptomic-metabolomic analysis. SmMYB36 directly down-regulate the key enzyme gene of primary metabolism, SmGAPC, up-regulate the tanshinones biosynthesis branch genes SmDXS2, SmGGPPS1, SmCPS1 and down-regulate the phenolic acids biosynthesis branch enzyme gene, SmRAS. Meanwhile, SmERF6, a positive regulator of tanshinone synthesis activating SmCPS1, was up-regulated and SmERF115, a positive regulator of phenolic acid biosynthesis activating SmRAS, was down-regulated. Furthermore, the seven acidic amino acids at the C-terminus of SmMYB36 are required for both self-activating domain and activation of target gene expression. As a consequence, this study contributes to reveal the potential relevance of transcription factors synergistically regulating the biosynthesis of tanshinone and phenolic acid.

Leveraging interindividual variability of regulatory activity for refining genetic regulation of gene expression in schizophrenia

Schizophrenia is a polygenic psychiatric disorder with limited understanding about the mechanistic changes in gene expression regulation. To elucidate on this, we integrate interindividual variability of regulatory activity (ChIP-sequencing for H3K27ac histone mark) with gene expression and genotype data captured from the prefrontal cortex of 272 cases and controls. By measuring interindividual correlation among proximal chromatin peaks, we show that regulatory element activity is structured into 10,936 and 10,376 cis-regulatory domains in cases and controls, respectively. The schizophrenia-specific cis-regulatory domains are enriched for fetal-specific (p = 0.0014, OR = 1.52) and depleted of adult-specific regulatory activity (p = 3.04 × 10−50, OR = 0.57) and are enriched for SCZ heritability (p = 0.001). By studying the interplay among genetic variants, gene expression, and cis-regulatory domains, we ascertain that changes in coordinated regulatory activity tag alterations in gene expression levels (p = 3.43 × 10−5, OR = 1.65), unveil case-specific QTL effects, and identify regulatory machinery changes for genes affecting synaptic function and dendritic spine morphology in schizophrenia. Altogether, we show that accounting for coordinated regulatory activity provides a novel mechanistic approach to reduce the search space for unveiling genetically perturbed regulation of gene expression in schizophrenia.

In preprints: morphogens in motion.

In preprints: morphogens in motion.

Evolution and developmental expression of the sodium-iodide symporter (NIS, slc5a5) gene family: Implications for perchlorate toxicology.

The vertebrate sodium–iodide symporter (NIS or SLC5A5) transports iodide into the thyroid follicular cells that synthesize thyroid hormone. The SLC5A protein family includes transporters of vitamins, minerals, and nutrients. Disruption of SLC5A5 function by perchlorate, a pervasive environmental contaminant, leads to human pathologies, especially hypothyroidism. Perchlorate also disrupts the sexual development of model animals, including threespine stickleback (Gasterosteus aculeatus) and zebrafish (Danio rerio), but the mechanism of action is unknown. To test the hypothesis that SLC5A5 paralogs are expressed in tissues necessary for the development of reproductive organs, and therefore are plausible candidates to mediate the effects of perchlorate on sexual development, we first investigated the evolutionary history of Slc5a paralogs to better understand potential functional trajectories of the gene family. We identified two clades of slc5a paralogs with respect to an outgroup of sodium/choline cotransporters (slc5a7); these clades are the NIS clade of sodium/iodide and lactate cotransporters (slc5a5, slc5a6, slc5a8, slc5a8, and slc5a12) and the SGLT clade of sodium/glucose cotransporters (slc5a1, slc5a2, slc5a3, slc5a4, slc5a10, and slc5a11). We also characterized expression patterns of slc5a genes during development. Stickleback embryos and early larvae expressed NIS clade genes in connective tissue, cartilage, teeth, and thyroid. Stickleback males and females expressed slc5a5 and its paralogs in gonads. Single‐cell transcriptomics (scRNA‐seq) on zebrafish sex‐genotyped gonads revealed that NIS clade‐expressing cells included germ cells (slc5a5, slc5a6a, and slc5a6b) and gonadal soma cells (slc5a8l). These results are consistent with the hypothesis that perchlorate exerts its effects on sexual development by interacting with slc5a5 or its paralogs in reproductive tissues. These findings show novel expression domains of slc5 genes in stickleback and zebrafish, which suggest similar functions across vertebrates including humans, and provide candidates to mediate the effects of perchlorate on sexual development.

Mcrs1 is required for branchial arch and cranial cartilage development

Mcrs1 is a multifunctional protein that is critical for many cellular processes in a wide range of cell types. Previously, we showed that Mcrs1 binds to the Six1 transcription factor and reduces the ability of the Six1-Eya1 complex to upregulate transcription, and that Mcrs1 loss-of-function leads to the expansion of several neural plate genes, reduction of neural border and pre-placodal ectoderm (PPR) genes, and pleiotropic effects on various neural crest (NC) genes. Because the affected embryonic structures give rise to several of the cranial tissues affected in Branchio-otic/Branchio-oto-renal (BOR) syndrome, herein we tested whether these gene expression changes subsequently alter the development of the proximate precursors of BOR affected structures - the otic vesicles (OV) and branchial arches (BA). We found that Mcrs1 is required for the expression of several OV genes involved in inner ear formation, patterning and otic capsule cartilage formation. Mcrs1 knockdown also reduced the expression domains of many genes expressed in the larval BA, derived from either NC or PPR, except for emx2, which was expanded. Reduced Mcrs1 also diminished the length of the expression domain of tbx1 in BA1 and BA2 and interfered with cranial NC migration from the dorsal neural tube; this subsequently resulted in defects in the morphology of lower jaw cartilages derived from BA1 and BA2, including the infrarostral, Meckel's, and ceratohyal as well as the otic capsule. These results demonstrate that Mcrs1 plays an important role in processes that lead to the formation of craniofacial cartilages and its loss results in phenotypes consistent with reduced Six1 activity associated with BOR.

Highlighted articles for June 2022

Highlighted articles for June 2022