Computational Fluid Dynamic Simulations of Maternal Circulation: Wall Shear Stress in the Human Placenta and Its Biological Implications

Introduction

In the human placenta the maternal blood circulates in the intervillous space (IVS). The syncytiotrophoblast (STB) is in direct contact with maternal blood. The wall shear stress (WSS) exerted by the maternal blood flow on the STB has not been evaluated. Our objective was to determine the physiological WSS exerted on the surface of the STB during the third trimester of pregnancy.

Material and Methods

To gain insight into the shear stress levels that the STB is expected to experience in vivo, we have formulated three different computational models of varying levels of complexity that reflect different physical representations of the IVS. Computations of the flow fields in all models were performed using the CFD module of the finite element code COMSOL Multiphysics 4.4. The mean velocity of maternal blood in the IVS during the third trimester was measured in vivo with dynamic MRI (0.94±0.14 mm.s^-1). To investigate if the in silico results are consistent with physiological observations, we studied the cytoadhesion of human parasitized (Plasmodium falciparum) erythrocytes to primary human STB cultures, in flow conditions with different WSS values.

Results

The WSS applied to the STB is highly heterogeneous in the IVS. The estimated average values are relatively low (0.5±0.2 to 2.3±1.1 dyn.cm^-2). The increase of WSS from 0.15 to 5 dyn.cm^-2 was associated with a significant decrease of infected erythrocyte cytoadhesion. No cytoadhesion of infected erythrocytes was observed above 5 dyn.cm^-2 applied for one hour.

Conclusion

Our study provides for the first time a WSS estimation in the maternal placental circulation. In spite of high maternal blood flow rates, the average WSS applied at the surface of the chorionic villi is low (<5 dyn.cm^-2). These results provide the basis for future physiologically-relevant in vitro studies of the biological effects of WSS on the STB.     

Characterizaton of the Vessel Geometry,Flow Mechanics and Wall Shear Stress in the Great Arteries of Wildtype Prenatal Mouse

Introduction

Abnormal fluid mechanical environment in the pre-natal cardiovascular system is hypothesized to play a significant role in causing structural heart malformations. It is thus important to improve our understanding of the prenatal cardiovascular fluid mechanical environment at multiple developmental time-points and vascular morphologies. We present such a study on fetal great arteries on the wildtype mouse from embryonic day 14.5 (E14.5) to near-term (E18.5).

Methods

Ultrasound bio-microscopy (UBM) was used to measure blood velocity of the great arteries. Subsequently, specimens were cryo-embedded and sectioned using episcopic fluorescent image capture (EFIC) to obtain high-resolution 2D serial image stacks, which were used for 3D reconstructions and quantitative measurement of great artery and aortic arch dimensions. EFIC and UBM data were input into subject-specific computational fluid dynamics (CFD) for modeling hemodynamics.

Results

In normal mouse fetuses between E14.5–18.5, ultrasound imaging showed gradual but statistically significant increase in blood velocity in the aorta, pulmonary trunk (with the ductus arteriosus), and descending aorta. Measurement by EFIC imaging displayed a similar increase in cross sectional area of these vessels. However, CFD modeling showed great artery average wall shear stress and wall shear rate remain relatively constant with age and with vessel size, indicating that hemodynamic shear had a relative constancy over gestational period considered here.

Conclusion

Our EFIC-UBM-CFD method allowed reasonably detailed characterization of fetal mouse vascular geometry and fluid mechanics. Our results suggest that a homeostatic mechanism for restoring vascular wall shear magnitudes may exist during normal embryonic development. We speculate that this mechanism regulates the growth of the great vessels.     

Agrobacterium-mediated transformation of tomato elicits unexpected flower phenotypes with similar gene expression profiles

Background

Genetic transformation mediated by Agrobacterium tumefaciens is known to cause unexpected phenotypes. Mutations of a specific set of homeotic genes can result in alterred floral structure.

Methodology/Principal Findings

Previously we identified two genes (LeTGA1 and SOLly GLB1) induced by nutrient availability in tomato. To further elucidate their function, we sought to knock out the genes using antisense RNAi. When antisense constructs for the two different tomato genes were each transformed into Micro-Tina tomato plants, one primary transformant with similar mutant flower phenotypes was identified from transformation of each construct. Microarray analysis shows that a similar set of genes were up- or downregulated in both mutants. Sequencing of insertion sites indicates that each is inserted into a repetitive region which could impact expression of affected genes but direct alteration of floral homeotic gene sequences was not detected.

