A simplistic approach of algal biofuels production from wastewater using a Hybrid Anaerobic Baffled Reactor and Photobioreactor (HABR-PBR) System

The current technologies of algal biofuels production and wastewater treatment (e.g., aerobic) process are still in question, due to the significant amount of fresh water and nutrients requirements for microalgae cultivation, and negative energy balance in both processes, especially when considered in the context of developing counties around the world. In this research, a simplistic sustainable approach of algal biofuels production from wastewater was proposed using a Hybrid Anaerobic Baffled Reactor (HABR) and Photobioreactor (PBR) system. The study suggests that the HABR was capable of removing most of the organic and solid (>90% COD and TSS removal) from wastewater, and produced a healthy feedstock (high N: P = 3:1) for microalgae cultivation in PBRs for biofuels production. A co-culture of Chlorella vulgaris, Chlorella sorokiniana, and Scenedesmus simris002 showed high lipid content up to 44.1%; and the dominant FAMEs composition (C16-C18) of 87.9% in produced biofuels. Perhaps, this proposed low-cost technological approach (e.g., HABR-PBR system) would connect the currently broken link of sustainable bioenergy generation and wastewater treatment pathway for developing countries.


Introduction
Currently, the world is facing a global sanitation crisis concerning wastewater management, about 70% of wastewater is treated in high-income countries, 38% in upper-middle-income, 28% in lower-middle-income, and only 8% in low-income countries [1]. On the contrary, most of these developing (low-income and lower-middle-income) counties are located either in a tropical or subtropical region with a warm climate (15-35˚C), which is favorable for biological wastewater treatment. Moreover, most of these countries also have electricity scarcity, which makes it challenging to promote aerobic wastewater treatment options [2]. At the same time, there is a global energy crisis due to the rapid utilization of fossil fuels in recent years for industrialization and urbanization. Renewable bioenergy is viewed as one of the ways to improve the current global energy crisis [3]. The energy demand and wastewater pollution PLOS ONE | https://doi.org/10.1371/journal.pone.0225458 December 5, 2019 1 / 24 a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 Although much research has been conducted on microalgae cultivation using raw wastewater or secondary/tertiary treated wastewater [8,9,[18][19][20][21][22][23], few studies have been carried out using anaerobic bioreactor effluent for microalgae cultivation which would be more cost-effective and sustainable for biofuels production and wastewater treatment. Therefore, this research work aims to determine the feasibility of biofuels production using effluent from the proposed HABR followed by microalgae cultivation in a photobioreactor (PBR) in the tropical climate of Bangladesh (e.g., developing country). The main objective of the research was to develop a low-cost, sustainable bioenergy generation corridor and to solve wastewater management/ treatment issues (i.e., a problem is a solution or resource) using this HABR-PBR system (Fig 1).
Also, the temperature has a significant effect on reactor treatment efficiency. Researchers [24][25][26] had shown that treatment efficiencies of the ABR deteriorated significantly when the temperature dropped below 15˚C. To overcome temperature effects, decreasing hydraulic retention times (HRT) or heating of wastewater could achieve higher removal efficiency [12], which also involves cost and energy. Therefore, the present study also examined the effect of the insulation of HABR (e.g., HABR (I)) in comparison with uninsulated HABR (e.g., HABR (U)) on effluent quality to cultivate microalgae in PBR for biofuels production.