Conclusion

This is the first report that dominant flower mutations could be caused by genetic transformation designed to knock out two nutrient stress related genes.     

Self-Gated CINE MRI for Combined Contrast-Enhanced Imaging and Wall-Stiffness Measurements of Murine Aortic Atherosclerotic Lesions

Background

High-resolution contrast-enhanced imaging of the murine atherosclerotic vessel wall is difficult due to unpredictable flow artifacts, motion of the thin artery wall and problems with flow suppression in the presence of a circulating contrast agent.

Methods and Results

We applied a 2D-FLASH retrospective-gated CINE MRI method at 9.4T to characterize atherosclerotic plaques and vessel wall distensibility in the aortic arch of aged ApoE^−/− mice after injection of a contrast agent. The method enabled detection of contrast enhancement in atherosclerotic plaques in the aortic arch after I.V. injection of micelles and iron oxides resulting in reproducible plaque enhancement. Both contrast agents were taken up in the plaque, which was confirmed by histology. Additionally, the retrospective-gated CINE method provided images of the aortic wall throughout the cardiac cycle, from which the vessel wall distensibility could be calculated. Reduction in plaque size by statin treatment resulted in lower contrast enhancement and reduced wall stiffness.

Conclusions

The retrospective-gated CINE MRI provides a robust and simple way to detect and quantify contrast enhancement in atherosclerotic plaques in the aortic wall of ApoE^−/− mice. From the same scan, plaque-related changes in stiffness of the aortic wall can be determined. In this mouse model, a correlation between vessel wall stiffness and atherosclerotic lesions was found.     

A Novel Micro-Linear Vector for In Vitro and In Vivo Gene Delivery and Its Application for EBV Positive Tumors

Background

The gene delivery vector for DNA-based therapy should ensure its transfection efficiency and safety for clinical application. The Micro-Linear vector (MiLV) was developed to improve the limitations of traditional vectors such as viral vectors and plasmids.

Methods

The MiLV which contained only the gene expression cassette was amplified by polymerase chain reaction (PCR). Its cytotoxicity, transfection efficiency in vitro and in vivo, duration of expression, pro-inflammatory responses and potential application for Epstein-Barr virus (EBV) positive tumors were evaluated.

Results

Transfection efficiency for exogenous genes transferred by MiLV was at least comparable with or even greater than their corresponding plasmids in eukaryotic cell lines. MiLV elevated the expression and prolonged the duration of genes in vitro and in vivo when compared with that of the plasmid. The in vivo pro-inflammatory response of MiLV group was lower than that of the plasmid group. The MEKK1 gene transferred by MiLV significantly elevated the sensitivity of B95-8 cells and transplanted tumor to the treatment of Ganciclovir (GCV) and sodium butyrate (NaB).

Conclusions

The present study provides a safer, more efficient and stable MiLV gene delivery vector than plasmid. These advantages encourage further development and the preferential use of this novel vector type for clinical gene therapy studies.     

Evidence for selection on a CONSTANS-like gene between two red oak species

Background and Aims

Hybridizing species such as oaks may provide a model to study the role of selection in speciation with gene flow. Discrete species'' identities and different adaptations are maintained among closely related oak species despite recurrent gene flow. This is probably due to ecologically mediated selection at a few key genes or genomic regions. Neutrality tests can be applied to identify so-called outlier loci, which demonstrate locus-specific signatures of divergent selection and are candidate genes for further study.

Methods

Thirty-six genic microsatellite markers, some with putative functions in flowering time and drought tolerance, and eight non-genic microsatellite markers were screened in two population pairs (n = 160) of the interfertile species Quercus rubra and Q. ellipsoidalis, which are characterized by contrasting adaptations to drought. Putative outliers were then tested in additional population pairs from two different geographic regions (n = 159) to support further their potential role in adaptive divergence.

Key Results

A marker located in the coding sequence of a putative CONSTANS-like (COL) gene was repeatedly identified as under strong divergent selection across all three geographically disjunct population pairs. COL genes are involved in the photoperiodic control of growth and development and are implicated in the regulation of flowering time.

Conclusions

The location of the polymorphism in the Quercus COL gene and given the potential role of COL genes in adaptive divergence and reproductive isolation makes this a promising candidate speciation gene. Further investigation of the phenological characteristics of both species and flowering time pathway genes is suggested in order to elucidate the importance of phenology genes for the maintenance of species integrity. Next-generation sequencing in multiple population pairs in combination with high-density genetic linkage maps could reveal the genome-wide distribution of outlier genes and their potential role in reproductive isolation between these species.     

Cross-platform microarray data normalisation for regulatory network inference

Background

Inferring Gene Regulatory Networks (GRNs) from time course microarray data suffers from the dimensionality problem created by the short length of available time series compared to the large number of genes in the network. To overcome this, data integration from diverse sources is mandatory. Microarray data from different sources and platforms are publicly available, but integration is not straightforward, due to platform and experimental differences.

Methods

We analyse here different normalisation approaches for microarray data integration, in the context of reverse engineering of GRN quantitative models. We introduce two preprocessing approaches based on existing normalisation techniques and provide a comprehensive comparison of normalised datasets.

Conclusions

Results identify a method based on a combination of Loess normalisation and iterative K-means as best for time series normalisation for this problem.     

Low WSS Induces Intimal Thickening,while Large WSS Variation and Inflammation Induce Medial Thinning,in an Animal Model of Atherosclerosis

Objective

Atherosclerotic plaque development in the arterial wall is the result of complex interaction between the wall’s endothelial layer and blood hemodynamics. However, the interaction between hemodynamic parameters and inflammation in plaque evolution is not yet fully understood. The aim of the present study was to investigate the relation between wall shear stress (WSS) and vessel wall inflammation during atherosclerotic plaque development in a minipig model of carotid stenosis.

Methods

A surgical procedure was performed to create left common carotid artery stenosis by placement of a perivascular cuff in minipigs under atherogenic diet. Animals were followed up on 3T MRI, 1 week after surgery and 3, 6, and 8 months after initiation of the diet. Computational fluid dynamics simulation estimated WSS distribution for the first imaging point. Vascular geometries were co-registered for direct comparison of plaque development and features (Gadolinium- and USPIO-Contrast Enhanced MRI, for permeability and inflammation respectively) with the initial WSS. Histological analysis was performed and sections were matched to MR images, based on spatial landmarks.

Results

Vessel wall thickening, permeability and inflammation were observed distally from the stenosis. They were eccentric and facing regions of normal wall thickness. Histological analysis confirmed eccentric plaque formation with lipid infiltration, intimal thickening and medial degradation. High phagocytic activity in the stenosis region was co-localized with high WSS, corresponding to intense medial degradation observed on histology samples.

Conclusion

Lower WSS promotes atherosclerotic plaque development distal to an induced stenosis. Vascular and perivascular inflammation locations were predominant in the high WSS stenosis segment, where medial thinning was the major consequence.     

Combining Organophosphate Treated Wall Linings and Long-lasting Insecticidal Nets for Improved Control of Pyrethroid Resistant Anopheles gambiae

Background

New approaches to delivering insecticides need to be developed to improve malaria vector control. Insecticidal durable wall lining (DL) and net wall hangings (NWH) are novel alternatives to indoor residual spraying which can be produced in a long-lasting format. Non-pyrethroid versions could be used in combination with long-lasting insecticidal nets for improved control and management of insecticide resistant vector populations.

Methods

Experimental hut trials were carried out in Valley du Kou, Burkina Faso to evaluate the efficacy of pirimiphos methyl treated DL and NWH either alone or in combination with LLINs against pyrethroid resistant Anopheles gambiae ss. Comparison was made with pyrethroid DL. Mosquitoes were genotyped for kdr and ace-1R resistant genes to investigate the insecticide resistance management potential of the combination.

Results

The overall kdr and ace-1^R allele frequencies were 0.95 and 0.01 respectively. Mortality with p-methyl DL and NWH alone was higher than with pyrethroid DL alone (>95% vs 40%; P<0.001). Combining pyrethroid DL with LLINs did not show improvement in mortality (48%) compared to the LLIN alone (44%) (P>0.1). Combining p-methyl DL or NWH with LLINs reduced biting rates significantly (8–9%) compared to p-methyl DL and NWH alone (>40%) and killed all An gambiae that entered the huts. Mosquitoes bearing the ace-1^R gene were more likely to survive in huts with p-methyl DL alone (p<0.03) whereas all resistant and susceptible genotypes were killed by the combination.

Conclusion