Wastewater treatment by HABRs
HABR configuration and operation. The schematic diagram of the experimental set up with both uninsulated HABR (HABR (U)) and insulated HABR (HABR (I)) is shown in Fig 2  (A). Both reactors were identical in the configuration, as summarized in Table 1. These HABRs were constructed using the acrylic sheet with external dimensions of 90, 20, and 30 cm for length, width, and depth, respectively. The effective volume of HABR (U) and HABR (I) were 36.38 L and 36.39 L, respectively. Each HABR consisted of a front sedimentation chamber (U-1 and I-1), four regular chambers (U-2 to U-5, and I-2 to I-5) followed by two floated filter media chambers (U-6 and U-7; and I-6 and I-7). The first chamber volume, designed as the settling chamber, was twice the subsequent chambers. Each chamber was again divided into two portions by hanging baffles, which separated the chamber in down-and up-flow zone. The ratio between down-flow and up-flow was 1:4, and the bottom portion of the baffle was  media in the last two chambers of each reactor (Table 1). These locally available materials were used due to their favorable physical properties that would act as floated filter media, and won't let the reactor failure for clogging during wastewater treatment. Polyurethane foam (Pu Foam, Boya, Korea) was used for insulating one HABR by applying a 2-inch liquid foam layer and let it dry at room temperature (21-25˚C). Arduino UNIO technology (each chamber equipped with a DS18B20 waterproof digital temperature sensor connected to the data logger) was also installed to monitor temperature (24/7-1 h interval) during operation. Each compartment also had a 3-mm vent pipe (located behind temperature sensors pipe) to exhaust gas (e.g., methane). Both HABRs were operated under the same ambient conditions to evaluate the effluent quality for microalgae cultivation. Domestic wastewater was collected from KUET (Khulna University of Engineering & Technology, Khulna, Bangladesh) campus residential area, and stored in a feed tank system. The feed tank system consisted of 5 500-L premium grade PVC tanks interconnect each other, one tank (feeder) equipped with mixing device for uniform feed strength. The characteristics of influent wastewater are presented in Table 2. The wastewater was then fed from feeder tank to both HABRs continuously (running system 24/7-10 min feeding in each h) using a peristaltic pump (WT600-1F, Longer pump Co., China). The mixing device in the feeder tank and peristaltic pump were also connected to a Sino-timer (Sino timer, China), which was programmed to run the system for 10 min/hr during feeding time throughout the entire experiment. The HRTs of both HABR (U) and HABR (I) were 30 h for the first 40 days and then 20 h for remaining 100 days, except the system was put idle for 42 days (between 60 and 102 days). HABR inoculation. Both HABRs were inoculated with septic sludge collected from KUET campus residential area. The stable septic sludge was sieved using a 2.0-mm mesh prior to adding into reactors. In each HABR, approximately 9.2 L (3.2 L in chamber 1, and 1.5 L in each chamber 2-5) of sludge was added to chamber 1 to 5, the remaining volume of the chamber being filled with septic tank effluent including chambers 6 and 7. This seeded sludge contributed substantially to the solid requirement in the reactor system after settling. The sieved sludge contained total solids (TS) of 8960±1824 mg/L and total volatile solids (TVS) of 6880 ±1137 mg/L. After inoculating, both HABRs were left at ambient temperature for 30 d without further modification.
Wastewater sampling and analysis. Nine wastewater samples from each HABR (U & I): raw (U-R and I-R), seven sampling ports of each HABR (U-1 to U-7, and I-1 to I-7), and effluent (U-E and I-E), were collected routinely for laboratory analysis. Raw and effluent samples were analyzed for pH, oxygen redox potential (ORP), turbidity, total chemical oxygen demand (COD), ammonia-N (NH 4 , total suspended solids (TSS), volatile suspended solids (VSS), total coliform (TC) and fecal coliform (FC) according to the standard methods [27]. Samples collected from reactor chambers were also analyzed for selected parameters.
Microalgae cultivation, harvesting, lipid extraction, and biofuels conversation PBR configuration. The microalgae cultivation, harvesting, lipid extraction, and biofuels conversion are challenging steps and are on-going research topics throughout the world. In the present study, priority was given on cost-effectiveness (low-cost) in PBR design; hence, the gravity settling was considered for microalgae harvesting after cultivation. Four (4) identical 3.5 L capacity photobioreactors (PBR-1 to PBR-4) were constructed using locally available 2-L clear coke soft drink bottles attached back-to-back and were placed vertically with support for microalgal cultivation. The conical bottom shape of the PBR allowed the microalgae to settle at the bottom for harvesting (Fig 2(B)). A flow control valve was also installed at the lower end of the PBR for the collection of concentrated wet biomass after cultivation.
Microalgae strains selection and cultivation. To avoid the major risk of contamination of microalgae monocultures [28], a co-culture of Chlorella vulgaris, Chlorella sorokiniana, and Scenedesmus simris002 (ratio = 1:2:1) were used in the present study for microalgae cultivation. All three freshwater microalgae have been reported as having high lipid content (14-56.7%) by researches [29][30][31]; moreover, they also have good tolerance in saline water [31]. The research area (KUET campus) is located in a coastal zone with a high salinity in supplied water, which also contributes to the domestic wastewater.
Four (4) photobioreactors (PBR-1 to PBR-4) were used for microalgae cultivation. Two (2) two-headers (i.e. 2 x 3 L/min) capacity aquarium air pump (Venusaqua, China) were also installed for airflow to supply CO 2 in PBR during cultivation. Each air pump header was also connected by a 5-mm airflow tube to an air stone. The air stone was placed approximately 4 cm above the bottom of the reactor because of allowing gravity sedimentation of microalgae for harvesting. The airflow system was operated by the above mentioned Sino-timer (Sino timer, China), which was programmed to run the air pumps (e.g., CO 2 supply) for 10 min/hr (running continuously 24/7) during eight days of cultivation period. These PBRs were placed in different locations of a shed room (with the only roof) to observe the effect of sunlight irradiation on microalgal growth.
When comparing the effluent quality of both reactors regarding nutrients content, it appeared that HABR (U) effluent was higher in nutrient (high N and P) than HABR (I) ( Table 2). In addition, the total coliform (TC) and fecal coliform (FC) results showed that the higher TC/FC removal was achieved in HABR (U) (e.g. 75% TC and 59% FC removal in HABR (U)) ( Table 2), which reduced the risk of bacterial contamination during microalgae cultivation. Therefore, the effluent of HABR (U) was considered a healthy feedstock (high N: P = 3:1) [32,33] and was used as an influent of PBRs for microalgae cultivation in this study. Besides, selected wastewater parameters pH, DO, EC, ORP, and water temperature; and light irradiation (AM, noon, and PM) were also monitored for each day of cultivation period (8 days).
Microalgae growth kinetics and biomass estimation. The microalgal growth kinetics and biomass production potential were assessed for four PBRs, especially for the consideration of the light irradiation pattern on growth rate during cultivation. The optical density of microalgae cell culture suspensions was observed at 680 nm every day during eight days of cultivation using a HACH DR 3900 spectrophotometer. Biomass productivity (g/L/d) was calculated from the variation in biomass concentration (g/L) within a cultivation time (d) according to the following Eq (1) [34,35]: Where 'X t ' is the biomass concentration at the time 't', while 'X 0 ' is the initial biomass concentration at inoculation time t 0 . P max (g/L/d) was designated to the maximum productivity. The specific growth rate (μ, d -1 ) was calculated using Eq (2): The maximum specific growth rate (μ max , d -1 ) was determined from the different μ values calculated, while the maximum biomass obtained was designated as X max (g/L). Cell doubling time (t d (μ avg , d) was estimated using Eq (3) [34]: Carbon dioxide uptake rate (P CO 2 ) (g/L/d) was measured using Eq (4), whereas 'P' is the productivity calculated above in Eq (4) [34].
Microalgae harvesting and lipid extraction. The concentrated wet microalgae were harvested from the bottom of each PBR after cultivation using the flow control valve for analyses. The harvested algal biomass was collected together, and then split into six (6) samples: four (4) samples were used to determine moisture content by keeping samples at 65˚C in oven for overnight [36], and then two (2) (microalgae dry cell 1 and microalgae dry cell 2) out of these four (4) samples were taken for Fourier transforms infrared (FTIR) analysis; and two (2) wet microalgal samples were used to determine lipid content (LP1 and LP2), which were then converted to biofuels (biodiesel 1 and biodiesel 2) after the transesterification process [34,36,37].
As mentioned above, two (2) concentrated wet samples were used for lipid extraction using a single-step lipid extraction procedure [36]: 8 ml of a 2:1 chloroform: methanol (v/v) mixture was added to fresh microalgal biomass paste in each centrifuge tube. The biomass was manually suspended by vigorously shaking the tubes for a few seconds, and 2 ml of a 0.73% NaCl water solution was added. Phase separation was facilitated by 2 min of centrifugation at 350 g using NUVE NF 800 centrifuge, and the lower phase was recovered using a micropipette and was placed in an aluminum foil cup for overnight solvent evaporation at room temperature followed by gravimetrical determination of the lipid content (LP1 and LP2).
Transesterification process. Methylation of extracted lipids (LP1 and LP2) was conducted for converting all fatty acids to their corresponding methyl esters through the transesterification process, and the profile was then analyzed using FTIR and Gas Chromatography-Flame Ionization Detector (GC-FID) as the discussed in following sections. Also, a vegetable oil sample was also converted to biodiesel (e.g., biodiesel (veg. oil)) through the transesterification process and was analyzed using FTIR for comparison.
Fatty acid methyl esters (FAMEs) were prepared using 30 mg of total lipid dissolved in 1 ml of methanol, was then mixed with 1 ml of 1% NaOH solution prepared in methanol [38]. To this solution, equal volume (i.e., 2 ml) of 5% HCl (13.5 ml concentrated HCl in methanol) solution was added and then heated at 75˚C for 15 min. This solution was allowed to cool at room temperature, and 1 ml of distilled water was added and shaken. The organic layer containing FAMEs was carefully transferred to a new clean vial for FTIR and GC-FID analysis [39].
Fourier transform infrared spectroscopy (FTIR) analysis. The composition of microalgae dry cells and biofuels, and the type of functional groups of the algal dry cell and biofuels were assessed through FTIR spectroscopy study [40]. FTIR analyses were conducted on microalgae dry cell (1 & 2), biodiesel (veg. oil), and biodiesel (1 & 2) at room temperature using Shimadzu (IRTracer-100) FTIR spectrophotometer [34]. The dried algal biomass samples were further broken into powder. Dried algal cells were pressed against the diamond cell before scanning. The extracts from these samples were observed for their functionalities in the spectrogram. The spectra were collected in the mid-IR range from 4000 to 800 cm -1 (at a spectral resolution of 2 cm -1 ), and data were analyzed using Microsoft Excel, irAnalyze-RAMalyze (LabCognition GmbH & Co. KG) and KnowItAll (Bio-Rad Laboratories Inc., Pennsylvania, USA).
Gas chromatography-flame ionization detector (GC-FID) analysis. The FAMEs composition of the biodiesel sample (biodiesel 2) was analyzed using gas chromatography (GC-2010 Plus, Shimadzu, Japan) system with a flame ionization detector (FID) and equipped with a capillary column (MEGA SE-52 25 m × 0.25 mm × 0.25μm). 1.0 μL of methyl ester sample solution was injected for FAME analysis. Nitrogen was used as a carrier gas. The injector temperature was 180˚C, which was increased to 250˚C at a temperature gradient of 15˚C/min. The identification of FAMEs was made by comparing the retention times with those of the standard compounds [39].
Wastewater analysis after microalgae harvesting. After harvesting the concentrated wet microalgae from the PBRs, the remaining wastewater was analyzed for pH, ORP, Turbidity, COD, NO 3 --N, NO 2 --N, PO 4 3parameters.

Statistical analysis
Data analysis was performed with Microsoft Excel. The one-way ANOVA was used to determine the significance of the analytical results and difference between groups, and P<0.05 was considered as significant (S1 File).

Results and discussion
In the present study, pH and ORP were monitored for influent (U-R and I-R), effluent (U-E and I-E), and samples from each chamber of HABRs (U-1 to U-7, and I-1 to I-7) as presented in Table 2 (S1 and S2 Tables), and as shown in Fig 3(A) and 3(B). The pH was 8.1±0.2 for influent wastewater for both HABRs; and ORPs were 44.2±85.5 mV and 42.5±58.5 mV for influent wastewater of HABR (U) and HABR (I), respectively. The ORP ranges from -315 to 255 mV in the chambers of HABR(U), and from -137 to 117 mV in the chambers of HABR (I), respectively. This indicates a favorable anaerobic/anoxic/oxic condition existed in both HABRs for organics biodegradation as well as nitrification/denitrification/anammox processes. The temperature was also monitored using the Arduino UNIO temperature data logger as shown in Fig 3(C) (S2 File). It appeared that the insulation provided a better temperature control in HABR (I) during operation, whereas there was a significant temperature variation observed within the chambers of HABR (U) (e.g., mostly followed influent wastewater and ambient air temperature patterns). This variation of temperature within the reactor's chambers ultimately affected the treatment efficiency of the HABRs.

Wastewater treatment by HABRs
COD and solids removal. COD, turbidity, TSS, and VSS concentration of influent and effluent; and their removal efficiencies for both HABR (U) and HABR (I) are presented in Table 2. As the actual domestic wastewater was used for the experiments, the influent COD concentrations were observed to be varying [41]. The results indicated overall 93±7% and 89 ±9% COD removal by HABR (U) and HABR (I), respectively. It appeared that the COD removal efficiencies for both reactors fluctuated (e.g., 79-100% in HABR (U) and 86-100% in HABR (I)) during these experiments. The COD removal was 91±7% in HABR (U) and 86±8% in HABR (I) for 30 h HRT; and 94±7% in HABR (U) and 91±10% in HABR (I) for 20 h HRT, respectively, which indicates the COD removal efficiencies increase with lower HRT (e.g., 20 h). The organic loading rate (OLR) were 0.61±0.30 kgCOD/m 3 .d for HABR (U), and 0.67 ±0.31 kgCOD/m 3 .d for HABR (I), respectively. Also, the organic removal rate (ORR) were 0.57±0.29 kgCOD/m 3 .d for HABR (U) and 0.61±0.31 kgCOD/m 3 .d for HABR (I). The results indicate the COD removal is directly influenced by OLR [41,42]. No significant influence on COD removal efficiency was observed because of the insulation of the HABR. The average effluent COD was 39±37 mg/L for HABR (U), and 64±50 mg/ for HABR (I). As mentioned above, the ORP data indicated that there were favorable oxic/anoxic/anaerobic conditions existed in both the reactor's chambers for biological organic matter degradation [43].
During the experiments, turbidity was measured for influent and effluent samples for both reactors ( Table 2). The turbidity reduced significantly from 490±377 NTU to 11±7 NTU of in HABR (U), and from 563±343 NTU to 13±8 NTU in HABR (I); resulting 98±1% and 97±2% turbidity reduction in HABR (U) and HABR (I), respectively. Both HABRs showed superior performance for TSS and VSS removal: 91±9% for both TSS and VSS removal in HABR (U), and 96±5% of TSS and 95±9% of VSS in HABR (I). The effluent TSS was 15±5 mg/L and 20 ±17 mg/L in HABR (U) and HABR (I), respectively. Also, the effluent VSS was 9±3 mg/L and 12±8 mg/L in HABR (U) and HABR (I), respectively, suggests minimum biomass washout from both reactors. Feng et al. [24] have studied a bamboo carrier ABR and reported effluent TSS 14.35±3.01 mg/L (e.g., TSS removal of 81.92±3.53%) when operating at constant temperature 28±1˚C for 48 h HRT. In this study, both HABRs suggested higher TSS removal efficiency in comparison with their research. The average VSS/TSS ratio of influent wastewater was 0.55 for HABR (U), and 0.59 for HABR (I) suggested high VSS/TSS ratio for successfully anaerobic digestion [44]. A minor increase (e.g., 4-5%) of TSS or VSS removal was observed in HABR (I) because of insulation, which was perhaps due to ideal intermediate dispersion (e.g., weaker) existed in HABR (I) as suggested from a hydrodynamic study of the HABR [17].
Nitrogen and phosphate removal.  -N loading rate was also 0.05±0.02 kg/m 3 .d in both HABRs. The results suggested that at 30 h HRT, based on ORP (+49 to +96 mV in HABR (U), and +49 to +114 mV in HABR (I)), NH 4 + -N was primly removed by nitrification because of oxic/anoxic condition existed in both HABRs; however, at low 20 h HRT, NH 4 + -N was removed by both nitrification and perhaps anammox processes (-315 to +234 mV ORP in HABR (U), and -137 to +117 mV ORP in HABR (I)). Fig 4 also shows that the NO 3 --N removal rate followed its loading rate. An average NO 3 --N removal efficiency due to denitrification was 63±30% in HABR (U), and 75±26% in HABR (I). It appeared that NO 3 --N removal rate was not depended on NO 3 --N loading rate for both HABRs (not significant, p = 0.46>0.05 for HABR(U); p = 0.60>0.05 for HABR(I)). However, it was actually depended on NH 4 + -N loading rate (significant, p = 0.04<0.05 for HABR(U); p = 0.03<0.05 for HAB(I)). The influent ORP ranges (-80 to +230 mV in HABR (U), and -46 to +112 mV in HABR (I)) suggested that oxic/anoxic favorable condition for nitrification and denitrification [46]. However, these processes were not stable because of the significant variation of NH  Promising algal biofuels production and wastewater treatment: A simplistic approach for developing countries efficiency was 33±22% in HABR (U), and 38±24% in HABR (I in HABR (U)) because of biological phosphorus release by fermentative bacteria by producing fatty acids resulting higher phosphate concentration in the effluent. As a consequence, it decreased PO 4 3removal efficiency, which concerned with the finding reported by Schön et al. [47]. However, once HABR (I) recovered this situation and ran for an addition 100 d at 20 h HRT (except system idle for 42 days), higher PO 4 3removal 50±22% was achieved by HABR (I) compare to 40±23% in HABR (U). It indicated that the insulation of HABR also increased PO 4 3removal efficiency when operated for a longer time (140 d), which opposed to recent findings operating for a shorter time (e.g., 50 d) [45].

Microalgae cultivation, harvesting, lipid extraction, and biofuels conversion
Microalgae cultivation and biomass growth. Fig 5(A) shows the co-culture of Chlorella vulgaris, Chlorella sorokiniana, and Scenedesmus simris002, which was cultivated in four PBRs (PBR-1 to PBR-4) for eight days (S3 Table). During microalgae cultivation pH, DO, ORP and EC on a daily basis; and light irradiation and culture media temperature in morning (e.g., AM), noon, and afternoon (e.g., PM) were monitored as presented in Table 3. The average pH was 9.5±0.2 in PBR-1, 9.6±0.3 in PBR-2, 9.6±0.4 in PBR-3 and 9.7±0.4 in PBR during eight days of cultivation. The airflow to supply CO 2 in PBRs was 10 min/hr continuously during the entire cultivation period. The average pH of the effluent from HABR (U) was 8.0±0.2 raised in PBRs due to air sparging. As reported by Hu et al. [48] that high pH (8.5-10) resulted in intense chemical precipitation of PO 4 3-, and also decreased the bioavailability of inorganic carbon resulting lower NO 3 --N removal from wastewater. They suggested that the problem could be resolved by controlling pH or CO 2 sparging. In this study, it could be resolved by only air sparging during day time (during light irradiation) instead of 24/7. The DO, EC and ORP results suggested favorable culture media existed throughout the cultivation period. The light irradiation during the morning (e.g., AM), noon and afternoon (e.g., PM) were monitored as presented in Table 3. The microalgae concentration (g/L), productivity (g/L/d), specific growth rate (/d) during cultivation are presented in Table 4 and shown in Fig 5(B). The irradiation patterns on PBRs during cultivation are also shown in Fig 6 to observe the effect of light irradiation pattern on algal growth.
The results show that similar microalgal concentration in most PBRs except PBR-2 (e.g. low X max = 0.37 g/L). This was perhaps due to low light irradiation in PBR-2 (morning: 33 ±5 μmol/m 2 /s, noon: 23±7 μmol/m 2 /s, afternoon: 6±8 μmol/m 2 /s, respectively) compare to other PBRs (Table 4). The microalgal productivity (e.g. P max ) were observed high in PBR-3 (0.57 g/L/d) and PBR-4 (0.49 g/L/d) compare to PBR-1 (0.23 g/L/d) and PBR-2 (0.13 g/L/d). Both in PBR-1 and PBR-2, the light irradiation was high in the morning which then gradually decreased at noon and in the afternoon; however, in PBR-3 it was low in the morning then significantly increased at noon and decreased in the afternoon; and in PBR-4 it was mostly similar during the morning and noon but decreased in the afternoon.
The highest productivity 0.57 g/L/d in PBR-3 was perhaps due to this light irradiation pattern, which was low in the morning (110±31 μmol/m 2 /s) then gradually raised at noon (644 ±485 μmol/m 2 /s) and then dropped in the afternoon (38±63 μmol/m 2 /s). This suggested that the light irradiation pattern had a significant impact on microalgal growth. This indicates that the maximum microalgal growth could be achieved if the light irradiation pattern is similar to PBR-3 (e.g., low at the beginning, then gradually raises and then drops).
The results indicate the successful transesterification conversion of the microalgal dry cell lipid to biodiesel [34,40,51].

The overall performance of HABR-PBR system
The study suggests that both HABRs were capable of removing most of the organic (COD removal of 93±7% by HABR (U), and 89±9% by HABR (I)) and suspended solid (TSS removal of 91±9% by HABR(U), and 96±5% by HABR (I)) from domestic wastewater. Also, the effluent from HABR (U) contained higher nutrients (high N and P), and lower TC/FC (e.g., 75% TC and 59% FC removal in HABR (U)) than HABR (I), which reduced the risk of bacterial contamination during microalgae cultivation. Hence, the effluent of HABR (U) was considered a healthy feedstock (high N: P = 3:1) [32,33] and was used for microalgae cultivation.
A co-culture of Chlorella vulgaris, Chlorella sorokiniana, and Scenedesmus simris002 showed high lipid content up to 44.1%. The study also suggests that the sunlight irradiation pattern also has a significant influence on the productivity of algal biomass. The maximum microalgal growth could be achieved if the light irradiation pattern is similar to PBR-3 (e.g., low at the beginning then gradually raises and then drops). The results of the GC-FID analysis showed that the FAME profile was similar to the finding reported by several researchers [36,39,40,52]:  Promising algal biofuels production and wastewater treatment: A simplistic approach for developing countries dominant of palmitic acid (C16:0) of 20.5%, heptadecanoic acid (C17:1) of 8.6%, oleic acid (C18:1) of 12.1%, linoleic acid (C18:2) of 25.6%, and linolenic acid (C18:3) of 16.3%, which resulted 87.9% of the major FAME components (C16-C18). The quality of biodiesel would provide several advantages: quality ignition, low viscosity resulting lubricity, and higher oxidative stability for longer storage [39,40,[52][53][54].

Conclusion
In this research, a simplistic sustainable approach of algal biofuels production from wastewater was proposed using the HABR-PBR system. The study suggests that the HABR was capable of removing most of the organic and solid (>90% COD and TSS removal) from wastewater, and produced a healthy feedstock (high N: P = 3:1) for microalgae cultivation in PBRs for biofuels production. A co-culture of Chlorella vulgaris, Chlorella sorokiniana, and Scenedesmus sim-ris002 showed high lipid content up to 44.1%; and higher 87.9% of dominant FAMEs composition (C16-C18) in biodiesel.
The results suggest significant quality improvement of wastewater after microalgae cultivation and harvesting; however, there is still addition scopes to improve this HABR-PBR system, perhaps by reducing air sparging during cultivating, post flocculation of algal biomass before discharge. The study also suggests that there is a potential opportunity to re-use this treated water with further improvement to reduce water footprint for these developing countries. The HABR-PBR technological approach (i.e., a problem turned into a solution and resource) will work as a double-edged solution to mitigating wastewater problems and cogenerating algal biomass to produce bioenergy (as biofuels) to overcome energy demand especially for those countries located in a subtropical/tropical region with a warm climate.
Supporting information S1 Table.