P-methyl DL and NWH outperformed pyrethroid DL. Combining p-methyl DL and NWH with LLINs could provide significant epidemiological benefits against a vector population which is resistant to pyrethroids but susceptible to organophosphates. There was evidence that the single intervention would select kdr and ace-1^R resistance genes and the combination intervention might select less strongly. Technology to bind organophosphates to plastic wall lining would be worth developing.     

Assessing signatures of selection through variation in linkage disequilibrium between taurine and indicine cattle

Background

Signatures of selection are regions in the genome that have been preferentially increased in frequency and fixed in a population because of their functional importance in specific processes. These regions can be detected because of their lower genetic variability and specific regional linkage disequilibrium (LD) patterns.

Methods

By comparing the differences in regional LD variation between dairy and beef cattle types, and between indicine and taurine subspecies, we aim at finding signatures of selection for production and adaptation in cattle breeds. The VarLD method was applied to compare the LD variation in the autosomal genome between breeds, including Angus and Brown Swiss, representing taurine breeds, and Nelore and Gir, representing indicine breeds. Genomic regions containing the top 0.01 and 0.1 percentile of signals were characterized using the UMD3.1 Bos taurus genome assembly to identify genes in those regions and compared with previously reported selection signatures and regions with copy number variation.

Results

For all comparisons, the top 0.01 and 0.1 percentile included 26 and 165 signals and 17 and 125 genes, respectively, including TECRL, BT.23182 or FPPS, CAST, MYOM1, UVRAG and DNAJA1.

Conclusions

The VarLD method is a powerful tool to identify differences in linkage disequilibrium between cattle populations and putative signatures of selection with potential adaptive and productive importance.     

Differential Gene Repertoire in Mycobacterium ulcerans Identifies Candidate Genes for Patho-Adaptation

Background

Based on large genomic sequence polymorphisms, several haplotypes belonging to two major lineages of the human pathogen Mycobacterium ulcerans could be distinguished among patient isolates from various geographic origins. However, the biological relevance of insertional/deletional diversity is not understood.

Methodology

Using comparative genomics, we have investigated the genes located in regions of difference recently identified by DNA microarray based hybridisation analysis. The analysed regions of difference comprise ∼7% of the entire M. ulcerans genome.

Principal Findings

Several different mechanisms leading to loss of functional genes were identified, ranging from pseudogenization, caused by frame shift mutations or mobile genetic element interspersing, to large sequence polymorphisms. Four hot spot regions for genetic instability were unveiled. Altogether, 229 coding sequences were found to be differentially inactivated, constituting a repertoire of coding sequence variation in the rather monomorphic M. ulcerans.

Conclusions/Significance

The differential gene inactivation patterns associated with the M. ulcerans haplotypes identified candidate genes that may confer enhanced adaptation upon ablation of expression. A number of gene conversions confined to the classical lineage may contribute to particular virulence of this group comprising isolates from Africa and Australia. Identification of this spectrum of anti-virulence gene candidates expands our understanding of the pathogenicity and ecology of the emerging infectious disease Buruli ulcer.     

The Neovolcanic Axis Is a Barrier to Gene Flow among Aedes aegypti Populations in Mexico That Differ in Vector Competence for Dengue 2 Virus

Background

Aedes aegypti is the main mosquito vector of the four serotypes of dengue virus (DENV). Previous population genetic and vector competence studies have demonstrated substantial genetic structure and major differences in the ability to transmit dengue viruses in Ae. aegypti populations in Mexico.

Methodology/Principal Findings

Population genetic studies revealed that the intersection of the Neovolcanic axis (NVA) with the Gulf of Mexico coast in the state of Veracruz acts as a discrete barrier to gene flow among Ae. aegypti populations north and south of the NVA. The mosquito populations north and south of the NVA also differed in their vector competence (VC) for dengue serotype 2 virus (DENV2). The average VC rate for Ae. aegypti mosquitoes from populations from north of the NVA was 0.55; in contrast the average VC rate for mosquitoes from populations from south of the NVA was 0.20. Most of this variation was attributable to a midgut infection and escape barriers. In Ae. aegypti north of the NVA 21.5% failed to develop midgut infections and 30.3% of those with an infected midgut failed to develop a disseminated infection. In contrast, south of the NVA 45.2% failed to develop midgut infections and 62.8% of those with an infected midgut failed to develop a disseminated infection.

Conclusions

Barriers to gene flow in vector populations may also impact the frequency of genes that condition continuous and epidemiologically relevant traits such as vector competence. Further studies are warranted to determine why the NVA is a barrier to gene flow and to determine whether the differences in vector competence seen north and south of the NVA are stable and epidemiologically significant.     

Combined Effects of Flow Diverting Strategies and Parent Artery Curvature on Aneurysmal Hemodynamics: A CFD Study

Purpose

Flow diverters (FD) are increasingly being considered for treating large or giant wide-neck aneurysms. Clinical outcome is highly variable and depends on the type of aneurysm, the flow diverting device and treatment strategies. The objective of this study was to analyze the effect of different flow diverting strategies together with parent artery curvature variations on altering intra-aneurysmal hemodynamics.

Methods

Four ideal intracranial aneurysm models with different parent artery curvature were constructed. Computational fluid dynamics (CFD) simulations of the hemodynamics before and after applying five types of flow diverting strategies (single FD, single FD with 5% and 10% packing density of coils, two FDs with 25% and 50% overlapping rate) were performed. Changes in pressure, wall shear stress (WSS), relative residence time (RRT), inflow velocity and inflow volume rate were calculated and compared.

Results

Each flow diverting strategy resulted in enhancement of RRT and reduction of normalized mean WSS, inflow volume rate and inflow velocity in various levels. Among them, 50% overlapped FD induced most effective hemodynamic changes in RRT and inflow volume rate. The mean pressure only slightly decreased after treatment. Regardless of the kind of implantation of FD, the mean pressure, inflow volume rate and inflow velocity increased and the RRT decreased as the curvature of the parent artery increased.

Conclusions

Of all flow diverting strategies, overlapping FDs induced most favorable hemodynamic changes. Hemodynamics alterations post treatment were substantially influenced by parent artery curvature. Our results indicate the need of an individualized flow diverting strategy that is tailored for a specific aneurysm.     

A feature selection method for classification within functional genomics experiments based on the proportional overlapping score

Background

Microarray technology, as well as other functional genomics experiments, allow simultaneous measurements of thousands of genes within each sample. Both the prediction accuracy and interpretability of a classifier could be enhanced by performing the classification based only on selected discriminative genes. We propose a statistical method for selecting genes based on overlapping analysis of expression data across classes. This method results in a novel measure, called proportional overlapping score (POS), of a feature’s relevance to a classification task.

Results

We apply POS, along‐with four widely used gene selection methods, to several benchmark gene expression datasets. The experimental results of classification error rates computed using the Random Forest, k Nearest Neighbor and Support Vector Machine classifiers show that POS achieves a better performance.

Conclusions

A novel gene selection method, POS, is proposed. POS analyzes the expressions overlap across classes taking into account the proportions of overlapping samples. It robustly defines a mask for each gene that allows it to minimize the effect of expression outliers. The constructed masks along‐with a novel gene score are exploited to produce the selected subset of genes.

Electronic supplementary material

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

Relocation of genes generates non-conserved chromosomal segments in Fusarium graminearum that show distinct and co-regulated gene expression patterns

Background

Genome comparisons between closely related species often show non-conserved regions across chromosomes. Some of them are located in specific regions of chromosomes and some are even confined to one or more entire chromosomes. The origin and biological relevance of these non-conserved regions are still largely unknown. Here we used the genome of Fusarium graminearum to elucidate the significance of non-conserved regions.

Results

The genome of F. graminearum harbours thirteen non-conserved regions dispersed over all of the four chromosomes. Using RNA-Seq data from the mycelium of F. graminearum, we found weakly expressed regions on all of the four chromosomes that exactly matched with non-conserved regions. Comparison of gene expression between two different developmental stages (conidia and mycelium) showed that the expression of genes in conserved regions is stable, while gene expression in non-conserved regions is much more influenced by developmental stage. In addition, genes involved in the production of secondary metabolites and secreted proteins are enriched in non-conserved regions, suggesting that these regions could also be important for adaptations to new environments, including adaptation to new hosts. Finally, we found evidence that non-conserved regions are generated by sequestration of genes from multiple locations. Gene relocations may lead to clustering of genes with similar expression patterns or similar biological functions, which was clearly exemplified by the PKS2 gene cluster.

Conclusions

Our results showed that chromosomes can be functionally divided into conserved and non-conserved regions, and both could have specific and distinct roles in genome evolution and regulation of gene expression.

Electronic supplementary material

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